Re: [R] panel ordering in nlme and augPred plots

From: Deepayan Sarkar <deepayan_at_stat.wisc.edu>
Date: Tue 26 Apr 2005 - 01:29:08 EST

On Monday 25 April 2005 09:40, Petr Pikal wrote:
> Dear all
>
> I am trying nlme together with Pinheiro/Bates book. I constructed
> grouped data object with suitable plotting layout (according to
> some common factor, panels from bottom to top are in increasing
> order).
>
> When I do nlme(... some stuff...) I get fitted object which I can
> plot with
>
> plot(augPred(fit.nlme6, level=0:1))
>
> but it results in completely different ordering. Is there any way
> how I can plot panels in some defined order e.g.

Could you give us a reproducible example? Following the example on the help page

fm <- lme(Orthodont)
plot(Orthodont)
plot(augPred(fm, level = 0:1))

gives me the same ordering on both plots.

> ord.f<-order(my.1.fac, my.2.fac)
> plot(augPred(fit.nlme6, level=0:1), ord.f)
>
> The only thing I found out is that if I order grouped.data object
>
> gr.dat<-gr.dat[ord.f,]
>
> and do nlme fit, then the ordering in augPred plot is OK.
>
> BTW
>
> is there a way how to specify in plot
>
> row1 -> 3 panels
> row2 -> 3 panels
> row3 -> 3 panels
> row4 -> 2 panels
> row5 -> 4 panels

Possibly. plot.augPred produces a Trellis plot, and usually arguments to the underlying plotting function can be passed on through the top-level call. e.g., with the Orthodont data

plot(augPred(fm1, level = 0:1), skip = rep(c(F,T), c(16, 2)))

or

p <- plot(augPred(fm1, level = 0:1))
update(p, skip = rep(c(F,T), c(16, 2)))

You would of course have to know what valid arguments are; for that see ?xyplot and ?update.trellis (in the lattice package).

Deepayan

---

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue Apr 26 01:33:47 2005