I've shared a document with you: plot.png http://docs.google.com/leaf?id=0B_J5_yv8GrERYWRhZjU4NDItMmMyOS00ZGIxLWE2ZDUtMDFiYTY3MzliMzhk&hl=en&invite=CIa4nqkK
It's not an attachment -- it's stored online at Google Docs. To open this document, just click the link above. Hi there, Using the code below I generated the plot attached. The error bars at day = -3, 2 and 4 appear are larger. I was thinking if there is a way I could make all the error bars to be of the same size. I don't know if that makes sense. If not, then, is there a way I can plot only these 3 largest errorbars at day axis (day = -3, 2 and 4) and omit the rest of the error bars. I would be much obliged for any help Warmest Ogbos temp<-split(y,x) means<-sapply(temp,mean) stderr<-function(x)sqrt(var( x)/length(x)) p<-sapply(temp,stderr) n<-sapply(temp,length) ciw<-qt(0.975,n)*p/sqrt(n) library(gplots) day = -5:5 par(mfrow = c(3,1)) plotCI(day,means, ciw,gap = 0.05, sfrac= 0.005, type = "l",cex = 1.0) [[alternative HTML version deleted]] ______________________________________________ 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.