The syntax below works, and gives the expected results, in R version 2.0.0,
provided that I have loaded packages Matrix, latticeExtra and lme4.
However, I cannot get it to work in version 2.2.1.   Can anyone tell me how
to fit this model in the more recent version?

 

Thanks in anticipation,

 

Nick Galwey

 

barleyprogeny.model1lme <- lme(yield_g_m2 ~ 1, 

   data = barleyprogeny.unbalanced, random = ~ 1|fline + fblock) 

 


	[[alternative HTML version deleted]]


