时间戳

首页 > 解决方案 > 在 R 中的计划比较中使用 summary() 和 summary.lm() - 为什么输出不同？

r - 在 R 中的计划比较中使用 summary() 和 summary.lm() - 为什么输出不同？

问题描述

我在 R 中运行了一个三向独立方差分析。（声音）操纵是我的自变量，具有三个级别：全等（KON）、不全等（INK）和无声音（控制）。此外，我已经构建了有计划的比较。第一个比较 c1 是 KON & 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

摘要.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() 输出中使用了我预定义的比较标签，即“控制与实验”和“INK 与 KON”，但它在第二个输出 summary.lm 中使用了其他东西()。为什么是这样？

此外，奇怪的是，第一次比较的 p 值在两个输出中不同，即在 summary() 的情况下为 0.0778，在 summary.lm() 的情况下为 0.0621。这种差异从何而来？

标签： rcomparisonsummaryanova

解决方案

你应该检查class(model)：

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

首先是“aov” - 所以summary(M)是一样的summary.aov(M)

推荐阅读