Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Econometrics

.pdf
Скачиваний:
41
Добавлен:
28.01.2022
Размер:
4.46 Mб
Скачать

Використовуючи оператор оцінювання

 

 

 

 

 

(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

 

 

 

Це означає, що параметр a1

у моделі (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 a0 a1xt .

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 — інструментальна змінна, порядковий номер, a0 y a1x 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2

 

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

 

 

 

Оскільки стандартна похибка оцінки параметра a1 становить

близько 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 р. включно).

a1Et 1

a2Vt a3Iu .

Et a0

 

 

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

Соседние файлы в предмете Моделирование