[R-sig-ME] hierarchical unbalanced design

Irene Mantzouni ima at difres.dk
Thu Nov 22 18:30:29 CET 2007


Hi all!
 
I have a mixed model with timeseries of y and x for populations:
lme(y~x,random=1|Population)
I am wondering if I can make it hierarchical by adding a random effect for year:
lme(y~x,random=1|Population/Year)
However the design is unbalanced: the time-series have different lengths for each population although the time-periods are overlapping in general. 
 
Would it be meaningful to include the Year effect this way? 
 
Thank you!
 
Irene




More information about the R-sig-mixed-models mailing list