Dear R-supporters, I have results of lme similar to those shown in Fig. 1.21, p.51 of Pinheiro/Bates. However, In my data set, one of the panels is missing, leading to an ugly frame shift of the following panels. How can I replace one of the panels by an empty one to restore the raster? Dieter Menne --------------------------------------- Dr. Dieter Menne Biomed Software 72074 T?bingen Tel (49) (7071) 52176 FAX (49) (7071) 55 10 46 dieter.menne at menne-biomed.de www.menne-biomed.de -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> Kjetil Halvorsen: > I dont have the bates&Pinheiro book here now, but I think there is a > skip= argument you can use to skip certain panels in the layout.Thanks, that gives a workaround, but it leaves a hole instead of an empty plot. I thought I was missing a standard mechanism that computes the position in the array from the lme result. Dieter -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._