R中线性模型的整体方差分析答案

【问题标题】：Overall ANOVA for Linear Models in RR中线性模型的整体方差分析
【发布时间】：2015-02-19 10:14:25
【问题描述】：

我可以使用anova 函数来获取线性模型的方差分析表。然而，这给出了每个解释变量对响应变量的单独影响。我想知道 R 中是否有任何函数可以获取overall ANOVA，我们可以在其中获取所有解释变量对响应变量的总体影响。在此先感谢您的帮助。

MWE

set.seed(12345)
X1 <- 10*abs(rnorm(n=10, mean = 5, sd = 1))
X2 <- 20*abs(rnorm(n=10, mean = 4, sd = 2))
Y <- 50*abs(rnorm(n=10, mean = 10, sd = 4))

df <- data.frame(X1, X2, Y)

fm1 <- lm(formula=Y~X1+X2, data=df)
anova(fm1)

Analysis of Variance Table

Response: Y
          Df Sum Sq Mean Sq F value Pr(>F)
X1         1  37842   37842  0.7815  0.406
X2         1 115774  115774  2.3910  0.166
Residuals  7 338941   48420

已编辑 (2017-10-23)

需要的输出是

Analysis of Variance Table

Response: Y
           Df Sum Sq      Mean Sq      F value  Pr(>F)
Regression 2  153616.587  76808.293    1.586    0.2703
Residuals  7  338940.570  48420.081

【问题讨论】：

你看过名字（fm1）吗？它列出了 anova 的所有部分，包括效果。
整体方差分析与拟合的 F 统计量有何不同，在 summary(fm1)?? 中给出？
@jlhoward：是的，整体方差分析将给出完全相同的拟合 F 统计量，由 summary(fm1) 给出。我想知道是否有任何R 函数来提供整体方差分析的每个组成部分。
@MYaseen208 如果你说 summary(fm1) 没有给你想要的东西，你的问题对我来说不是很清楚。
你找到这个问题的解决方案了吗？

标签： r statistics regression lm anova

【解决方案1】：

好的，让我们试试这个：拟合的 F 统计量将拟合与没有预测变量的拟合进行比较，例如适合Y~1。所以

anova(lm(Y~1,df),lm(Y~X1+X2,df))
# Analysis of Variance Table
#
# Model 1: Y ~ 1
# Model 2: Y ~ X1 + X2
#   Res.Df    RSS Df Sum of Sq      F Pr(>F)
# 1      9 492557                           
# 2      7 338941  2    153617 1.5863 0.2703

为您提供拟合的 F 统计量的组成部分：

summary(lm(Y~X1+X2,df))$fstatistic
#   value   numdf   dendf 
# 1.58629 2.00000 7.00000

【讨论】：