Re: [R] A model-building strategy in mixed-effects modelling

2010-01-19 Thread Ben Bolker
Stats Wolf stats.wolf at gmail.com writes: Dear all, Consider a completely randomized block design (let's use data(Oats) irrespoctive of the split-plot design it was arranged in). Look: library(nlme) fit - lme(yield ~ nitro, Oats, random = ~1|Block, method=ML) fit2 - lm(yield ~ nitro

[R] A model-building strategy in mixed-effects modelling

2010-01-18 Thread Stats Wolf
Dear all, Consider a completely randomized block design (let's use data(Oats) irrespoctive of the split-plot design it was arranged in). Look: library(nlme) fit - lme(yield ~ nitro, Oats, random = ~1|Block, method=ML) fit2 - lm(yield ~ nitro + Block, Oats) anova(fit, fit2) gives this: Model