汽车包装中不同的 III 型平方和

Question

我正在寻找 III 型平方和。我了解 Type III SS 默认情况下不同意 SAS。然而，我运行了以下两个模型并得到了不同的结果。我意识到我的缺乏理解可能在于我不理解 constr.poly() 语句。我在网上查了一下，似乎找不到直接的答案。感谢任何和所有的帮助！

install.packages("librarian")
librarian::shelf(car)

mtcars$cylC <- factor(mtcars$cyl)

model <- lm(mpg ~ cylC*disp, data=mtcars, contrasts=list(cylC=contr.poly))
car::Anova(model, type="III")
## Anova Table (Type III tests)
## 
## Response: mpg
##             Sum Sq Df  F value    Pr(>F)    
## (Intercept) 937.97  1 166.7687 8.101e-13 ***
## cylC        137.46  2  12.2201 0.0001814 ***
## disp         90.79  1  16.1415 0.0004465 ***
## cylC:disp    97.39  2   8.6574 0.0013133 ** 
## Residuals   146.23 26                       
## ---
## Signif. codes:  0 *** 0.001 ** 0.01 * 0.05 . 0.1   1


model <- lm(mpg ~ cylC*disp, data=mtcars)
car::Anova(model, contrasts=list(cylC=contr.poly), type="III")
## Anova Table (Type III tests)
## 
## Response: mpg
##              Sum Sq Df  F value    Pr(>F)    
## (Intercept) 1030.10  1 183.1477 2.791e-13 ***
## cylC         137.46  2  12.2201 0.0001814 ***
## disp         131.88  1  23.4472 5.096e-05 ***
## cylC:disp     97.39  2   8.6574 0.0013133 ** 
## Residuals    146.23 26                       
## ---
## Signif. codes:  0 *** 0.001 ** 0.01 * 0.05 . 0.1   1 

我在 Google 上搜索了 Type III SS 并查看了几页。我查看了 R 在线帮助（car、constr.poly、constrasts、Anova）。我查看了许多讨论 Type III SS 和 R 的页面，包括 Venables 的注释。

我希望两个 III 型 SS 都会同意。在这一点上，我显然没有从根本上理解。

Answer 1

在

?car::Anova

中说：

...请勿使用。

因此，您对

contrasts=list(cylC=contr.poly)

的调用中的

car::Anova()

将被忽略，且不会发出任何警告。

这是在同一模型上包含和排除该参数的示例。正如您所看到的，输出是相同的，这告诉我们包含参数没有区别。

mtcars$cylC <- factor(mtcars$cyl)

model <- lm(mpg ~ cylC*disp, data=mtcars)
car::Anova(model, type="III")
#> Anova Table (Type III tests)
#> 
#> Response: mpg
#>              Sum Sq Df  F value    Pr(>F)    
#> (Intercept) 1030.10  1 183.1477 2.791e-13 ***
#> cylC         137.46  2  12.2201 0.0001814 ***
#> disp         131.88  1  23.4472 5.096e-05 ***
#> cylC:disp     97.39  2   8.6574 0.0013133 ** 
#> Residuals    146.23 26                       
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

car::Anova(model, contrasts=list(cylC=contr.poly), type="III")
#> Anova Table (Type III tests)
#> 
#> Response: mpg
#>              Sum Sq Df  F value    Pr(>F)    
#> (Intercept) 1030.10  1 183.1477 2.791e-13 ***
#> cylC         137.46  2  12.2201 0.0001814 ***
#> disp         131.88  1  23.4472 5.096e-05 ***
#> cylC:disp     97.39  2   8.6574 0.0013133 ** 
#> Residuals    146.23 26                       
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

^{创建于 2024-04-22，使用 reprex v2.1.0}