Wie berechnet man das Vorhersageintervall für eine OLS-Mehrfachregression?


Antworten:


31

Nehmen Sie ein Regressionsmodell mit Beobachtungen und Regressoren: Nk

y=Xβ+u

Unter der Annahme eines Vektors der vorhergesagte Wert für diese Beobachtung E [y \ vert \ mathbf {x_0}] = \ hat y_0 = \ mathbf {x_0} \ hat \ beta. Ein konsistenter Schätzer für die Varianz dieser Vorhersage ist \ hat V_p = s ^ 2 \ cdot \ mathbf {x_0} \ cdot (\ mathbf {X'X}) ^ {- 1} \ mathbf {x'_0}, wobei s ^ 2 = \ frac {\ Sigma_ {i = 1} ^ {N} \ hat u_i ^ 2} {Nk}. Der Vorhersagefehler für ein bestimmtes y_0 ist \ hat e = y_0- \ hat y_0 = \ mathbf {x_0} \ beta + u_0- \ hat y_0. Die Null-Kovarianz zwischen u_0 und \ hat \ beta impliziert, dass \ Var [\ hat e] = \ Var [\ hat y_0] + \ Var [u_0] ist, und ein konsistenter Schätzer dafür ist x0

E[y|x0]=y^0=x0β^.
V p = s 2 ⋅ x 0 ⋅( X ' X ) - 1 x ' 0 , s 2 = Σ N i = 1 u 2 i
V^p=s2⋅x0⋅(X′X)−1x0′,
s2=Σi=1Nu^i2N−k.
y0 e =y0-
e^=y0−y^0=x0β+u0−y^0.
u0β^
Var[e^]=Var[y^0]+Var[u0],
V^f=s2⋅x0⋅(X′X)−1x0′+s2.

Das 1−α confidence ist:

y0±t1−α/2⋅V^p.
Das 1−α prediction ist breiter:
y0±t1−α/2⋅V^f.


Die obige Antwort ist sehr gut gemacht, aber ich denke, diese Quelle bietet einen Kontext für die Frage.
— June Skeeter

@Dimitriy Ich glaube, deine zweite Gleichung sollte eine Karotte / einen Hut haben, '^', über der . β
— Don Slowik

Ist der Prognosefehler nicht der Rest: ? e^=u^
— Don Slowik
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.