Displaying 1 result from an estimated 1 matches for "box105bi".
Did you mean:
box105bis
2007 Dec 20
1
hierarchical linear models, mixed models and lme
...mmary(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 2 and 6 instead of 3*2=6 and 6* 3=18 !
I have also tried the following formula :
> box105bis.lme<-lme(content~treatment, box105.gd,
random=~1|treatment/rats/prep
In this case, the variance components as well as the number of rats and
preps are correct. Nevertheless, I have two new problems : (1) the
treatment is first treated as a random effects although it is fixed !
(2) there is...