Question

我在R中运行了一个三向独立的ANOVA。（声音）操纵是我的三个级别的自变量：全等（KON），不一致（INK）和没有声音（控制）。此外，我已经构建了计划比较。第一个比较c1是KON＆amp;的对比度。 INK与对照组和第二比较c2是KON与INK的对比。输出看起来像这样：

摘要（模型）

                                         Df Sum Sq Mean Sq F value Pr(>F)
Manipulation                              2  11.97   5.985 2.388 0.0975 .
  Manipulation: control vs. Experimental  1   7.97   7.970 3.181 0.0778 .
  Manipulation: INK vs. KON               1   4.00   3.999 1.596 0.2097 
Residuals                                91 228.01   2.506                                           
---
Signif. codes:  
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

summary.lm（模型）

Residuals:
    Min      1Q  Median      3Q     Max 
-2.5062 -1.3333 -0.3333  1.1398  4.4111 

Coefficients:
               Estimate Std. Error t value Pr(>|t|)    
(Intercept)      3.0317     0.1647  18.411   <2e-16 *** 
Manipulationc1  -0.2214     0.1172  -1.889   0.0621 .  
Manipulationc2  -0.2531     0.2003  -1.263   0.2097    
---
Signif. codes:  
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 1.583 on 91 degrees of freedom
Multiple R-squared:  0.04988,   Adjusted R-squared:  0.02899 
F-statistic: 2.388 on 2 and 91 DF,  p-value: 0.0975

让我印象深刻的是R使用我预先定义的比较标签，即第一个summary（）输出中的“control vs. experimental”和“INK vs. KON”，但它在第二个输出中使用了其他内容summary.lm（）。这是为什么？

此外，似乎奇怪的是，第一次比较的p值在两个输出中不同，即在summary（）的情况下为0.0778，在summary.lm（）的情况下为0.0621。这种差异来自哪里？

Answer 1

您应该检查class(model)：

M <- aov(formula = Petal.Length ~ Species, data = iris)
summary(M)
summary.lm(M)
class(M)

首先是＆＃34; aov＆＃34; - 所以summary(M)与summary.aov(M)

相同

在R中的计划比较中使用summary（）和summary.lm（） - 为什么输出会有所不同？

1 个答案: