линия лёсса отображается неправильно

У меня возникли проблемы с подгонкой лёссового сглаживания плюс доверительные пределы к диаграмме рассеяния остатков.

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

Надеюсь, вы можете помочь. Большое спасибо,

Мат

5
задан plannapus 14 November 2013 в 07:39
поделиться