panel ordering in nlme and augPred plots

-- Mensaje Original --
From: "Petr Pikal" <petr.pikal at precheza.cz>
To: r-help at stat.math.ethz.ch
Date: Mon, 25 Apr 2005 16:40:40 +0200
Subject: [R] panel ordering in nlme and augPred plots


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.

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

Thank you
Best regards


Petr Pikal
petr.pikal at precheza.cz