使用 R 计算和 plot 每年回归的残差方差

Question

我用lm() function 计算了线性 model。 我知道如何使用 resid resid() function 提取残差，但我想按年份计算残差方差的值和 plot 它。

Answer 1

您可以使用残差和年份中的数据框，将其barplot aggreagte

fdf <- data.frame(resid=fit$residuals, year=fit$model$year)
res <- aggregate(resid ~ year, fdf, var)
barplot(resid ~ year, res)

---

玩具数据

set.seed(42)
dat <- transform(cbind.data.frame(year=1:5, x=rnorm(500)), y=.5 * x + rnorm(500))
fit <- lm(y ~ x + year, dat)