黄土の滑らかなプラスの信頼限界を残差の散布図に当てはめるのに問題があります。
私のモデルは身長～体重+胸囲です。胸囲の直線性を確認するために、胸囲のないモデル (身長 ~ 体重) を当てはめ、このモデルの残差を胸囲に対してプロットしました。ここまでは順調ですね。次に、黄土線と信頼限界を使用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)
お役に立てれば幸いです。どうもありがとう、
マット