[R] mixed model nested ANOVA (part two)

S Ellison S.Ellison at lgc.co.uk
Mon Feb 25 00:04:19 CET 2008


>> Also i have read in Quinn and Keough 2002, design and analysis of
>> experiments for
>> biologists, that a variance component analysis should only be conducted
>> after a rejection
>> of the null hypothesis of no variance at that level.

Hmmm...

This does rather assume that 'no significant result' means 'near-zero variance contribution'.

These are not identical statements if the anova has low power; an 'insignificant' term can conceal practically important sizes of effect. So if you have a smallish number of groups (say, ten or less) you might want to find out what that estimated between-group variance could have been before you throw it away. That's especially important if you're expecting to say something about standard errors or confidence intervals of fixed effects. 

I may well be biased, here, though. In the kinds of nested design I get involved in (often inter-laboratory or homogeneity studies in chemistry), there is nearly always a between-group  effect; the only question is its size. Under those circumstances, the null hypothesis is not a particularly compelling starting point.  I'd rather have a variance component estimate and know how vague it was than assume it wasn't there at all.

But if you have good power and a good reason for believing there's no case to answer, sure; assume zero unless proven otherwise.

S


*******************************************************************
This email and any attachments are confidential. Any use...{{dropped:8}}



More information about the R-help mailing list