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
______________________________________________
R-help@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.