Множественная линейная регрессия: построение прямой линии с доверительными интервалами

Question

Вот мой вопрос:

1) Я провел множественную линейную регрессию: предположим, что:

lm(attitude~quality+price+location+Income)

Меня больше всего интересует связь между отношением и качеством, а другие переменные являются контрольными переменными.

2) Тогда я хотел сделать разброс между отношением и качеством. Это просто:

Q <-ggplot(data=data, aes(x=quality, y=attitude)) 
Q + geom_point(size = 1)

3) Далее я хотел бы построить подгоночную линию между x и y, и наклон должен быть коэффициентом частичной регрессии из множественной линейной регрессии. То есть это должен быть b1 в следующей формуле: отношение = b1 * качество + b2 * цена + b3 * местоположение + b4 * доход, а не b в следующей формуле: отношение = b * качество. Следовательно, следующий код не будет работать правильно, так как будет отображать наклон b, а не b1.

g <- g + geom_smooth(method = lm)

Кто-то задал очень похожий вопрос, см. Здесь

Предоставленный ответ выглядит следующим образом (заменено моими переменными):

g <- g + geom_smooth(data=data, aes(x=quality, y=attitude, ymin=lcl, ymax=ucl))

Однако, это НИЗКИЙ график (как вы можете видеть на рисунке, размещенном в посте), а не линейный график.

Мой вопрос: как добавить прямую линию наклона b1 с полосой доверительного интервала?

Answer 1

Если вы хотите увидеть b1, вы должны нарисовать графики частичной регрессии , как я знаю.

В этом случае

1. регресс attitude по каждой переменной , за исключением quality
2. регресс quality на других предикторах

3. График остатков этих подгонок:

   X <-
     data_frame(
       x = lm(quality ~.-attitude, data = data)$resid,
       y = lm(attitude ~ .-quality, data = data)$resid
     )
   X %>%
     ggplot(aes(x, y)) +
     geom_smooth(method = "lm")
   

Эта строка может быть такой же, как b1, но не x, y точек.