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

Econometrics

.pdf
Скачиваний:
21
Добавлен:
28.01.2022
Размер:
4.46 Mб
Скачать
У лінійній моделі оцінки параметрів

3,806662

4,248495

3,401197

4,75359

 

 

 

 

3,871201

4,317488

3,526361

4,736198

 

 

 

 

3,931826

4,369448

3,555348

4,787492

 

 

 

 

3,89182

4,343805

3,496508

4,820282

 

 

 

 

3,988984

4,406719

3,610918

4,779123

 

 

 

 

4,007333

4,382027

3,610918

4,859812

 

 

 

 

4,043051

4,317488

3,663562

4,859812

 

 

 

 

4,025352

4,418841

3,637586

4,882802

 

 

 

 

3,988984

4,394449

3,583519

4,867534

 

 

 

 

4,077537

4,465908

3,688879

4,820282

 

 

 

 

4,110874

4,521789

3,73767

4,89784

 

 

 

 

4,127134

4,553877

3,7612

4,919981

 

 

 

 

4,158883

4,574711

3,73767

4,934474

 

 

 

 

Скориставшись функцією програми «Exсel» «ЛИНЕЙН», знайдемо оцінки параметрів:

ln a0

0,044037;

a1

0,495475;

a2

– 0,10808;

1,099585;

a3

 

anti ln a0

1,045.

Отже, економетрична модель прибутку у степеневій формі на-

бирає вигляду Y 1,045X10,495 X 2 0,108 X 31,10 .

Як бачимо, економетрична модель прибутку в лінійній формі відрізняється від моделі у степеневій формі. Ця різниця полягає передусім у тому, що оцінки параметрів в обох моделях мають різний економічний зміст.

a1,a2 ,a3 характеризують граничний приріст прибутку залежно від граничного приросту кожного ресурсу на одиницю (коли решта — сталі) і в тих одиницях, в яких вони подаються у вихідній інформації.

У степеневій моделі оцінки параметрів a1, a2 , a3 характеризу-

ють кількісний зв’язок між прибутком та відповідно кожним ресурсом у відносному (відсотковому) виразі — еластичність. Тому

108

їх потрібно тлумачити так: якщо інвестиції зростуть на 1%, а основні виробничі фонди і фонд робочого часу не зміняться, то прибуток зросте на 0,495% a1 0,495 ; якщо основні виробничі

фонди зростуть на 1%, а решта ресурсів буде сталою, то прибуток зменшиться на 0,108% a2 0,108 ; і, нарешті, якщо фонд

робочого часу зросте на 1%, і решта ресурсів буде сталою, то прибуток зросте на 1,1% a3 1,1 .

Зміну напряму взаємозв’язку між прибутком та основними виробничими фондами у степеневій моделі можна пояснити особливостями статистичної інформації (можлива мультиколінеарність, автокореляція, про які йтиметься далі).

4.6. Коваріаційна матриця оцінок параметрів моделі

За допомогою коваріаційної матриці розраховуються основні показники випадкового розсіювання оцінок aj навколо відповід-

них істинних значень параметрів, що аналізуються, а також харак-

теристики взаємозв’язків отриманих оцінок.

вектор

У

класичній

регресійній

моделі Y=XA+u;

u u1, u2 , ... , un

і залежний від нього вектор Y ( y1, y2 , ..., yn ) є ви-

 

 

 

 

 

 

 

 

падковими змінними. До оператора оцінювання A входить вектор

( X X )

1

X

 

Y

, а отже, оператор

 

Y A

 

 

A також можна вважати випа-

 

 

 

 

 

 

 

 

дковоюфункцією оцінювання параметрівмоделі.

 

Відомо,

що для характеристики випадкових змінних aj

поряд

із математичним сподіванням застосовуються також дисперсія2aj і коваріація ajak (j k). Істинні (справжні) значення цих па-

раметрів класичної економетричної моделі утворюють диспер-

сійно-коваріаційну матрицю

 

 

 

 

 

2

a1a2

...

a1a j

... a1am

 

 

 

 

 

 

 

a1

 

 

 

 

 

 

a a

a2

...

a

a

j

...

a

a

m

 

 

 

 

 

 

2 1

2

 

2

 

 

2

 

 

 

2

( X X )

1

...

... ... ...

 

... ...

 

 

. (4.12)

cov ( A)

u

 

 

a ja1

a ja2

...

a2 j

 

...

a jam

 

 

 

 

 

 

 

 

 

 

 

 

...

... ... ...

 

... ...

 

 

 

 

 

 

 

 

ama1

ama2

...

ama j

...

2

 

 

 

 

 

 

 

 

am

 

 

109

Оцінки коваріаційної матриці cov ( A) €u ( X X )

 

використо-

2

1

 

вуються для знаходження стандартних похибок та обчислення довірчих інтервалів оцінок параметрів aj . Вони використовують-

ся й під час перевірки їхньої статистичної значущості.

На головній діагоналі матриці

містяться оцінки диспер-

cov ( A)

сій 2aj j-ї оцінки параметрів, що ж до елементів ajak (j k), які розміщені поза головною діагоналлю, то вони є оцінками коварі-

ації між aj і ak . Отже,

 

 

 

 

2

a1a2

...

a1aj

...

 

 

 

 

 

a1

 

 

 

 

 

 

 

 

 

 

aa

a2

...

aa

j

...

 

 

 

 

 

2 1

 

 

2

 

 

 

2

 

2

u2 ( X X ) 1

...

 

... ... ... ...

 

Aj

 

 

aja1

aja2

...

 

 

2

 

 

 

 

 

aj ...

 

 

 

...

 

... ... ... ...

 

 

 

 

ama1

ama2

...

amaj

...

 

 

 

 

 

 

 

 

 

a1ama2am

... , (4.13)

ajam

...

2am

де u2 — незміщена оцінка дисперсії залишків;

u2

u' u

.

 

 

n m

Зауважимо, що матриця коваріації оцінок параметрів моделі характеризує також ступінь їх випадкового розсіювання, який обчислюється, як значення детермінанта коваріаційної матриці cov(A).

Детермінант коваріаційної матриці

є так званою

det(cov(A))

узагальненою дисперсією, яка кількісно характеризує ступінь випадкового розсіювання значень векторної випадкової величини навколо свого середнього у відповідному багатовимірному просторі:

 

 

 

Dуз ( A) det (cov( A)) .

(4.14)

Часто використовується й інша характеристика цього випадкового розсіювання значень багатовимірної випадкової величини —

слід коваріаційноїматриці tr:

110

 

a2

a2

... a2.

(4.15)

tr(cov( A)) a2

1

2

3

 

m

 

Виходячи із додатної визначеності матриці

і змісту ді-

cov ( A)

агональних елементів 2aj , можна стверджувати, що величини,

які визначені співвідношеннями (4.14) і (4.15), завжди додатні. Чим більше значення знайдених характеристик (детермінанта, сліду дисперсійно-кореляційної матриці), тим більша загальна варіація оцінок параметрів моделі).

Ступінь тісноти взаємозв’язку між окремими оцінками пара-

метрів моделі aі aвектора краще визначати на основі кое-

k j A

фіцієнта кореляції, який, у свою чергу, визначається через елемен-

ти коваріаційної матриці

 

 

 

 

 

 

cov ( A) :

 

 

 

 

 

r

 

уakaj

 

.

(4.16)

 

 

 

 

 

kj

 

a

a

 

 

 

 

 

 

jj

 

 

 

 

kk

 

 

Усе це говорить про те, що коваріаційна матриця вектора оцінок параметрів моделей містить досить важливу інформацію про їхню якість.

Незміщена оцінка дисперсії залишків розраховується так:

u2

u'u

,

(4.17)

n m

 

 

 

де n — кількість спостережень; m — кількість змінних у моделі.

 

 

 

, то добуток векто-

Оскільки вектор залишків u Y Y

Y XA

рів u'u можна записати так:

 

 

 

 

 

u' u Y

 

 

Y .

 

 

Y A X

 

Звідси маємо альтернативну форму запису дисперсії залишків:

u2

 

 

Y Y A X Y .

 

 

 

n m

 

 

Позначимо (k, j)-й елемент матриці ( X X ) 1

символом ckj , то-

 

 

 

 

обчислюєть-

ді j-й елемент по головній діагоналі матриці cov ( A)

ся за формулою:

 

 

) .

 

(4.18)

a2j

u2 c jj ( j

1, m

 

111

Коваріації

a

a

,

що містяться за межами головної діагоналі,

відповідно такі:

j k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у€

 

 

2 c

 

 

.

 

 

 

(4.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aa

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j k

 

 

 

 

 

 

 

 

 

 

 

Приклад 4.3. Для економетричної моделі Y=8,8+0,2Х1+

+6,97Х2+u, (приклад

 

4.1)

 

 

обчислимо

коваріаційну матрицю

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cov (A) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отже, маємо:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8,8

 

 

 

 

 

 

1

 

 

 

 

0,314

 

 

 

0,0017

0,0446

 

 

0,2

 

 

( X

 

X )

 

 

 

 

0,00017

 

 

 

 

0,00003

0,00012

 

A

 

;

 

 

 

 

 

 

 

 

 

;

 

 

6,97

 

 

 

 

 

 

 

 

 

 

 

 

0,0446

 

 

 

0,00012

0,0165

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1496

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

Y

 

520090

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n=16; m=2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8367,9

 

 

 

 

 

 

 

Обчислимо незміщену оцінку дисперсії залиш-

Розв’язання. 1.

ків €2 , скориставшись співвідношенням:

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u2

 

 

 

 

 

Y ;

 

 

 

 

 

 

 

 

 

 

 

 

Y Y A X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n m

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

342

... 1972 176394 ;

 

 

 

Y Y yi2 222

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1496

 

 

 

 

 

 

 

 

 

 

8,8

 

 

0,2

6,97

 

 

520090

 

 

 

 

 

 

 

A X Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8367,9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13156 104018 58324,0 175498,0 ;

 

 

 

 

u2 176394 175498

 

896 68,92 .

 

 

 

 

 

 

 

 

 

 

 

 

16 3

 

 

 

 

13

 

 

 

2. Визначимо дисперсії оцінок a2j

:

 

 

 

 

 

 

 

 

var (a€ ) €2

 

2c

 

= 68,92 0,314 = 21,64;

 

 

 

 

 

 

1

 

 

a

 

 

 

 

 

u

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

var (a€ )

2

2c

22

= 68,92 0,00003 = 0,00207;

 

 

 

 

2

 

 

a

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

112

var (a3 ) €2a3 u2c33 = 68,92 0,0165 = 1,137.

3. Обчислимо коваріації відповідних оцінок параметрів:

aa

 

2c

= 68,92 (–0,00017) = –0,0118;

 

 

 

u 12

 

 

1

 

2

 

 

 

у€

 

 

 

 

2c

= 68,92 (–0,0446) = –3,0738;

 

aa

 

u 13

 

 

 

1

 

3

 

 

aa

 

u2c23

= 68,92 (–0,00012) = –0,00827.

 

2

3

 

 

 

 

Знак «мінус» перед оцінками коваріацій у€ajak указує на те, що

зі збільшенням однієї оцінки параметрів інша зменшується в середньому і навпаки.

Отже, дістанемо дисперсійно-коваріаційну матрицю

 

21,64

0,0118

3,0738

 

 

0,0118

0,00207

0,00827

 

cov ( A)

 

.

 

 

3,0738

0,00827

1,137

 

4. Запишемо стандартні похибки оцінок параметрів моделі:

Saj cov (aj ) ;

Sa1 21,64 4,65 ;

Sa2 0,00207 0,0455 0,046 ;

Sa3 1,137 1,066 .

Стандартні похибки характеризують середні лінійні коливання оцінок параметрів моделі навколо свого математичного сподівання. Чим менші ці похибки, тим стійкіші оцінки параметрів моделі. Остаточні висновки стосовно стійкості оцінок можна зробити лише тоді, коли порівняти її з абсолютними значеннями оцінок параметрів моделі.

Порівняємо кожну стандартну похибку Saj з відповідним число-

вим значенням оцінки параметра, тобто знайдемо відношення Saj : aj

113

 

Sa1

100

 

4,65

100

52,8 % ;

 

a

8,8

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

Sa2

100

 

 

0,046

100

23 % ;

 

a

 

0,2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Sa3

100

 

1,066

100

15 % .

 

a

6,97

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

Отже, стандартні похибки оцінок параметрів щодо рівня самих оцінок становлять відповідно 52,8%, 23% і 15%, а це свідчить про зміщеність оцінок.

Це означає, що залишки можуть мати систематичну складову, яка зумовлюється неточною специфікацією моделі. Наприклад, не всі основні чинники, що впливають на тижневі витрати, пов’язані з харчуванням (скажімо, ціни на продукти харчування) внесено до моделі.

4.7. Прогноз залежної змінної

Економетричне моделювання зв’язку між економічними показниками завжди складаєтьмя з трьох етапів:

1)побудови економетричної моделі;

2)перевірки статистичної значущості моделі та оцінювання її параметрів;

3)прогнозування на основі моделі.

Використаємо модель (4.1) для знаходження прогнозного значення y0, яке відповідатиме очікуваним значенням матриці незалежних змінних X0.

Розглянемо спочатку точковий прогноз і припустимо, що ми визначили його як деяку лінійну функцію від yi:

n

(4.20)

y0

ci yi ,

 

i 1

 

 

 

); ci

 

де і — номер спостереження ( i 1, n

— вагові коефіцієнти

значень yi (їх потрібно вибрати так, щоб значення y0 було най-

кращим лінійним незміщеним прогнозом).

 

Оскільки Y XA u, то незміщена точкова оцінка прогнозу

M [ y0 ( X 0 )] X 0 A,

(4.21)

114

де Х0 — матриця очікуваних значень пояснювальних змінних. Задаючи X0, підставимо значення цього вектора в побудовану

економетричну модель

 

(4.22)

y0

 

X 0 A.

Щоб дістати інтервальний прогноз, необхідно розрахувати середню похибку прогнозу.

Вона зростає з віддаленням прогнозного значення x0 j від від-

повідного середнього значення вибірки. Розрахуємо спочатку дисперсію прогнозу.

Уматричномувиглядідисперсіяпохибкипрогнозуподаєтьсятак:

пр2

u2 X 0 ( X X ) 1 X 0 .

(4.23)

Середньоквадратична похибка прогнозу

s

y€ пр

X

(X X ) 1 X

0

.

(4.24)

 

u

0

 

 

 

Довірчий інтервал для прогнозних значень

yt

 

X ( X

X ) 1 X M ( y

( X

 

))

(4.25)

0 u

0

 

 

 

0

 

0

 

0

 

 

 

 

 

 

 

 

y

t

u

X

( X X ) 1

X

0

,

 

 

 

 

0

 

 

0

 

 

 

 

 

 

де t — критичне значення t-критерію при n–m ступенях свободи і рівні значущості .

Зауважимо, що y0

 

є точковою оцінкою як математично-

X 0 A

го сподівання прогнозного значення

y0 , так і його індивідуально-

го значення y0

для відповідних незалежних змінних X 0 , що ле-

жить за межами базового періоду.

 

 

 

 

 

 

Для визначення інтервального прогнозу індивідуального зна-

чення y0 необхідно знайти відповідну стандартну похибку:

2

 

2

2

2

2

X (X

X ) 1 X

0

 

пр(i)

 

u

np

 

u

 

np

0

 

 

 

 

 

 

2

(1

X

 

(X

X ) 1 X

0

).

 

 

 

 

 

u

 

 

0

 

 

 

 

 

 

Отже, інтервальний прогноз індивідуального значення визначається як

y0 t пр i y0 y0 t у€пр i

або

115

yt

1 X ( X X ) 1 X

y

(4.26)

0

u

 

 

0

 

0 0

 

 

 

 

y

t

 

у€

1

X

( X X ) 1 X

0

.

 

0

 

u

 

0

 

 

 

 

Приклад 4.4. Необхідно розрахувати для економетричної моделі (приклад 4.1) точковий та інтервальний прогнози математичого сподівання та індивідуального значення залежної змінної, коли для прогнозного періоду заданий вектор

1 . X 0 5006

Розв’язання. 1.Визначимо точкові прогнозні значення залежної змінної, коли

 

 

1

 

 

 

X0

 

500

:то y0

8,8 0,2x1 6,97x2

8,8 0,2 500 6,97 6 150,62 .

 

 

6

 

 

 

 

 

 

 

 

Отже, y0 можна інтерпретувати як точкову оцінку прогнозного значення математичного сподівання та індивідуального значення витрат на харчування, коли відомі загальні витрати x1=500 і розмір сім’ї становить x2=6.

2. Визначаємо прогнозний інтервал математичного сподівання

M y0 :

 

 

2

 

2

X

( X

X )

1

X

 

 

 

 

 

 

 

 

 

 

 

 

n

u

 

0

X var ( A)X

0

 

 

 

 

 

 

 

 

 

0

р

 

 

 

 

0

 

 

 

 

 

 

 

 

 

1 500

 

 

 

 

21,64

 

0,0118

 

3,0738

 

 

1

 

 

 

 

 

6

0,0118

 

0,00207

 

0,00827

 

 

500

 

 

 

 

 

 

3,0738

0,00827

 

1,137

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2,703

 

 

0,9736

0,3868

 

1

 

 

 

 

 

 

 

 

 

 

 

 

500

481,777.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Стандартна

похибка

прогнозу математичного

сподівання

M y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

пр

 

пр2

 

 

481,777

21,95 .

 

 

 

 

 

 

 

116

3. Знайдемо інтервальний прогноз для M y0 . При цьому не-

хай =0,05 і n–m = 13; тоді t0,05=2,160.

Отже,

y0 t у€пр M ( y0 ) y0 t пр

і

150,62 – 2,160 21,95 M ( y0 ) 150,62 + 2,160 21,95; 150,62 – 47,412 M ( y0 ) 150,62 + 47,412;

103,208 M ( y0 ) 198,032.

4. Обчислимо дисперсію і стандартну похибку прогнозу індивідуального значення y0 :

у2пр(i) 2n u2 481,777 68,92 550,697 .

Стандартнапохибка прогнозу індивідуального значення y0 така:

пр(i) 2пр(i) 550,697 23,467 .

5.Визначимо інтервальний прогноз індивідуального значення y0:

y0 t упр(i) y0 y0 t пр(i) ;

150,62 – 2,160 23,467 y0 150,62 + 2,160 23,467; 150,62 – 50,689 y0 150,62 + 50,689;

99,931 y0 201,309.

Значення t знаходимо в таблиці при = 0,05 і ступені свободи = 13. У такому разі t0,05=2,160.

Отже, з імовірністю р = 0,95 ( = 0,05) прогноз математичного сподівання М(y0) потрапляє в інтервал [103,208; 198,032], а прогноз індивідуального значення — в інтервал [99,931; 201,309].

Можна також сказати, що з імовірністю р = 0,95 знайдені прогнози покривають М(y0) і y0, коли взяти досить велику кількість вибірок і для кожної з них обчислити інтервальні прогнози.

Економічна інтерпретація: якщо у прогнозному періоді загальні витрати мають рівень 500 одиниць, а сім’я складається з шести осіб, то середні витрати на харчування потрапляють в інтервал

103,208 M ( y0 ) 198,032.

117

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