撇开关于第三类方差分析和边缘原则的辩论以及所有......
我已经设置了两个模型,其平方和应该是不同的(并且III型ANOVA将测试该差异)。这是代码:
library(car)
library(openintro)
data(hsb2)
hsb2$gender <- factor(hsb2$gender)
contrasts(hsb2$gender) <- "contr.sum"
contrasts(hsb2$ses) <- "contr.sum"
math_gender_int <- lm(math ~ gender + gender:ses, data = hsb2)
math_gender_ses_int <- lm(math ~ gender + ses + gender:ses, data = hsb2)
现在我应该能够看到这两个模型之间的平方和的差异。毕竟,&#34;完整&#34;模型还有一个术语:
anova(math_gender_int, math_gender_ses_int)
但输出结果显示:
Analysis of Variance Table
Model 1: math ~ gender + gender:ses
Model 2: math ~ gender + ses + gender:ses
Res.Df RSS Df Sum of Sq F Pr(>F)
1 194 15858
2 194 15858 0 -1.819e-12
这里发生了什么?