残差与杠杆关系图(使用R语言)
残差与杠杆关系图是统计学中常用的可视化工具,用于识别在线性回归模型中具有异常杠杆或较大残差的数据点。本文将介绍如何使用R语言绘制残差与杠杆关系图,并提供相应的源代码。
在R中,我们可以使用lm()
函数来拟合线性回归模型,并使用plot()
函数来绘制残差与杠杆关系图。下面是一个示例代码,演示了如何进行这些操作:
# 载入所需的包
library(ggplot2)
# 创建示例数据
x <- c(1, 2, 3, 4, 5)
y <- c(2, 4, 6, 8, 10)
# 拟合线性回归模型
model <- lm(y ~ x)
# 提取残差和杠杆值
residuals <- residuals(model)
leverage <- hatvalues(model)
# 创建残差与杠杆关系图
plot(residuals ~ leverage, xlab = "杠杆", ylab = "残差", main = "残差与杠杆关系图")
# 添加标识异常值的标签
abline(h = 0, col = "red", lty = 2)
text(x = leverage, y = residuals, labels = ifelse(abs(residuals) > 2 * sd(residua