[R] lme plot

[Leonardo D Bacigalupe]

Hi all,

I have a questions regarding mixed effects models:

I'm trying to plot e.g. fitted vs residuals for each level of the 
random effects, and i'm getting the same error.
I guess this might be a problem of the graphic capabilities of R.

Is there any way to obtain those plots?

library(nlme)
attach(ergoStool)
names(ergoStool)
[1] "effort"  "Type"    "Subject"
m1<-lme(effort~Type,random=~1|Subject)
plot(m1,form=resid(.,type="p")~fitted(.)|Subject,abline=0)#resid and 
fitted for each level of subject
Error in as.vector(x, "list") : cannot coerce to vector


Thanks
leonardo


Leonardo D. Bacigalupe

Department of Animal and Plant Sciences
University of Sheffield
Sheffield S10 2TN
United Kingdom

L.Bacigalupe at sheffield.ac.uk