R - 如何在 ANOVA 输出中显示空模型

Question

假设您在R 中拟合了一个带有截距项的模型。当您在模型上调用 anova 函数时，它不会在添加截距项之前显示空模型，甚至不会显示仅包含截距的模型。下面是一个使用简单线性回归的示例。

#Generate random dataset
set.seed(1)
n   <- 100
b0  <- 10
b1  <- 3
sig <- 4
xx  <- 10*runif(n)
ee  <- sig*rnorm(n)
DATA <- data.frame(x = xx, y = b0 + b1*xx + ee)

#Fit a linear model
MODEL <- lm(y ~ 1 + x, data = DATA)

#Show the ANOVA
anova(MODEL)

Analysis of Variance Table

Response: y
          Df Sum Sq Mean Sq F value    Pr(>F)    
x          1 6922.1  6922.1  488.52 < 2.2e-16 ***
Residuals 98 1388.6    14.2                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

我想为模型生成带有两个额外初步行的 ANOVA 表，从使用空模型的第一行开始没有截距项，然后显示使用模型的第二行只有拦截词。然后上面显示的剩余输出将从 ANOVA 表的第三行开始。

最简单的方法是什么？

Answer 1

MODEL1 <- lm(y ~ -1, data = DATA)
MODEL2 <- lm(y ~  1, data = DATA)
MODEL3 <- lm(y ~  x, data = DATA)
anova(MODEL1, MODEL2, MODEL3)