アンバランスな分散分析では平方和の計算方法が複数ある。
普通のaov関数()では、タイプⅠの平方和しか算出できない。
タイプⅡ、Ⅲの平方和を求めるときは、library(car)のAnova関数を使う。
fit <- lm( y ~ x1 + x2 ) Anova(fit, type=2) Anova(fit, type=3)
fit <- lm( y ~ x1 + x2 )
Anova(fit, type=2)
Anova(fit, type=3)