Message-ID: <200504251029.08685.deepayan@stat.wisc.edu>
Date: 2005-04-25T15:29:08Z
From: Deepayan Sarkar
Subject: panel ordering in nlme and augPred plots
In-Reply-To: <426D1D88.10548.1ECE11F@localhost>
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