我已绘制此图形以识别线性模型中的图形高杠杆点。
鉴于我的模型使用的数据集的变量“NOMBRES”,我试图绘制我的图形的所有点,但它变得难以辨认。这是我运行的代码:
> plot(hatvalues(tmodel),residuals(tmodel))
> text(hatvalues(tmodel),residuals(tmodel),labels=DSET$NOMBRES)
所以我想使用标签“DSET $ NOMBRES”来绘制杠杆(帽子值)高于0.05的点数。
答案 0 :(得分:5)
根据您的定义确定高杠杆点:
hlev <- which(hatvalues(tmodel)>0.05)
为图表添加数字标签:
text(hatvalues(tmodel)[hlev], residuals(tmodel)[hlev],
labels=DSET$NOMBRES[hlev])