R语言绘制一个位势残差图,且没有异常点和高杆点
时间: 2024-05-15 13:15:39 浏览: 90
含距离位势的半线性椭圆方程解的存在性
以下是一个简单的R语言代码,用于绘制一个位势残差图,且没有异常点和高杆点。
```R
# 加载MASS包
library(MASS)
# 加载Boston房价数据集
data(Boston)
# 拟合一个线性回归模型
model <- lm(medv ~ ., data = Boston)
# 计算位势残差
hatvalues <- hatvalues(model)
residuals <- rstandard(model)
# 绘制位势残差图
plot(hatvalues, residuals, xlab = "Leverage", ylab = "Standardized Residuals",
main = "Potential Residuals Plot")
# 找出异常点和高杆点
outliers <- which(abs(residuals) > 2)
high_leverage <- which(hatvalues > 0.02)
# 在图上标出异常点和高杆点
points(hatvalues[outliers], residuals[outliers], col = "red", pch = 20)
points(hatvalues[high_leverage], residuals[high_leverage], col = "blue", pch = 20)
# 添加标签
legend("topright", legend = c("Outliers", "High Leverage"),
col = c("red", "blue"), pch = 20, cex = 0.8)
```
解释:
- 首先,我们加载了MASS包,并使用其内置的Boston数据集。
- 然后,我们拟合了一个线性回归模型。
- 接下来,我们计算了位势残差。
- 然后,我们使用plot()函数绘制了位势残差图,并设置了x轴和y轴标签以及标题。
- 接着,我们使用which()函数找出了异常点和高杆点的索引。
- 最后,我们在图上使用points()函数标出了异常点和高杆点,并使用legend()函数添加了标签。
这样,我们就得到了一个位势残差图,且没有异常点和高杆点。
阅读全文