每年使用R计算并绘制回归的残差方差

Question

我使用lm()函数计算了线性模型。我知道如何使用resid()函数提取残差，但是我想按年份计算残差方差的值并将其绘制出来。

Answer 1

您可以使用残差和年份中的数据框，aggreagte并制成barplot。

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)