[R] help- subtitles for multiple charts

[David Winsemius]

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