[R] hi friends, is there any wait function in R
Petr PIKAL
petr.pikal at precheza.cz
Wed Jul 15 14:39:55 CEST 2009
Hi
deepak m r <deepaknbr at gmail.com> napsal dne 14.07.2009 16:07:16:
> Hi,
> I am not an expert to debug the R can u please help.
> best regards
> deepak
See
fortune("brain")
You does not need to debug R. I believe R-Core will do it for you (If you
really find some bug in R, which seems to be highly improbable).
Instead see which plots are "wrong" and try to look at your data by any
suitable functions provided by R itself (str(), ls(), head().
In the first sight
1:length(paste("alp",s[t],"_mean",sep=""))
I expect to get 1. What do ***you*** expect to get?
Regards
Petr
>
> On Tue, Jul 14, 2009 at 10:06 AM, deepak m r<deepaknbr at gmail.com> wrote:
> > Hi,
> >
> >
> > On Tue, Jul 14, 2009 at 9:38 AM, Duncan Murdoch<murdoch at stats.uwo.ca>
wrote:
> >> On 7/14/2009 8:56 AM, deepak m r wrote:
> >>>
> >>> Hi,
> >>> Empty plot is getting i dont know why. can u please clarify how
> >>> can i use Print function instead of plot function.
> >>
> >> You need print() if you are using grid-based graphics (lattice,
ggplot2,...)
> >> in a script. You are using classic graphics so it should not be
necessary.
> >>
> >> If you are getting blank plots when you shouldn't, that's a bug. If
you can
> >> put together a reproducible example that shows its a bug in R, rather
than a
> >> bug in your script, it will likely be fixed fairly quickly.
> >>
> >> Duncan Murdoch
> >>
> >>
> >>> best regards
> >>> deepak
> >>>
> >>> On Tue, Jul 14, 2009 at 7:48 AM, Petr PIKAL<petr.pikal at precheza.cz>
wrote:
> >>>>
> >>>> Hi
> >>>>
> >>>> For this type of problems I do multipage pdf.
> >>>>
> >>>> pdf("file", ....)
> >>>> for (i in ...) {
> >>>> do all stuff including plot
> >>>> }
> >>>> dev.off()
> >>>>
> >>>> and then check the plots afterwards. Recently there was some post
about
> >>>> how to wait but you do not want only wait you want also to
interactively
> >>>> change plotting parameters, won't you.
> >>>>
> >>>> cat("\n","Enter x","\n") # prompt
> >>>> x=scan(n=1) # read 1 line from console
> >>>>
> >>>> this construction print something on console and reads one line
from
> >>>> console. There are also some packages which leave you choose from
several
> >>>> options. I think in car and randomForest are such routines so you
could
> >>>> check how it is done.
> >>>>
> >>>> Regards
> >>>> Petr
> >>>>
> >>>> r-help-bounces at r-project.org napsal dne 14.07.2009 13:17:01:
> >>>>
> >>>>> hi,
> >>>>> is there any wait function in R. I am running one R script to
plot
> >>>>> many graphs it is in the for loop. its showing no error but its
not
> >>>>> plotting well I think i can solve this problem with a wait
function.
> >>>>> Please help me in this regards. If u need any clarification about
> >>>>> programme. u can find the script below.
> >>>>>
> >>>>> best regards,
> >>>>> Deepak.M.R
> >>>>> Biocomputing Group
> >>>>> University of Bologana.
> >>>>>
> >>>>>
> >>>>> #!/usr/bin/R
> >>>>> s<-c
> >>>>>
> >>>>
> >>>>
>
("GG","GA","GV","GL","GI","GM","GF","GW","GP","GS","GT","GC","GY","GN","GQ","GD","GE","GK","GR","GH","AA","AV","AL","AI","AM","AF","AW","AP","AS","AT","AC","AY","AN","AQ","AD","AE","AK","AR","AH","VV","VL","VI","VM","VF","VW","VP","VS","VT","VC","VY","VN","VQ","VD","VE","VK","VR","VH","LL","LI","LM","LF","LW","LP","LS","LT","LC","LY","LN","LQ","LD","LE","LK","LR","LH","II","IM","IF","IW","IP","IS","IT","IC","IY","IN","IQ","ID","IE","IK","IR","IH","MM","MF","MW","MP","MS","MT","MC","MY","MN","MQ","MD","ME","MK","MR","MH","FF","FW","FP","FS","FT","FC","FY","FN","FQ","FD","FE","FK","FR","FH","WW","WP","WS","WT","WC","WY","WN","WQ","WD","WE","WK","WR","WH","PP","PS","PT","PC","PY","PN","PQ","PD","PE","PK","PR","PH","SS","ST","SC","SY","SN","SQ","SD","SE","SK","SR","SH","TT","TC","TY","TN","TQ","TD","TE","TK","TR","TH","CC","CY","CN","CQ","CD","CE","CK","CR","CH","YY","YN","YQ","YD","YE","YK","YR","YH","NN","NQ","ND","NE","NK","NR","NH","QQ","QD","QE","QK","QR","QH","DD","DE"!
> >>
> >>>>>
,"DK","DR","DH","EE","EK","ER","EH","KK","KR","KH","RR","RH","HH")
> >>>>> for(t in 1:length(s))
> >>>>> {
> >>>>>
> >>>>
> >>>>
a<-read.table(paste("../All_alpha_proteins/alp",s[t],"mean.sat",sep=""),header=T)
> >>>>>
> >>>>> attach (a)
> >>>>> names(a)
> >>>>> al<-1:length(paste("alp",s[t],"_mean",sep=""))
> >>>>>
> >>>>
> >>>>
b<-read.table(paste("../All_beta_proteins/bet",s[t],"mean.sat",sep=""),header=T)
> >>>>>
> >>>>> attach(b)
> >>>>> names(b)
> >>>>> bl<-1:length(paste("bet",s[t],"_mean",sep=""))
> >>>>> p<-read.table(paste("../Alpha_and_beta_proteins_a+b/apb",s
> >>>>> [t],"mean.sat",sep=""),header=T)
> >>>>> attach(p)
> >>>>> names(p)
> >>>>> pl<-1:length(paste("apb",s[t],"_mean",sep=""))
> >>>>> o<-read.table(paste("../Alpha_and_beta_proteins_aorb/aob",s
> >>>>> [t],"mean.sat",sep=""),header=T)
> >>>>> attach(o)
> >>>>> names(o)
> >>>>> ol<-1:length(paste("aob",s[t],"_mean",sep=""))
> >>>>> postscript(file=paste("Mean_",s[t],"_Plot.ps",sep=""))
> >>>>> plot(0,xlim=c(0,160),ylim=c(0,70),main=paste("Mean Distance Plot
for
> >>>>> ",s[t]," Rsidue pair",sep=""),ylab="Mean Distance in
> >>>>> Angstrom",xlab="Residue Seperation Number")
> >>>>> lines(al,paste("alp",s[t],"_mean",sep=""),col="blue",lty = 2)
> >>>>> lines(bl,paste("bet",s[t],"_mean",sep=""),col="yellow",lty = 2)
> >>>>> lines(pl,paste("apb",s[t],"_mean",sep=""),col="red",lty = 2)
> >>>>> lines(ol,paste("aob",s[t],"_mean",sep=""),col="green",lty = 2)
> >>>>> legend("topleft",c("Alpha Proteins","Beta Proteins","Alpha +
> >>>>> Beta","Alpha or Beta"),lty =
> >>>>> c(2,2,2,2),col=c("blue","yellow","red","green"))
> >>>>> dev.off()
> >>>>> }
> >>>>>
> >>>>> ______________________________________________
> >>>>> R-help at r-project.org mailing list
> >>>>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>>>> PLEASE do read the posting guide
> >>>>
> >>>> http://www.R-project.org/posting-guide.html
> >>>>>
> >>>>> and provide commented, minimal, self-contained, reproducible code.
> >>>>
> >>>>
> >>>
> >>> ______________________________________________
> >>> R-help at r-project.org mailing list
> >>> https://stat.ethz.ch/mailman/listinfo/r-help
> >>> PLEASE do read the posting guide
> >>> http://www.R-project.org/posting-guide.html
> >>> and provide commented, minimal, self-contained, reproducible code.
> >>
> >>
> >
More information about the R-help
mailing list