Econometrics
.pdfВикористовуючи оператор оцінювання |
|
|
|
|
|
(9.15) |
|||||
|
|
A (Z X ) |
1 |
Z Y , |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
маємо: |
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
0 |
1 |
|
|
ny |
|
|
||
A |
|
n |
x2 x1 |
|
|
y |
|
|
|
, |
(9.16) |
0 |
|
n |
|
y |
|
||||||
|
|
2 |
|
|
2 |
|
2 |
1 |
|
|
|
де x2 і x1 характеризують середні відхилення значень X відпові-
дно вгору і вниз від медіани, а |
y2 і y1 — середні значення залеж- |
|||||||||
ної змінної, які відповідають середнім x2 |
і x1 . Звідси |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
||
|
|
|
|
|
|
|
. |
|
(9.17) |
|
|
|
|
A |
|
y |
|
||||
|
|
|
|
y |
|
|
|
|||
|
|
|
|
|
2 |
1 |
|
|
|
|
|
|
|
|
x |
x |
|
|
|||
|
|
|
|
|
2 |
1 |
|
|
|
|
Це означає, що параметр a€1 |
у моделі (9.14) подається у такому |
|||||||||
вигляді: |
|
|
|
|
|
|
|
|
|
|
a€ |
|
y2 |
y1 |
, причому a€ |
y a x. |
(9.18) |
||||
|
|
|||||||||
1 |
|
x2 |
x1 |
|
|
|
|
0 |
1 |
|
|
|
|
|
|
|
|
|
|
||
Коли вибіркова сукупність містить непарну кількість спостережень, то перш ніж розпочинати обчислення, необхідно відкинути середнє спостереження.
За загальних припущень оцінка, яку здобуто методом Вальда, є обґрунтованою, але її вибіркова дисперсія може бути досить великою, тобто оцінка є неефективною.
9.3.2. Особливості оцінювання методом Бaртлета.
Бартлет показав, що ефективність оцінки можна збільшити, якщо розбити впорядковані значення змінної X на три групи однакового розміру. Перша з них містить найменші значення X, друга — середні, а третя — найбільші. Вилучивши середню групу — n/3 спостережень, дістанемо оцінку для параметра:
a€ |
|
y3 |
y1 |
, |
(9.19) |
|
|
||||
1 |
|
x3 |
x1 |
|
|
|
|
|
|||
348
де x3 , x1 , y3 , y1 — середні величини двох крайніх груп спостере-
жень. Вільний член оцінюється так само, як і в (9.18).
Поділ вибіркової сукупності спостережень на три рівні групи доцільно використовувати у прикладних дослідженнях, оскільки немає змоги дістати точну інформацію про закон розподілу значень X.
9.3.3. Оператор оцінювання Дарбіна. Дарбін запропо-
нував упорядковувати значення вектора X в порядку зростання, і ввів як інструментальну змінну порядковий номер (ранг), тобто числа 1, 2, 3, 4, ... n. Учений показав, що для великих вибіркових сукупностей ефективність оцінок параметрів моделі становитиме майже 96% від ефективності оцінок 1МНК, а для сукупностей n=20 ефективність оцінок згідно з оператором Дарбіна становить близько 86%.
Модель Дарбіна не має вільного члена. Щоб застосувати цей метод для оцінювання всіх параметрів моделі, зокрема й для вільного члена, матриці змінних Z' i X' подамо у вигляді:
1 |
1 |
1 |
1 ... |
1 |
|
Z |
2 |
3 |
4 ... |
n |
|
1 |
|
||||
і
|
1 |
1 |
1 |
1 |
... |
1 |
|
, |
||||
X x |
x |
2 |
x |
3 |
x |
4 |
... |
x |
|
|
||
|
1 |
|
|
|
|
|
n |
|
||||
де xi характеризують відхилення від середнього значення |
x, які |
|||||||||||||||||
впорядковуються за зростанням. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Оператор оцінювання |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a€ |
|
iy |
i |
, |
|
|
|
|
|
|
|
|
|
(9.20) |
|||
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|||||||
|
n |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ixi |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
де і — порядковий номер, причому |
|
€ |
|
y |
€ |
|
. Коваріація оці- |
|||||||||||
a0 |
|
a1x |
||||||||||||||||
нок параметрів запишеться: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
€ |
2 |
|
|
X ) |
1 |
(Z |
|
Z )( X |
|
Z ) |
1 |
. |
|
(9.21) |
||||
asycov(a) |
u (Z |
|
|
|
|
|
|
|||||||||||
Метод Дарбіна можна застосовувати і тоді, коли модель містить декілька пояснювальних змінних. У такому разі спочатку
349
знаходяться відхилення значень кожної змінної од відповідного середнього значення. Потім ці відхилення упорядковуються за зростанням і кожному з них присвоюється порядковий номер.
9.4. Помилки вимірювання змінних
Раніше ми припускали, що змінні вимірюються без помилок, і лише відхилення u — це єдина припустима форма помилок. Останнє було пов’язане з наміром врахувати вплив різних пояснювальних змінних, які не входять до економетричної моделі в явному вигляді.
Проте досить часто під час вимірювання змінних, які належать до економетричної моделі, припускаються помилок. Тоді постає запитання, як наявність помилок змінних може вплинути на оцінку параметрів моделі?
Щоб відповісти на це запитання, розглянемо матрицю незалежних змінних X, елементи якої містять помилки.
Нехай
~ |
V , |
(9.22) |
X X |
де ~ — матриця розміром n×m справжніх (фактичних) значень,
X
а V — матриця помилок вимірювання. Тоді модель має вигляд
~
Y Xa u,
або з (9.22) |
~ |
X V , тоді |
|
X |
|
||
|
|
Y XA (u VA). |
(9.23) |
Оцінка параметрів для цієї моделі 1МНК матиме вигляд
A A ( X X ) 1 X (u VA),
де (X X ) 1 X (u VA) — величина зміщення оцінки.
Обґрунтованість цієї оцінки залежить від того, чи дорівнює границя за ймовірністю нулю
|
1 |
|
|
p lim |
|
X (u VA) . |
|
n n |
|
|
|
350
Запишемо
1 |
|
|
1 |
|
1 |
|
|
|
p lim |
|
|
p lim |
plim |
|
|||
X (u VA) |
X u |
X V A. |
||||||
n n |
|
|
n n |
|
n n |
|
|
|
Якщо «справжні» значення змінних X і їх помилки не корелюють із залишками, то
|
1 |
|
|
0 . |
(9.24) |
p lim |
|
X u |
|||
n |
|
|
|
|
|
Але насправді границя за ймовірністю (9.24) дорівнює:
1 |
|
|
1 |
|
|
1 |
|
||
p lim |
|
p lim |
|
|
p lim |
|
VV ' . |
||
X V |
n |
X V |
|||||||
n n |
|
|
n |
|
|
n n |
|
||
Отже, навіть тоді, коли помилки вимірювання змінних X не корелюють зі справжніми значеннями цих змінних, і перший доданок у правій частині дорівнює нулю, другий доданок, який характеризує матрицю коваріацій помилок, здебільшого не дорівнює нулю. А це означає, що за наявності помилок вимірювання змінних оцінка параметрів моделей 1МНК є необґрунтованою і асимптотичне зміщення визначається формулою
|
1 |
|
1 |
1 |
|
|
|
p lim ( A A) p lim |
|
p lim |
|
|
|||
X X |
V V A. |
||||||
n |
n n |
|
|
n n |
|
|
|
Наприклад, якщо ми оцінюємо параметри моделі з двома змінними 1МНК, то зміщення параметра a:
|
|
p lim (a€ |
a) |
|
2a |
|
, |
|
|
||
|
|
|
v |
|
|
|
|||||
|
|
2 |
|
2 |
|
|
|||||
|
|
n |
|
|
|
~ |
v |
|
|
|
|
або |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
||
|
|
p lim a€ |
|
, |
|
|
|
(9.25) |
|||
|
|
|
2 |
|
|
|
|||||
|
|
n |
|
1 2v |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
~ |
|
|
|
|
|
де |
|
|
|
|
|
x |
|
|
|
~ |
— дисперсія |
2 |
— дисперсія помилки вимірювання X, а |
||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
v |
|
|
|
|
|
|
|
|
x |
|
справжніх значень X, причому припускаємо, що помилки вимірювання не корелюють із цими значеннями X.
Рівняння (9.25) показує, що оцінка справжнього значення параметра моделі занижена. Наприклад, хоч би як ми збільшували
351
сукупність спостережень, коли |
|
2 |
=10% від |
~ , то оцінка пара- |
|
|
|
|
|
2 |
|
метра a€ |
|
|
v |
|
x |
відрізнятиметься від справжнього значення також майже |
|||||
на 10%, тобто за наявності помилок вимірювання змінних збільшення сукупності спостережень не компенсує зміщення.
Тому при оцінюванні параметрів економетричної моделі, коли трапляються помилки вимірювання змінних, доцільно застосувати метод інструментальних змінних, який ми розглянули раніше.
Приклад 9.1. Побудувати економетричну модель методом інструментальних змінних, яка характеризує залежність між зайнятістю населення і виробництвом продукції, скориставшись даними, наведеними в табл. 9.1.
Таблиця 9.1
Рік |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
|
|
|
|
|
|
|
|
|
|
Виробництво, |
130 |
128 |
194 |
157 |
195 |
205 |
142 |
225 |
168 |
133 |
гр. од. |
|
|
|
|
|
|
|
|
|
|
Зайнятість, чол. |
114 |
96 |
134 |
112 |
113 |
144 |
105 |
150 |
109 |
110 |
|
|
|
|
|
|
|
|
|
|
|
Розв’язання
1. Ідентифікація змінних і специфікація моделі.
Нехай xt — виробництво продукції, незалежна (пояснювальна) змінна, yt— зайнятість населення, залежна змінна.
Економетрична модель має вигляд yt=a0+a1хt+ut;
yt a€0 a€1xt .
2. Оцінка параметрів моделі.
Оскільки вихідні дані можуть мати помилки вимірювання, то для оцінювання параметрів моделі застосуємо метод інструментальних змінних.
2.1. Оператор оцінювання за методом інструментальних змінних
€ |
|
X ) |
1 |
|
A (Z |
|
Z Y , |
||
352
де Z — матриця інструментальних змінних; X — матриця пояснюючих змінних.
2.2. Визначимо матрицю інструментальних змінних за методом Дарбіна. Для цього впорядкуємо значення вектора X від меншого до більшого і надамо кожному елементу цього вектора порядковий номер. Запишемо ці дані в табл.9.2.
|
|
|
|
|
Таблиця 9.2 |
|
|
|
|
|
|
Y |
X |
Інструментальна |
€ |
€ |
€2 |
|
|
змінна Z |
Y |
u |
u |
|
|
|
|
|
|
96 |
128 |
1 |
91,6398 |
4,360 |
19,01 |
|
|
|
|
|
|
114 |
130 |
2 |
93,0030 |
20,997 |
440,87 |
|
|
|
|
|
|
110 |
133 |
3 |
95,0478 |
14,95 |
223,5 |
|
|
|
|
|
|
105 |
142 |
4 |
101,1822 |
3,81 |
14,51 |
|
|
|
|
|
|
112 |
157 |
5 |
111,4062 |
0,5938 |
0,352 |
|
|
|
|
|
|
109 |
168 |
6 |
118,9000 |
–9,9 |
98,01 |
|
|
|
|
|
|
134 |
194 |
7 |
136,6254 |
–2,625 |
6,89 |
|
|
|
|
|
|
113 |
195 |
8 |
137,3000 |
–24,3 |
590,49 |
|
|
|
|
|
|
144 |
205 |
9 |
144,1230 |
–0,123 |
0,015 |
|
|
|
|
|
|
150 |
225 |
10 |
157,7550 |
–7,755 |
60,14 |
|
|
|
|
|
|
|
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
1 |
|
Z |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
. |
Оскільки економетрична модель має вільний член, перший рядок матриці інструментальних змінних складається з одиниць,
353
а другий є інструментальною змінною замість вектора пояснювальної змінної X.
1 1281 1301 133
1 142
Матриця X має такий вигляд: X 1 157 .
1 1681 1941 1951 2051 225
2.3. Знайдемо оцінки параметрів моделі
|
10 |
|
|
|
|
|
|
|
i y |
|
|
6926 0,6816, |
|
a€ |
|
t 1 |
t |
t |
|
|
|
|
|
||||
1 |
10 |
|
|
|
10161 |
|
|
|
it xt |
|
|||
|
|
t 1 |
|
|
|
|
де іt — інструментальна змінна, порядковий номер, a€0 y a€1x 118,7 0,6816 167,8 4,395.
Звідси економетрична модель запишеться так: yt 4,395 0,681xt .
2.4. Підставимо значення змінної xt в модель і знайдемо розрахункові значення залежної змінної — зайнятості населення ( yt ).
Ці значення і відхилення їх від фактичних наведено в табл. 9.2. 3. Визначимо коваріаційну матрицю оцінок параметрівмоделі:
|
|
|
|
|
|
|
€ |
|
|
2 |
|
1 |
(Z Z )( X Z ) |
1 |
. |
|
|
|
|
|
|
asy cov ( A) |
u (Z X ) |
|
|
|
|||||||||
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u€2 |
|
1453,787 |
|
|
|
|
|
|
|
|
||||
3.1. |
u2 |
|
t 1 |
t |
|
|
181,73. |
|
|
|
|
|
|
||||
n m |
8 |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
3.2. |
Z |
|
10 |
|
1677 |
; |
|
(Z X ) |
1 |
|
|
1,08384 |
|
0,17888 |
|||
X |
55 |
|
|
|
|
|
0,0058666 |
0,001060 |
. |
||||||||
|
|
|
|
10161 |
|
|
|
|
|
|
|
||||||
354
3.3. |
Z Z |
10 |
55 |
|
|
|
|
|
|
|
|
|
|
|
|
55 |
. |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
385 |
|
|
|
|
|
|
|
|
|
|
|
3.4. |
|
|
10 |
55 |
|
|
; ( X Z ) |
1 |
|
1,08384 |
0,0058666 |
|||
X Z |
|
|
|
|
|
|
0,17888 |
0,001066 |
. |
|||||
|
|
1677 10161 |
|
|
|
|
|
|
||||||
3.5. |
(Z X ) |
1 |
(Z |
|
|
|
1 |
9,26 |
|
|
|
|
||
|
Z ) |
|
|
|
0,0877 |
. |
|
|
|
|||||
|
|
|
|
0,000036 |
|
|
|
|
||||||
3.6. |
(Z X ) |
1 |
(Z Z )( X Z ) |
1 |
|
2,74 |
0,004 |
|
|
|||||
|
|
|
|
|
|
|
0,00009 |
. |
|
|
||||
|
|
|
|
|
|
|
0,0157 |
|
|
|
||||
3.7. |
|
|
€ |
|
|
|
2,74 |
|
0,004 |
|
|
|
||
asy cov ( A) 181,73 |
|
0,0157 |
|
|
. |
|
|
|
||||||
|
|
|
|
|
|
|
0,00009 |
|
|
|
||||
4. Знайдемо стандартні похибки оцінок параметрів моделі
Sa€ |
181,73 2,74 |
|
497,9402 22,3; |
0 |
|
|
|
Sa€ |
181,73 0,00009 |
|
0,0163557 0,13. |
1 |
|
|
|
Оскільки стандартна похибка оцінки параметра a€1 становить
близько 17% до її абсолютного значення, а стандартна похибка вільного члена перевищує його абсолютну величину в кілька разів, то можна стверджувати про наявність зміщення цих оцінок параметрів моделі та їх неефективність.
9.5. Економетрична модель експорту продукції із стохастичними пояснювальними змінними
Розглянемо динамічну економетричну модель експорту продукції, яка в загальному вигляді запишеться так:
Et f Et 1,Vt , Iut ,ut ,
де Et — експорт продукції в період t;
Et-1 — експорт продукції в період t – 1;
Vt — обсяг виробництва продукції в період t;
Iut — індекс співвідношення цін на продукцію на зовнішньому та внутрішньому ринках;
355
ut — стохастична складова, яка характеризує вплив на експорт продукції випадкових чинників у період t.
Специфікуємо цю модель у лінійній формі:
Et=a0+a1Et–1+a2Vt+a3Iut+ut. (9.26)
Пояснювальні змінні даної моделі є стохастичними, бо експорт продукції в період t (Et) залежить від стохастичної складової ut, а експорт продукції в період t – 1 (Et–1) залежить від стохастичної складової ut–1. Звідси, оскільки Et залежить від Et–1, то Et–1 залежить від ut.
Цей взаємозв’язок вимагає розглядати всю матрицю пояснювальних змінних як стохастичну і спонукає висунути гіпотезу щодо залежності матриці пояснювальних змінних і залишків. Тому для оцінювання параметрів даної моделі потрібно застосувати метод інструментальних змінних або двокроковий метод найменших квадратів, який буде розглянуто пізніше.
Побудуємо цю економетричну модель, аналізуючи експорт цукру з України в Росію на основі статистичної інформації за 12 років (з 1990 по 2001 р. включно).
€ |
a€1Et 1 |
a€2Vt a€3Iu . |
Et a€0 |
||
|
|
t |
У даній моделі замість лагової залежної змінної Et–1 використаємо інструментальну змінну Zt.
Для визначення цієї інструментальної змінної знайдемо відхилення змінної Et–1 від своєї середньої, а потім проранжуємо величину відхилень від меншого до більшого значення. Пронумеруємо ранжовані значення відхилень змінної Et–1 і цей порядковий номер використаємо в економетричній моделі експорту цукру як інструментальну змінну (метод Дарбіна).
Запишемо економічну інформацію в табл. 9.3.
|
|
|
Таблиця 9.3 |
|
|
|
|
|
|
Роки |
Обсяг експорту цукру |
Обсяг виробництва |
Індекс співвідношен- |
|
(млн т) |
(млн т) |
ня цін на цукор |
||
|
||||
|
|
|
|
|
1 |
3200 |
4900 |
0,122 |
|
|
|
|
|
|
2 |
3500 |
4200 |
0,675 |
|
|
|
|
|
|
3 |
2500 |
3500 |
1,205 |
|
|
|
|
|
356
4 |
2475 |
3600 |
1,274 |
|
|
|
|
5 |
2633 |
3900 |
0,851 |
|
|
|
|
6 |
1800 |
3400 |
1,052 |
|
|
|
|
7 |
1700 |
3900 |
1,231 |
|
|
|
|
8 |
1054 |
3300 |
1,18 |
|
|
|
|
9 |
740 |
2000 |
1,201 |
|
|
|
|
10 |
157 |
1900 |
1,131 |
|
|
|
|
11 |
69 |
1600 |
1,401 |
|
|
|
|
12 |
–400 |
1540 |
1,425 |
Ідентифікуємо змінні моделі:
Et — обсяг експорту цукру в період t, залежна змінна;
на;
Et–1 — обсяг експорту цукру в період t–1, пояснювальна змін-
Vt — обсяг виробництва цукру в період t, пояснювальна змінна; Iut — співвідношення цін на цукор на російському та україн-
ському ринках.
Сукупність спостережень зменшиться на одиницю, бо за умови лагової змінної Et–1 експорт в період t потрібно подати без першого спостереження, а експорт у період t–1 — без останнього. Звідси обсяг виробництва цукру в період t та індекс цін потрібно розглядати без першого спостереження.
Запишемо ці змінні в табл. 9.4.
Таблиця 9.4
Рік |
Et |
Et–1 |
Vt |
Iut |
Рік |
Et |
Et–1 |
Vt |
Iut |
2 |
3500 |
3200 |
4200 |
0,675 |
8 |
1054 |
1700 |
3300 |
1,18 |
|
|
|
|
|
|
|
|
|
|
3 |
2500 |
3500 |
3500 |
1,205 |
9 |
740 |
1054 |
2000 |
1,202 |
|
|
|
|
|
|
|
|
|
|
4 |
2475 |
2500 |
3600 |
1,274 |
10 |
157 |
740 |
1900 |
1,131 |
|
|
|
|
|
|
|
|
|
|
5 |
2633 |
2475 |
3900 |
0,851 |
11 |
69 |
157 |
1600 |
1,401 |
|
|
|
|
|
|
|
|
|
|
6 |
1800 |
2633 |
3400 |
1,052 |
12 |
400 |
69 |
1540 |
1,425 |
|
|
|
|
|
|
|
|
|
|
7 |
1700 |
1800 |
3900 |
1,231 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
357
