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

【问题标题】：Calculate and plot residuals variance of a regression, per year, using R使用 R 计算和绘制每年回归的残差方差
【发布时间】：2020-07-24 21:13:59
【问题描述】：

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

【问题讨论】：

【解决方案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)

【讨论】：