Как рассчитать интервал прогнозирования для множественной регрессии МНК?


24

Что такое алгебраическая запись для вычисления интервала прогнозирования множественной регрессии?

Это звучит глупо, но у меня возникают проблемы с нахождением четкой алгебраической записи этого.

Ответы:


31

Возьмем регрессионную модель с наблюдениями и k регрессорами: \ mathbf {y = X \ beta + u} \ newcommand {\ Var} {\ rm Var}Nk

y=Xβ+u

Для данного вектора x0 прогнозируемое значение для этого наблюдения будет

E[y|x0]=y^0=x0β^.
Последовательная оценка дисперсии этого предсказания:
V^p=s2x0(XX)1x0,
где
s2=Σi=1Nu^i2Nk.
Ошибка прогноза для конкретного y0 :
e^=y0y^0=x0β+u0y^0.
Нулевая ковариация между u0 и β^ подразумевает, что
Var[e^]=Var[y^0]+Var[u0],
и непротиворечивая оценка этого
V^f=s2x0(XX)1x0+s2.

1α confidence интервал будет:

y0±t1α/2V^p.
1α prediction интервал будет шире:
y0±t1α/2V^f.


Ответ выше очень хорошо сделан, но я думаю, что этот источник помогает обеспечить некоторый контекст для вопроса.
Июнь Скитер

@Dimitriy Я полагаю, что у вашего второго eqn должна быть морковь / шляпа, '^', над . β
Дон Словик

Не является ли ошибка прогноза остаточной: ? e^=u^
Дон Словик
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.