Displaying 1 result from an estimated 1 matches for "tryg".
Did you mean:
try
2003 Feb 13
1
fixed and random effects in lme
...o ask a question on fixed and random effecti in lme. I am
fiddlying around Mick Crawley dataset "rats" :
http://www.bio.ic.ac.uk/research/mjcraw/statcomp/data/
The advantage is that most work is already done in Crawley's book (page 361
onwards) so I can check what I am doing.
I am tryg to reproduce the nested analysis on page 368:
model<-aov(Glycogen~Treatment/Rat/Liver + Error(Treatment/Rat/Liver), rats)
using lme.
The code:
model1<- lme(Glycogen~Treatment, random = ~1|Rat/Liver, data=rats)
VarCorr(model1)
Variance StdDev
Rat = pdLogChol(1)...