[R-sig-ME] [R] lme nesting/interaction advice

Federico Calboli f.calboli at imperial.ac.uk
Mon May 12 11:20:47 CEST 2008


On 11 May 2008, at 22:45, Andrew Robinson wrote:
>>
>> lme(y ~ selection * males, random = ~1|replica/selection/males,  
>> mydata)
>
> forgive me, but I seem to see nesting in the random statement.   
> That is
> what happens when we separate factors with a '/'; they are nested.  We
> would expect that statement to not provide the correct df for the
> bog-standard fully crossed design.

Please read page 23/24 of the Pinheiro and Bates book, "Mixed-Effects  
Models in S and S-PLUS". It might prove enlightening.

F

--
Federico C. F. Calboli
Department of Epidemiology and Public Health
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG

Tel +44 (0)20 75941602   Fax +44 (0)20 75943193

f.calboli [.a.t] imperial.ac.uk
f.calboli [.a.t] gmail.com




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