使用 R 每年计算并绘制回归的残差方差
Calculate and plot residuals variance of a regression, per year, using R
我用 lm() 函数计算了一个线性模型。我知道如何使用 resid() 函数提取残差，但我想按年份计算残差方差值并绘制它。
您可以使用残差和年份中的数据框，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)
我用 lm() 函数计算了一个线性模型。我知道如何使用 resid() 函数提取残差，但我想按年份计算残差方差值并绘制它。
您可以使用残差和年份中的数据框，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)