如果从残差分析图中获得异常值

残差对拟合值的那个图里的异常 就是残差绝对值最大的3个~
which(rank(-abs(fit$res))<4)

其他作图函数里标注的应该也是 头三个
用which就可以都找出来~

十里春风 发表于 2014-6-10 00:44
或者是绝对值大于两倍方差(标准差)？
which(abs(residuals(fit))>=2*var(residuals(fit))))

异常值outlier：一组测定值中与平均值的偏差超过两倍标准差的测定值。
与平均值的偏差超过三倍标准差的测定值，称为高度异常的异常值。which(abs(residuals(fit))>=2*sqrt(var(residuals(fit))))

事实上找异常值有个很常用的方法就是3σ准则
因为数值分布在（μ-3σ,μ+3σ)中的概率为0.9974
因此，超出此范围的数是异常值的概率就是为99.74%

，再来检查如何去除异常值。

异常观测值的判断：
远离回归线的离群点  ：学生化残差rstudent
远离自变量均值的杠杆点  ：杠杆值hatvalues
对回归线有重要影响的影响点。  ：考科斯距离cooks.distance
car包中influencePlot函数能将三个指标用气泡图形式将它们绘制出来，一个粗糙的判断准则：标准化残差值大于2或者小于-2的点可能是离群点。