fixed (and this is where i get confused) how should
be the model:
(1) y~A,random=~1|A/B (in this option A is also random and should be
only fixed e.g. treatment)
(2) y~A,random=~A|B (does the A|B define the nested structure?)
thanks
leonardo
Leonardo D. Bacigalupe
Department of Animal and
3
63.82 4
65.82 4
56.63 1
57.53 1
77.83 2
79.23 2
69.93 3
69.23 3
62.13 4
64.53 4
Thanks
Leonardo
Leonardo D. Bacigalupe
Department of Animal and Plant Sciences
University of Sheffield
Sheffield S1
brary(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&