У меня возникли проблемы с подгонкой лёссового сглаживания плюс доверительные пределы к диаграмме рассеяния остатков.
Моя модель: рост ~ вес + окружность груди. Чтобы проверить линейность окружности груди, я подогнал модель без окружности груди (т. е. рост ~ вес) и нанес на график невязки этой модели в зависимости от окружности груди. Все идет нормально. Затем я попытался использовать 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)
Надеюсь, вы можете помочь. Большое спасибо,
Мат