Re: [R] lme nesting/interaction advice

From: Federico Calboli <>
Date: Mon, 12 May 2008 10:20:47 +0100

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.


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]
f.calboli [.a.t]

______________________________________________ mailing list
PLEASE do read the posting guide
and provide commented, minimal, self-contained, reproducible code.
Received on Mon 12 May 2008 - 09:28:20 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 12 May 2008 - 09:30:37 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive