lme nesting/interaction advice
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