Hello All: Can anybody tell me what is the problem with my program please. I have an error message as appears below. My program is: ifn <- "Jul08_09.LM" data <- read.table(ifn) ofn <- "Jul.png" bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg "white",res=50) par(mar=c(5, 5, 3, 2),lwd=5) par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) par(mfrow = c(3,4)) ifn <- "T.dat" trait <- read.table(ifn) i <- 1 j <- 1 for(k in 1:12) { dat <- data[i:(i+2), ] colnames(dat)<-c("SM","E","NP) Ymin <- if( min(dat$SM) < 0.0 ) (dat$SM - dat$E) else 0.0 Ymax <- if( max(dat$SM) > 0.0 ) (dat$SM + dat$E) else 0.0 Graph<-barplot(dat$SM, names.arg=dat$NP, main = trait[j:j,], xlab data[i:i,1:1], ylim = c(Ymin,Ymax) ) segments(Graph, dat$SM + dat$E, Graph, dat$SM - dat$E) i <- i + 3 j <- j + 1 } dev.off() And the error message is: Error in plot.window(xlim, ylim, log = log, ...) : invalid 'ylim' value Calls: barplot -> barplot.default -> plot.window Execution halted Thanks very much, Mohsen [[alternative HTML version deleted]]
Hi, On Aug 6, 2009, at 4:48 PM, Mohsen Jafarikia wrote:> Hello All: > Can anybody tell me what is the problem with my program please. I > have an > error message as appears below. > > My program is: > > ifn <- "Jul08_09.LM" > data <- read.table(ifn) > > ofn <- "Jul.png" > > bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = > 30, bg > "white",res=50) > par(mar=c(5, 5, 3, 2),lwd=5) > par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) > > par(mfrow = c(3,4)) > > ifn <- "T.dat" > trait <- read.table(ifn) > > i <- 1 > j <- 1 > > for(k in 1:12) > { > dat <- data[i:(i+2), ] > > colnames(dat)<-c("SM","E","NP) > > Ymin <- if( min(dat$SM) < 0.0 ) (dat$SM - dat$E) else 0.0 > Ymax <- if( max(dat$SM) > 0.0 ) (dat$SM + dat$E) else 0.0 > > Graph<-barplot(dat$SM, names.arg=dat$NP, main = trait[j:j,], xlab > data[i:i,1:1], ylim = c(Ymin,Ymax) ) > > segments(Graph, dat$SM + dat$E, Graph, dat$SM - dat$E) > i <- i + 3 > j <- j + 1 > } > > dev.off() > > And the error message is: > > Error in plot.window(xlim, ylim, log = log, ...) : invalid 'ylim' > value > Calls: barplot -> barplot.default -> plot.window > Execution haltedDid you take a look at what c(Ymin, Ymax) is? -steve -- Steve Lianoglou Graduate Student: Computational Systems Biology | Memorial Sloan-Kettering Cancer Center | Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact
Mohsen Jafarikia wrote:> Hello All: > Can anybody tell me what is the problem with my program please. I have an > error message as appears below. > > My program is: > > ifn <- "Jul08_09.LM" > data <- read.table(ifn) > > ofn <- "Jul.png" > > bitmap(ofn, type = "png256", width = 30, height = 30, pointsize = 30, bg > "white",res=50) > par(mar=c(5, 5, 3, 2),lwd=5) > par(cex.main=1.6,cex.lab=1.6,cex.axis=1.6) > > par(mfrow = c(3,4)) > > ifn <- "T.dat" > trait <- read.table(ifn) > > i <- 1 > j <- 1 > > for(k in 1:12) > { > dat <- data[i:(i+2), ] > > colnames(dat)<-c("SM","E","NP) > > Ymin <- if( min(dat$SM) < 0.0 ) (dat$SM - dat$E) else 0.0 > Ymax <- if( max(dat$SM) > 0.0 ) (dat$SM + dat$E) else 0.0 > > Graph<-barplot(dat$SM, names.arg=dat$NP, main = trait[j:j,], xlab > data[i:i,1:1], ylim = c(Ymin,Ymax) ) > > segments(Graph, dat$SM + dat$E, Graph, dat$SM - dat$E) > i <- i + 3 > j <- j + 1 > } > > dev.off() > > And the error message is: > > Error in plot.window(xlim, ylim, log = log, ...) : invalid 'ylim' value > Calls: barplot -> barplot.default -> plot.window > Execution halted > > >Hi Mohsen, Try this: Ymin<-ifelse(min(dat$SM)<0.0,dat$SM-dat$SE,0) Ymax<-ifelse(min(dat$SM)>0.0,dat$SM+dat$SE,0) cat(Ymin,Ymax,"\n") Jim