Displaying 1 result from an estimated 1 matches for "box105".
Did you mean:
box100
2007 Dec 20
1
hierarchical linear models, mixed models and lme
...able ; 2 readings of the glycogen content / preparations
are realised. Treatment is fixed whereas Rats (nested in Treatment) and
Prep (nested in Rats) are random effects.
According to a previous discussion found in the R-help archives (January
2007), I have tried the following formula :
> box105.lme<-lme(content~treatment, box105.gd, random=~1|rats/prep)
However, the formula summary(box105.lme) gives wrong estimates for the
variance components (rats within treatments, prep within rats, readings
within preps) ! Moreover the numbers of rats and preps are also wrong,
with respectively...