Skip to content
Prev 14255 / 20628 Next

Replicating type III anova tests for glmer/GLMM

On Tue, Feb 23, 2016 at 01:06:18PM +0100, Francesco Romano wrote:
? On the other hand, I don't understand how Cai et al. (2012) p.842,
? "combined analysis experiments 1 and 2", reported the main effect
? of a factor with 4 levels via a single estimate, SE, z, p coefficient.
? How did they obtain this and is this the right way?

It's just a guess, but any sum-of-square can be seen as a particular
contrast, that is a particular combination of the coefficients in the
model (or of the different means, expressed another way) that is
tested against 0. So I guess this single estimate is the value of the
contrast associated to the corresponding sum-of-squares, and SE/z/p
are derived similarly.

You can play with multcomp::glht to test this, but knowing which
contrast is tested by which sum of square in a specific desing may be
tricky: it depends on the coding, on the (un)balance...

Kowing if this is the ? right ? way is I think the same debate that
knowing which kind of sum-of-square should be used and the question is
very application dependent. Just, if you don't know what this single
estimate estimates really, interpretation is at best difficult...