[R] help- subtitles for multiple charts

[David Winsemius]

On Jun 23, 2011, at 10:40 PM, David Winsemius wrote:

>
> 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*):
>>
> To get the "Block and Trt text as in the original request (which you  
> didn't copy for reference) you would use bquote:
>
>> 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; curlev <- strsplit(as.character(i), ".",  
>> fixed=TRUE)
# also add           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

>> plot(pMBE$Year[ss], pMBE$DBH[ss], xlab="Year", ylab="DBH [in]",
>> ylim=range(0:12), col=28, sub= bquote(Block~.(curlev[[1]] 
>> [[1]])~Trt~.(curlev[[1]][[2]])) )
> # 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 at 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.

David Winsemius, MD
West Hartford, CT