On Jun 23, 2011, at 4:41 PM, jalen wrote:
http://r.789695.n4.nabble.com/file/n3620982/mbeFORUM.csv mbeFORUM.csv I uploaded my data and one more time the script (*adjusted version*):
Try:
pMBE<- MBE[MBE$left!=0,] pMBE$bt<- interaction(pMBE$Block,pMBE$trt) par(mfrow=c(3,3), oma=c(2,0,2,0)) for(i in unique(pMBE$bt)){ ss <- pMBE$bt==i plot(pMBE$Year[ss], pMBE$DBH[ss], xlab="Year", ylab="DBH [in]", ylim=range(0:12), col=28, sub=as.character(i) )
# add ^^^^^^^^^^^^^^^^^^
sm<-loess(DBH~Year, data=pMBE[ss,]) x=seq(2004, 2010, by=1) points(x, predict(sm, data.frame("Year"=x)), type="l", col=650, title("Diameter", adj=0.5, outer=T)) }
David Winsemius, MD West Hartford, CT ______________________________________________ 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.