SZACOWANIE PARAMETRÓW STRUKTURALNYCH MODELU EKONOMETRYCZNEGO Z JEDNĄ ZMIENNĄ OBJAŚNIAJĄCĄ METODĄ NAJMNIEJSZYCH KWADRATÓW
Lp. ( to np. t) |
yt |
x1t |
yt* |
---|---|---|---|
1 | 5 | 3 | 4,745 |
2 | 5 | 3 | 4,745 |
3 | 8 | 5 | 7,995 |
4 | 6 | 4 | 6,37 |
5 | 7 | 5 | 7,995 |
6 | 10 | 6 | 9,62 |
7 | 10 | 6 | 9,62 |
8 | 10 | 6 | 9,62 |
9 | 11 | 7 | 11,245 |
Rysujemy wykres , sprawdzamy zależność
Model będzie miał taką postać
Yyt = α1Xx1t + α0 + ξt
Empiryczny wykres rozrzutu .
Z wykresu można wywnioskować, że:
- zależność jest liniowa ,
- zależność jest dodatnia , bo wraz ze wzrostem x wzrasta y,
- siła zależności ( gdy punkty zbite to duża , gdy rozrzucone mała lub brak).
Dodatnia zależność mówi, że α1 musi być też dodatnie ( spodziewamy się + dla α1). W przypadku gdyby zależność była dodatnia , a α1 ujemne to oznaczałoby to brak zależności , a musi być zależność przyczynowo – skutkowa.
dana jest formuła a = ( xTx)−1xTy
Aby wyliczyć pkt 2. sporządzamy macierze y, x, transponujemy macierz x, liczymy iloczyn (xTx), liczymy wyznacznik (xTx) , liczymy macierz (xTx) −1 , liczymy iloczyn (xTy) , (zaokrąglać do 4 miejsca po przecinku)
ślepa zmienna
$Y = \begin{bmatrix}
5 \\
5 \\
8 \\
6 \\
7 \\
10 \\
10 \\
10 \\
11 \\
\end{bmatrix}$ $X = \begin{bmatrix}
3 & 1 \\
3 & 1 \\
5 & 1 \\
4 & 1 \\
5 & 1 \\
6 & 1 \\
6 & 1 \\
6 & 1 \\
7 & 1 \\
\end{bmatrix}$ $X^{T}\ = \begin{bmatrix}
3 & 3 & 5 & 4 & 5 & 6 & 6 & 6 & 7 \\
1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 \\
\end{bmatrix}$ $X^{T}X\ = \begin{bmatrix}
241 & 45 \\
45 & 9 \\
\end{bmatrix}$
det ( XTX) = 241 × 9 − ,45 × 45 = 144
${{(X}^{T}X)\ }^{- 1} = \ \frac{1}{144}\ D^{T}$
D11 = ( − 1) × 9 = 9
D12 = ( − 1) × 45 = −45 zatem $D = \ \begin{bmatrix} 9 & - 45 \\ - 45 & 241 \\ \end{bmatrix} = \ D^{T}\ $
D21 = ( − 1) × 45 = −45
D22 = ( − 1)4 × 241 = 241
Uwaga xTx = $\begin{bmatrix} 241 & 45 \\ 45 & 9 \\ \end{bmatrix}$ oraz D = DT = $\begin{bmatrix} 9 & - 45 \\ - 45 & 241 \\ \end{bmatrix}$ są odwrotnościami
i należy bez liczenia główną przekątną napisać na odwrót , a do liczb na drugiej przekątnej (-) dodać
zatem
${{(X}^{T}X)\ }^{- 1} = \ \frac{1}{144}\ \times \ \begin{bmatrix} 9 & - 45 \\ - 45 & 241 \\ \end{bmatrix}$ ${{(X}^{T}X)\ }^{- 1} = \ \begin{bmatrix} 0,0625 & - 0,3125 \\ - 0,3125 & 1,6736 \\ \end{bmatrix}$
Dalej (xTy)
$X^{T}\ = \ \begin{bmatrix} 3 & 3 & 5 & 4 & 5 & 6 & 6 & 6 & 7 \\ 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 & 1 \\ \end{bmatrix}$ $Y = \ \begin{bmatrix} 5 \\ 5 \\ 8 \\ 6 \\ 7 \\ 10 \\ 10 \\ 10 \\ 11 \\ \end{bmatrix}$ $(X^{T}Y) = \begin{bmatrix} 386 \\ 72 \\ \end{bmatrix}$
wszystkie uzyskane dane wstawiamy do równania a = ( xTx)−1xTy
$a = \begin{bmatrix} 0,0625 & - 0,3125 \\ - 0,3125 & 1,6736 \\ \end{bmatrix}\text{\ \ \ } \times \begin{bmatrix} 386 \\ 72 \\ \end{bmatrix}\ = \ \text{\ \ }\begin{bmatrix} 1,625 \\ - 0,13 \\ \end{bmatrix}$ więc $a = \begin{bmatrix} 1,625 \\ - 0,13 \\ \end{bmatrix}\ $ $\begin{matrix} a_{1} \\ a_{0} \\ \end{matrix}$
czyli model możemy zapisać 1, 625 x1t−0, 13 + ut ( NIE ZAPOMINAĆ DOPISAĆ ut)
Dalej pytamy czy parametr przy x1t jest przyczynowo skutkowy.
Koincydencja x1t
sgn(ryx1t) ≡ sgn(a1) → ≡a1 > 0
Interpretacja :
Wzrost x1t o 1 ( jedną jednostkę) spowoduje wzrost y o 1,625 jednostki.
Interpretujemy parametr wolny a0
a0 = − 0, 13
Taką średnią wartość przyjmuje zmienna endogeniczna w przypadku, gdy zmienna objaśniająca będzie równa 0
X1t = 0 to Yt = − 0, 13 + ut (z dokładnością do składnika resztowego ut)
Można to również odczytać z wykresu
Przykłady linii regresji:
Może się też zdarzyć , że jakieś punkty będą leżeć daleko do linii regresji , to można je pominąć lub spróbować linię regresji zbliżyć do odszczepieńców).
dalej liczymy wartości teoretyczne
Yt* = 1, 625X1t − 0, 13
I wstawiamy kolejne x np.
Y1* = 1, 625 × 3 − 0, 13 = 4, 745
wstawiamy te wartości do 3 kolumny z początku ćwiczeń , rysujemy linię teoretyczną ( regresji)na wykresie z początku ćwiczeń , tzn. weryfikujemy model, czyli sprawdzamy czy jest to świat coś warty, czy spełnia kryteria estymatora (zauważamy , że linie prawie się pokrywają) .
Jak dopasował się Y teoretyczny do rzeczywistośći
Zadanie domowe
(wykres rozrzutu, wnioski, szacowanie parametrów modelu, sprawdzenie koincydencji parametru α1, interpretacja wartości teoretycznych, wykres)
Lp | yt |
xt |
yt* |
---|---|---|---|
1 | 2 | 2 | 3 |
2 | 5 | 4 | 5 |
3 | 4 | 3 | 4 |
4 | 4 | 2 | 3 |
5 | 7 | 6 | 7 |
6 | 2 | 1 | 2 |
Rysujemy wykres , sprawdzamy zależność
Model będzie miał taką postać Yt = α1X1t + α0 + ξt
yt , yt*
x1t
zależność jest liniowa , zależność jest dodatnia , bo wraz ze wzrostem x wzrasta y
dana jest formuła a = ( XTX)−1XTY
sporządzamy macierze Y , X , transponujemy macierz x, liczymy iloczyn (XTX), liczymy wyznacznik (XTX) , liczymy macierz (XTX) −1 , liczymy iloczyn (XTY)
$Y = \begin{bmatrix} 2 \\ 5 \\ 4 \\ 4 \\ 7 \\ 2 \\ \end{bmatrix}$ $X\ = \begin{bmatrix} 2 & 1 \\ 4 & 1 \\ 3 & 1 \\ 2 & 1 \\ 6 & 1 \\ 1 & 1 \\ \end{bmatrix}$ $X^{T}\ = \ \begin{bmatrix} 2 & 4 & 3 & 2 & 6 & 1 \\ 1 & 1 & 1 & 1 & 1 & 1 \\ \end{bmatrix}$ $X^{T}X\ = \begin{bmatrix} 70 & 18 \\ 18 & 6 \\ \end{bmatrix}$
det ( XTX)= 70 × 6 − 18 × 18 = 96 ${{(X}^{T}X)\ }^{- 1} = \ \frac{1}{96}D^{T}$ $D = \ \begin{bmatrix} 6 & - 18 \\ - 18 & 70 \\ \end{bmatrix} = \ D^{T}$
zatem
${{(X}^{T}X)\ }^{- 1} = \ \frac{1}{96}\ \times \ {{(X}^{T}X)\ }^{- 1}\begin{bmatrix} 6 & - 18 \\ - 18 & 70 \\ \end{bmatrix}\ \ = \ \begin{bmatrix} 0,0625 & - 0,1875 \\ - 0,1875 & 0,7291 \\ \end{bmatrix}\ $
Dalej (xTy)
$X^{T}\ = \begin{bmatrix} 2 & 4 & 3 & 2 & 6 & 1 \\ 1 & 1 & 1 & 1 & 1 & 1 \\ \end{bmatrix}$ $Y = \ \begin{bmatrix} 2 \\ 5 \\ 4 \\ 4 \\ 7 \\ 2 \\ \end{bmatrix}$ $X^{T}Y\ = \ \begin{bmatrix} 88 \\ 24 \\ \end{bmatrix}$
wstawiamy do równania a =( XTX)−1XTY
$a = \begin{bmatrix} 0,0625 & - 0,1875 \\ - 0,1875 & 0,7291 \\ \end{bmatrix} \times \ \begin{bmatrix} 88 \\ 24 \\ \end{bmatrix}\ = \ \text{\ \ }\begin{bmatrix} 1 \\ 1 \\ \end{bmatrix}\ $ więc $a = \begin{bmatrix} 1 \\ 1 \\ \end{bmatrix}\ $ $\begin{matrix} a_{1} \\ a_{0} \\ \end{matrix}$
czyli model możemy zapisać 1x1t+ 1 + ut
Dalej pytamy czy parametr przy x1t jest przyczynowo skutkowy.
Koincydencja x1t Sgn(γ yx1t) ≡ sgn (a1)
yt ,
≡a1 >0
Interpretacja :
wzrost X1t o 1 ( jedną jednostkę) spowoduje wzrost y o 1jednostke.
Interpretujemy parametr wolny a0 ∖ na0 = 1
Taką średnią wartość przyjmuje zmienna endogeniczna w przypadku, gdy zmienna objaśniająca będzie równa 0
X1t = 0 to Yt = 1 + ut
dalej liczymy wartości teoretyczne
Yt* = 1X1t + 1
I wstawiamy kolejne x np.
Y1* = 1 × 2 + 1 = 3
wstawiamy te wartości do 3 kolumny z początku zadania , rysujemy linię teoretyczną regresji, tzn. weryfikujemy model, czyli sprawdzamy czy jest to świat coś warty, czy spełnia kryteria estymatora.
MODEL Z DWIEMA ZMIENNYMI OBJAŚNIAJĄCYMI
Lp | Yt |
X1t |
X2t |
---|---|---|---|
1 | 2 | 0 | 0 |
2 | 3 | 1 | 0 |
3 | 2 | 0 | 0 |
4 | 4 | 0 | 1 |
5 | 4 | 0 | 1 |
6 | 3 | 1 | 0 |
7 | 5 | 1 | 2 |
8 | 5 | 1 | 2 |
Yt = α0 + α1X1t + α2X2t + ξt
Sprawdzamy czy możemy szacować , rysujemy dwa wykresy x1, x2 oba muszą być dodatnie
Dalej a = ( X′X)−1X′Y
y= $\begin{bmatrix} 2 \\ 3 \\ 2 \\ 4 \\ 4 \\ 3 \\ 5 \\ 5 \\ \end{bmatrix}$ x =$\ \begin{bmatrix} 1 & 0 & 0 \\ 1 & 1 & 0 \\ 1 & 0 & 0 \\ 1 & 0 & 1 \\ 1 & 0 & 1 \\ 1 & 1 & 0 \\ 1 & 1 & 2 \\ 1 & 1 & 2 \\ \end{bmatrix}$