r - 黄土线没有正确绘制

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])
...