Zagadnienia

9. Własności estymatorów MNK

9.1. Warunkowa wartość oczekiwana jako rzut ortogonalny

Definicja 9.1

Przypomnienie

Warunkowa wartość oczekiwana:

Niech Y będzie całkowalną zmienną losową w przestrzeni probabilistycznej Ω,F,P, M σ-ciałem takim, że M∈F. Warunkową wartością oczekiwaną Y pod warunkiem M nazywamy zmienną losową E(Y|M), że:

  1. E(Y|M) jest M-mierzalna,

  2. ∀A∈M∫AYdP=∫AE(Y|M)dP.

Załóżmy, że dla przestrzeni probabilistycznej Ω,F,P określone zostały zmienne losowe całkowalne z kwadratem: X,Y ∈ L2⁢Ω,F,P. Zdefiniujmy Y^=E(Y|X) (=E(Y|σ(X))).

Stwierdzenie 9.1

Y^ jest rzutem ortogonalnym Y na L2⁢Ω,σ⁢X,PX, gdzie σ⁢X jest σ-ciałem generowanym przez X, a PX to miara prawdopodobieństwa warunkowego pod warunkiem zmiennej losowej X.

Załóżmy, że Z jest rzutem ortogonalnym Y na σ⁢X. Wtedy ∀ A ∈ σ⁢X:

E⁢1⁢A⁢Y=E⁢1⁢A⁢Z+E⁢1⁢A︸∈σ⁢X⁢Y-Z︸⟂σ⁢X=E⁢1⁢A⁢Z+0

Z definicji warunkowej wartości oczekiwanej Y^=Z p.n.

∎

9.2. Twierdzenie Pitagorasa

Niech X,Y oznaczają zmienne losowe, X,Y ∈ L2⁢Ω,F,P, X ∈ L0=L2⁢Ω,σ⁢X,PX. Zdefiniujmy iloczyn skalarny jako <X,Y≥Ep(XY).

Twierdzenie 9.1

Pitagorasa

Y-X2=Y-Y^2+Y^-X2.
Y-X2=Y-Y^+Y^-X2=Y-Y^2+Y^-X2+2<Y-Y^︸⟂σ⁢X,Y^-X︸∈σ⁢X≥
=Y-Y^2+Y^-X2+0.
∎
Wniosek 9.1

Rzut ortogonalny jest więc najlepszym przybliżeniem Y w klasie f ∈ L0 w sensie:

Y-Y^2=minf∈L0⁡Y-f2
Stwierdzenie 9.2

Jeśli f ∈ L0, to<Y,f≥<Y^,f>.

<Y,f≥<Y-Y^︸⟂σ⁢X,f︸∈σ⁢X>+<Y^,f≥0+<Y^,f>.
∎
Stwierdzenie 9.3

Niech Y^^ oznacza rzut Y^ na L1⊆L0. Y^^ jest rzutem Y na przestrzeń L1. Rzut rzutu jest rzutem.

∀f∈L1<Y-Y^^,f≥<Y,f>-<Y^^,f≥

korzystając ze stwierdzenia 9.2,

=<Y^,f>-<Y^^,f≥<Y^-Y^^,f≥0

z założenia.

∎
Stwierdzenie 9.4

Oznaczmy Y¯ jako rzut ortogonalny Y na lin1, Y^ rzut Y na L=lin⁢X,1. Wtedy:

<Y-Y¯,Y^-Y¯≥||Y^-Y¯||2.

Ponieważ Y^-Y¯ ∈ L0, równość wynika łatwo ze stwierdzenia 9.2.

∎
Definicja 9.2

Przypomnienie

Korelacja:

Cor⁢X,Y=Cov⁢X,Yσ⁢X⁢σ⁢Y=E⁢X-E⁢X⁢Y-E⁢Y Var⁢X⁢Var⁢Y.

Współczynnik dopasowania R2 to część zmienności Y wyjaśnionej przez zmienność Y^:

R2=Var⁢Y^Var⁢Y.

Błąd średniokwadratowy między X i Y:

E⁢X-Y2.
Twierdzenie 9.2
Cor2⁢Y,Y^=R2=Y^-Y¯2Y-Y¯2=1-Y^-Y2Y-Y¯2.
Cor2⁢Y,Y^=<Y-Y¯,Y^-Y¯>2<Y-Y¯,Y-Y¯><Y^-Y¯,Y^-Y¯>=

Korzystając z 9.4:

=Y^-Y¯2Y-Y¯2=

Korzystając z 9.2 dla X=Y¯:

=1-Y^-Y2Y-Y¯2.
∎
Twierdzenie 9.3

Pitagrorasa dla korelacji

∀f∈σ⁢X⁢⁢Cor2⁢Y,f=Cor2⁢Y,Y^⁢Cor2⁢Y^,f.

Załóżmy, że Y=Y-Y¯, Y^=Y^-Y¯, f=f-f¯ (centrujemy zmienne losowe).

Cor⁢Y,f=<Y,f><Y,Y>12<f,f>12=

Korzystając z 9.2 i 9.4, mamy:

=<Y,Y^><Y,Y>12<Y^,Y^>12⋅<Y^,f><Y^,Y^>12<f,f>12=
=Cor2⁢Y,Y^⁢Cor2⁢Y^,f.
∎
Wniosek 9.2

Największą korelację ze wszystkich f∈σ⁢X, Y ma ze swoim rzutem ortogonalnym na przestrzeń rozpiętą przez X:

maxf∈σ⁢X⁡Cor2⁢Y,f=Cor2⁢Y,Y^.
Wniosek 9.3

Patrząc na wnioski 9.1 i 9.2 oraz twierdzenie 9.2 zauważmy, że zachodzi zależność:

Minimalizacja błędu średniokwadratowego ⇔ Maksymalizacja kwadratu korelacji. Dla równoważnych problemów optymalnym jest rzut ortogonalny.

Obrazuje tą zależność także kolejne stwierdzenie:

Stwierdzenie 9.5

Niech X,Y będą wystandaryzowanymi zmiennymi losowymi (E⁢X=0, E⁢Y=0, σ⁢X=1, σ⁢Y=1 ). Wtedy:

E⁢X-Y2=2-2⁢Cor⁢X,Y.
E⁢X-Y2=E⁢X2︸=1+E⁢Y2︸=1-2⁢E⁢X⁢Y=2-2⁢Cor⁢X,Y.
∎

Z twierdzenia 9.1 wynika jeszcze bardzo ważna zależność znana z rachunku prawdopodobieństwa:

Wniosek 9.4
Var(Y)=Var(E(Y|X))+E(Var(Y|X)),

gdzie Var(Y|X)=E((Y-EY)2|X).

W twierdzeniu 9.1 za X podstawmy E⁢Y.

∎

Treść automatycznie generowana z plików źródłowych LaTeXa za pomocą oprogramowania wykorzystującego LaTeXML.

Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego.

Projekt współfinansowany przez Ministerstwo Nauki i Szkolnictwa Wyższego i przez Uniwersytet Warszawski.