[R] Missing panels in multipanel lme lattice/trellis

Dieter Menne dieter.menne at menne-biomed.de
Fri Nov 23 08:55:34 CET 2001


Paul,

> I don't have Pinheiro/Bates handy;  would you be able to send me the code
> that produces the output so I can take a look?

Here is a the example, with a simulated missing block/variety.
Following Kjetil's suggestion, I manually added a Skip-list.
Maybe a more general way to create the skip-list would be good 
enough, but possible one needs a dummy-panel that is inserted
for missing matrix-combinations.

# Trellis with missing graphs
library(nlme)
library(lattice)
data(Oats)
# Remove one set to simulate missing data
Oats1<-Oats[(Oats$Block!="III") | (Oats$Variety!="Victory"),]
# attempted correction, stupid, but good enough for ONE plot
SkipIt<-rep(F,18)
SkipIt[7]<-T # just "manual"
fm4Oats<-lme( yield ~ nitro, data = Oats1, random = ~ 1 | Block/Variety )
# produces Figure 1.21 of Pinheiro/Bates, 
# with one missing. 
plot( augPred( fm4Oats ), aspect = 2.5, layout = c(6, 3),
      between = list( x = c(0, 0, 0.5) ),skip=SkipIt ) 
# The plot is not really nice, I would prefer an empty
# box instead of the hole

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list