黄土线没有正确绘图

Question

我无法在一个残差的散点图上拟合黄土平滑加信心限制.

我的模型是身高〜体重+胸围.为了检查胸围的线性,我拟合了一个没有胸围的模型(即身高〜体重),并绘制了该模型对胸围的残差.到现在为止还挺好.然后我尝试使用loess()和predict()绘制黄土线,加上置信限制.结果看起来像这样(在图片中我只绘制了中心线,但CI线看起来相同):

散点图中的黄土问题

这些点是正确的(当我将黄土拟合为它看起来正确的点时),但由于某种原因,线条没有被绘制出我的预期.我的代码如下:

# bf.red = data set; mod.nch = model; chestc = chest circumference;
# loess = loess model; lo.pred = predict loess

plot(bf.red$chestc           #Chest circumference
 ,residuals(mod.nch))    #Residuals from height ~ weight model

loess <- loess(mod.nch$residuals ~ bf.red$chestc)
lo.pred <- predict(loess, se=T)

lines(bf.red$chestc,lo.pred$fit,pch=2) #Main line
lines(bf.red$chestc,lo.pred$fit+2*lo.pred$s, lty=2) #rough & ready CI
lines(bf.red$chestc,lo.pred$fit-2*lo.pred$s, lty=2)

希望你能帮忙.非常感谢,

垫

Answer 1

lines按照它们出现的顺序连接点,这有时是不合需要的.您可以按如下方式排序:

i <- order(bf.red$chestc)
lines(bf.red$chestc[i], lo.pred$fit[i])
...