残差诊断

残差图诊断

par(mfrow = c(2,2))
# 普通残差图
plot(fitted(model), residuals(model), xlab="拟合值", ylab="普通残差",pch=19)  
abline(h=0, lty=2)
# 删除残差图
plot(fitted(model), residuals(model)/(1-hatvalues(model)), xlab="拟合值", ylab="删除残差", pch=19)   
abline(h=0,lty=2)
# 学生化残差图
plot(fitted(model), rstandard(model), xlab="拟合值", ylab="学生化残差", ylim=c(-3.5,3.5), pch=19)  
abline(h=0, lty=2)
# 删除学生化残差图
plot(fitted(model), rstudent(model), xlab="拟合值", ylab="删除学生化残差", ylim=c(-3.5,3.5), pch=19)   
abline(h=0, lty=2)

使用Plot做回归诊断图

• 有无非线性关系

  plot(model,which=1)

• 残差是否正态分布

  plot(model, which=2)

• 是否有异方差

  plot(model, which=3)

• 查看强影响点

  plot(model, which=4)