检查III型ANOVA结果

时间:2017-03-31 17:20:35

标签: r lm anova

撇开关于第三类方差分析和边缘原则的辩论以及所有......

我已经设置了两个模型,其平方和应该是不同的(并且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      

这里发生了什么?

0 个答案:

没有答案