2014-02-11 117 views
5

我已经绘制了此图形来识别线性模型中图形化的高杠杆点。 enter image description here在散点图中标出离群点

鉴于我的模型使用的数据集的变量“NOMBRES”,我试图绘制我的图形的所有点,但难以辨认。下面是我跑的代码:

> plot(hatvalues(tmodel),residuals(tmodel)) 
> text(hatvalues(tmodel),residuals(tmodel),labels=DSET$NOMBRES) 

所以我想使用标签“DSET $ NOMBRES”上述0.05与杠杆(帽子值)来绘制只是点。

+3

'离群值< - 其中[hatvalues(tModel的)> 0.05]; text(hatvalues(tmodel)[outliers],residuals(tmodel)[outliers],labels = DSET $ NOMBRES [outliers])' –

+0

@BenBolker作为回答? – Thomas

回答

5

按照自己的定义确定高杠杆点:

hlev <- which(hatvalues(tmodel)>0.05) 

数字标签添加到所述图形:

text(hatvalues(tmodel)[hlev], residuals(tmodel)[hlev], 
    labels=DSET$NOMBRES[hlev])