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

2010-01-19 Thread Ben Bolker
Stats Wolf 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: Mo