Econometrics
.pdft |
|
|
rx j ,u j |
|
n 2 |
. |
(7.8) |
|
|
||||||
|
|
|
|
|
|||
|
|
|
1 r2 |
|
|||
|
|
|
|
|
x j ,u j |
|
Доведено, що ця характеристика має закон розподілу Стьюдента з кількістю ступенів свободи n 2 .
Якщо розраховане значення t-статистики перевищує критичне значення при ступені свободи n–2 та вибраному рівні значущості, то гіпотезу про наявність гетероскедастичності потрібно прийняти. У протилежному випадку вона відхиляється.
Розглянемо застосування тесту рангової кореляції для визначення наявності гетероскедастичності у статистичній інформації.
Приклад 7.5. Необхідно дослідити наявність гетероскедастичностіу статистичній інформації, поданій в табл. 4.2 (приклад 4.2).
Розв’язання
1.Побудуємо прості економетричні моделі, що описують залежність між прибутком та кожним з чинників. Наведемо ці моделі:
1)Y 3,07 0,624X1 ;
2)Y 19,90 0,935X 2 ;
3)Y 28,96 0,658X3 .
2.Розрахуємо залишки за цими економетричними моделями
(табл. 7.7).
Таблиця 7.7
Мі- |
U1 |
U2 |
U3 |
Мі- |
U1 |
U2 |
U3 |
сяць |
|
|
|
сяць |
|
|
|
1 |
– |
– |
– |
11 |
– |
– |
– |
|
2,7 |
1,458 |
0,4 |
|
2,162 |
1,740 |
3,6 |
|
935 |
65 |
783 |
|
116 |
203 |
400 |
|
15 |
|
58 |
|
|
|
65 |
2 |
– |
– |
– |
12 |
– |
– |
4,6 |
|
2,6 |
2,262 |
1,7 |
|
0,284 |
0,478 |
503 |
|
672 |
71 |
687 |
|
983 |
95 |
62 |
|
35 |
|
84 |
|
|
|
|
3 |
– |
2,214 |
1,8 |
13 |
1,964 |
0,521 |
– |
|
0,6 |
783 |
120 |
|
164 |
05 |
0,9 |
|
706 |
|
69 |
|
|
|
304 |
|
48 |
|
|
|
|
|
91 |
4 |
– |
– |
– |
14 |
7,080 |
0,651 |
1,0 |
16
|
2,2 |
|
3,132 |
|
4,0 |
|
|
|
|
|
|||
|
918 |
|
083 |
|
592 |
|
|
09 |
|
|
|
11 |
|
5 |
– |
|
– |
|
– |
|
|
2,1 |
|
2,066 |
|
3,3 |
|
|
655 |
|
77 |
|
753 |
|
|
29 |
|
|
|
82 |
|
6 |
9,7 |
|
10,41 |
|
5,5 |
|
|
047 |
|
072 |
|
731 |
|
|
78 |
|
|
|
3 |
|
|
|
|
|
|
|
|
7 |
– |
|
– |
|
– |
|
|
4,0 |
|
5,805 |
|
5,3 |
|
|
392 |
|
516 |
|
496 |
|
|
49 |
|
|
|
38 |
|
8 |
– |
|
– |
|
– |
|
|
1,7 |
|
2,936 |
|
2,3 |
|
|
901 |
|
143 |
|
753 |
|
|
02 |
|
|
|
82 |
|
9 |
– |
|
– |
|
1,9 |
|
|
1,9 |
|
3,674 |
|
407 |
|
|
129 |
|
89 |
|
89 |
|
|
69 |
|
|
|
|
|
10 |
– |
|
– |
|
0,9 |
|
|
1,4 |
|
1,609 |
|
922 |
|
|
112 |
|
576 |
|
77 |
|
|
63 |
|
|
|
|
|
|
31 |
|
677 |
|
695 |
|
|
|
|
|
|||
|
|
|
|
|
09 |
|
|
|
|
|
|
|
|
15 |
1,090 |
|
0,586 |
|
– |
|
|
444 |
|
364 |
|
1,9 |
|
|
|
|
|
|
047 |
|
|
|
|
|
|
47 |
|
16 |
0,339 |
|
0,455 |
|
– |
|
|
59 |
|
737 |
|
2,5 |
|
|
|
|
|
|
885 |
|
|
|
|
|
|
77 |
|
17 |
1,592 |
|
1,716 |
|
6,3 |
|
|
15 |
|
991 |
|
599 |
|
|
|
|
|
|
35 |
|
|
|
|
|
|
|
|
18 |
0,469 |
|
1,847 |
|
1,7 |
|
|
283 |
|
617 |
|
790 |
|
|
|
|
|
|
82 |
|
|
|
|
|
|
|
|
19 |
– |
|
1,912 |
|
0,8 |
|
|
0,404 |
|
931 |
|
048 |
|
|
437 |
|
|
|
26 |
|
|
|
|
|
|
|
|
20 |
0,346 |
|
4,847 |
|
1,4 |
|
|
416 |
|
617 |
|
886 |
|
|
|
|
|
|
55 |
|
|
|
|
|
|
|
|
Таблиця 7.8
М |
|
|
|
|
|
|
|
|
|
|
М |
|
|
|
|
|
|
|
|
|
|
іс |
R |
|
R |
R |
X2 |
R |
R |
X3 |
RU |
|
іс |
R |
|
RU |
R |
X2 |
R |
R |
X3 |
RU |
|
я |
|
X |
U |
|
U2 |
|
|
я |
|
X |
|
U2 |
|
|
|||||||
ц |
|
|
|
|
|
|
|
|
|
3 |
ц |
|
|
1 |
|
|
|
|
|
|
3 |
ь |
|
|
|
|
|
|
|
|
|
|
ь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
1 |
3 |
|
7 |
3 |
7 |
3 |
7 |
|
1 |
1 |
|
13 |
16 |
12 |
11 |
10 |
|||||
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
|
|
|
|
|
|
|
|
|
2 |
6 |
|
9 |
6 |
1 |
1 |
4 |
|
1 |
1 |
|
15 |
12 |
16 |
17 |
14 |
|||||
|
|
|
|
|
|
|
|
|
|
|
2 |
3 |
|
|
|
|
|
|
|
|
|
3 |
1 |
|
4 |
1 |
2 |
2 |
11 |
1 |
1 |
|
14 |
13 |
20 |
13 |
20 |
||||||
|
|
|
|
|
|
|
|
|
|
|
3 |
6 |
|
|
|
|
|
|
|
|
|
4 |
2 |
|
8 |
2 |
4 |
6 |
5 |
|
1 |
1 |
|
17 |
15 |
18 |
14 |
18 |
|||||
|
|
|
|
|
|
|
|
|
|
|
4 |
2 |
|
|
|
|
|
|
|
|
|
5 |
4 |
|
2 |
4 |
5 |
4 |
16 |
1 |
1 |
|
18 |
14 |
15 |
16 |
3 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
5 |
5 |
|
|
|
|
|
|
|
|
|
6 |
5 |
|
5 |
5 |
11 |
9 |
8 |
|
1 |
1 |
|
19 |
17 |
17 |
15 |
9 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
6 |
7 |
|
|
|
|
|
|
|
|
|
7 |
7 |
|
1 |
7 |
8 |
8 |
2 |
|
1 |
1 |
|
3 |
18 |
13 |
18 |
12 |
17
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
8 |
9 |
1 |
11 |
10 |
7 |
3 |
|
|
|
|
|
|
|
9 |
1 |
1 |
9 |
3 |
12 |
1 |
|
4 |
2 |
|
|
|
|
1 |
1 |
1 |
10 |
19 |
12 |
19 |
0 |
1 |
6 |
|
|
|
|
|
7 |
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
1 |
20 |
20 |
14 |
19 |
6 |
|
8 |
|
9 |
|
|
|
|
|
|
1 |
|
2 |
6 |
19 |
6 |
20 |
7 |
|
9 |
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3. Визначимо ранги для кожної з пояснювальних змінних та залишків. Для цього кожну змінну та кожний із залишків uij
пронумеруємо порядковим номером від одиниці до 20. Потім впорядковуємо кожну пояснювальну змінну та залишки, посортувавши їх від меншого до більшого (Exel, меню «Данные», розділ «Сортировка»). У результаті порядковий номер чисел цих показників зміниться і характеризуватиме ранг його в масиві.
Запишемо ранги пояснювальних змінних та ранг залишків у табл. 7.8.
4. Знайдемо різниці між рангами j-ї змінної та залишками, отриманими на основі цієї змінної, піднесемо їх до квадрата
(табл. 7.9).
|
|
Таблиця 7.9 |
|
|
|
d12 |
d22 |
d32 |
|
|
|
16 |
16 |
16 |
25 |
9 |
9 |
1 |
9 |
81 |
4 |
36 |
1 |
1 |
4 |
144 |
36 |
0 |
1 |
1 |
9 |
100 |
1 |
9 |
36 |
1 |
100 |
36 |
121 |
9 |
121 |
64 |
36 |
81 |
4 |
9 |
1 |
9 |
9 |
9 |
16 |
1 |
49 |
36 |
4 |
16 |
0 |
16 |
169 |
0 |
4 |
36 |
25 |
225 |
36 |
18
|
25 |
|
0 |
|
169 |
|
|
|
|||
|
196 |
|
169 |
|
9 |
|
|
|
|
|
|
d12 |
1120 |
d22 |
1120 |
d32 |
1120 |
i |
|
i |
|
i |
|
5. Визначимо коефіцієнти рангової кореляції (простий мішаний момент):
|
|
20 |
|
|
|
|
|
|
6 d 2 |
|
6 582 |
|
|
rX u |
1 |
i 1 i |
1 |
0,562 ; |
||
|
203 20 |
|||||
1 1 |
|
n3 n |
|
|
||
rX u |
1 |
6 674 |
0,493 ; |
|
||
203 20 |
|
|||||
2 2 |
|
|
|
|
||
r |
1 |
6 1120 |
0,158 . |
|
||
203 20 |
|
|||||
X3u3 |
|
|
|
|
6. Розрахуємо t-критерії для визначення статистичної значущості коефіцієнтів рангової кореляції Спірмена за формулою:
|
|
|
t j |
rX ju j |
|
n 2 |
; |
|
||||
|
|
|
1 r |
2 |
|
|
||||||
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
X ju j |
|
|
|||
t |
0,562 18 |
2,886 ; |
t |
2 |
|
|
0,493 18 |
2,886 ; |
||||
|
|
|||||||||||
1 |
1 |
0,5622 |
|
|
|
|
|
|
|
1 0,4932 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
t3 |
|
0,158 18 |
|
0,678 . |
|
|||||
|
|
1 0,1582 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
Критичне значення t-критерію для =0,05 і ступенів свободи n–2=18 дорівнює: t0,05крит 2,10 .
Порівнюючи розраховані t-критерії з критичним значенням, робимо висновок, що коефіцієнти рангової кореляції rX1u1 та rX2u2 є
статистично значущими. Звідси можна зробити висновок, що пояснювальнізмінні Х1 та Х2 можуть викликати гетероскедастичність.
7.3.6. Тест Парка. Для визначення наявності гетероскедастичності у статистичній інформації Р. Парк запропонував параметричний тест, в основі якого лежить визначення кількісної залежності між дисперсією пояснювальної змінної, яка може
19
викликати гетероскедастичність, та значеннями цієї змінної за функцією:
2x |
ij |
u2 xij j uvj j , |
(7.9) |
|
|
|
де 2xij —дисперсія j-ї пояснювальної змінної;
u2 —дисперсія залишків;
xij —i-те значення j-ї пояснювальної змінної;
uij —стохастичні залишки для i-го спостереження j-ї поясню-
вальної змінної.
Прологарифмувавши вираз, дістанемо:
ln 2xij |
ln u2 j ln xij |
v j , |
де вираз vj ln(uj) замінено на vj.
Оскільки 2xij , як правило, невідомі, то Парк запропонував за-
мінити їх квадратами залишків. Розглянемо алгоритм тесту Парка.
Крок 1. Побудова рівняння регресії між змінною (Y) та відповідною пояснювальною змінною (Xj)
€ |
a€1 X j . |
|
|||
Y a€0 |
|
||||
Крок 2. Визначення залишків |
u Y Y |
та піднесення їх до |
|||
квадрата u2 . |
|
|
|
|
|
Крок 3. Знаходження логарифмів |
u2 |
та x . |
|||
|
|
i |
|
ij |
|
Крок 4. Побудова рівняння регресії |
|
|
|
||
ln u2 |
x |
v |
|
; |
(7.10) |
i |
j ij |
ij |
|
|
де ln u2 .
Якщо досліджуються кілька пояснювальних змінних, то таке рівняння регресії будується для кожної з них.
Крок 5. Перевірка статичної значущості оцінки параметра j
|
|
|
|
|
j |
|
|
на основі t-критерію |
t |
j |
|
|
|
. Якщо фактичне значення t- |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
j |
|
20
статистики більше за критичне значення t-критерію, для вибраного рівня значущості і ступеня свободи n–2, то є статистично
значущий. А це означає, що гіпотезу про наявність гетероскедастичності, що викликається Xj, не можна відхиляти. Зауважимо, що тест Парка не вичерпує проблему визначення гетероскедастичності, бо вона може існувати навіть якщо врахувати ті пояснювальні змінні, які в даному випадку не досліджуються, але вони впливають на залежну змінну. Тому поряд із цим тестом доцільно використовувати інші. Зокрема, тест Глейзера доповнює тест Парка аналізом гетероскедастичності, що викликається тими змінними, які не включено до економетричної моделі. Він може виявити наявність мішаної гетероскедастичності.
Приклад 7.6. Необхідно перевірити наявність гетероскедастичності у статистичній інформації, наведеній у табл. 4.2 (приклад 7.2), використавши параметричний тест Парка.
Розв’язання
1.Побудуємо прості економетричні моделі, які описують зв’язок залежної змінної з кожною з пояснювальних.
2.Визначимо залишки за кожною з моделей та піднесемо їх до квадрата:
u2 |
u2 |
u2 |
1 |
2 |
3 |
7,803728 |
2,12766 |
0,228826 |
7,114145 |
5,119856 |
3,128598 |
0,449769 |
4,905264 |
3,283594 |
5,252388 |
9,809945 |
16,47719 |
4,689516 |
4,271538 |
11,3932 |
94,18272 |
108,3832 |
31,05978 |
16,31553 |
33,70402 |
28,61863 |
3,204467 |
8,620936 |
5,642439 |
3,659451 |
13,50481 |
3,766662 |
1,991663 |
2,590736 |
0,984613 |
4,674746 |
3,028306 |
13,25007 |
21
0,081215 |
0,229393 |
21,62587 |
3,85794 |
0,271494 |
0,865814 |
50,226 |
0,424683 |
1,143849 |
1,189067 |
0,343823 |
3,628063 |
0,115322 |
0,207696 |
6,700729 |
2,534942 |
2,948057 |
40,44878 |
0,220227 |
4,41369 |
3,165133 |
0,163569 |
3,659304 |
0,647745 |
0,120004 |
23,49939 |
2,216094 |
3. Прологарифмуємо квадрати залишків та пояснювальні змінні і подамо їх у табл 7.10.
Таблиця 7.10
№ |
ln u2 |
ln u2 |
ln u2 |
ln X |
1 |
|
ln X |
2 |
|
ln X |
3 |
|
з/п |
1 |
2 |
3 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
1 |
2,054602 |
0,755023 |
–1,474793 |
4,127134 |
3,091042 |
4,644391 |
||||||
|
|
|
|
|
|
|
||||||
2 |
1,962085 |
1,633126 |
1,140585 |
4,174387 |
3,218876 |
4,691348 |
||||||
|
|
|
|
|
|
|
||||||
3 |
–0,79902 |
1,590309 |
1,188939 |
4,043051 |
2,833213 |
4,59512 |
||||||
|
|
|
|
|
|
|
||||||
4 |
1,658683 |
2,283397 |
2,801977 |
4,189655 |
3,295837 |
4,736198 |
||||||
|
|
|
|
|
|
|
||||||
5 |
1,545329 |
1,451974 |
2,433017 |
4,234107 |
3,332205 |
4,75359 |
||||||
|
|
|
|
|
|
|
||||||
6 |
4,545237 |
4,685673 |
3,435914 |
4,060443 |
2,995732 |
4,70048 |
||||||
|
|
|
|
|
|
|
||||||
7 |
2,792118 |
3,517617 |
3,354058 |
4,276666 |
3,465736 |
4,779123 |
||||||
|
|
|
|
|
|
|
||||||
8 |
1,164546 |
2,154194 |
1,730316 |
4,248495 |
3,401197 |
4,75359 |
||||||
|
|
|
|
|
|
|
||||||
9 |
1,297313 |
2,603046 |
1,326189 |
4,317488 |
3,526361 |
4,736198 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
22
10 |
0,68897 |
0,951942 |
–0,015506 |
4,369448 |
3,555348 |
4,787492 |
|
|
|
|
|
|
|
11 |
1,542175 |
1,108004 |
2,584003 |
4,343805 |
3,496508 |
4,820282 |
|
|
|
|
|
|
|
12 |
–2,510652 |
–1,47232 |
3,07389 |
4,406719 |
3,610918 |
4,779123 |
|
|
|
|
|
|
|
13 |
1,350133 |
–1,303817 |
–0,144085 |
4,382027 |
3,610918 |
4,859812 |
|
|
|
|
|
|
|
14 |
3,916533 |
–0,856412 |
0,134399 |
4,317488 |
3,663562 |
4,859812 |
|
|
|
|
|
|
|
15 |
0,173169 |
–1,06763 |
1,288699 |
4,418841 |
3,637586 |
4,882802 |
|
|
|
|
|
|
|
16 |
–2,16003 |
–1,571678 |
1,902216 |
4,394449 |
3,583519 |
4,867534 |
|
|
|
|
|
|
|
17 |
0,930171 |
1,081146 |
3,700036 |
4,465908 |
3,688879 |
4,820282 |
|
|
|
|
|
|
|
18 |
–1,513097 |
1,227794 |
1,152192 |
4,521789 |
3,73767 |
4,89784 |
|
|
|
|
|
|
|
19 |
–1,810519 |
1,297273 |
-0,434259 |
4,553877 |
3,7612 |
4,919981 |
|
|
|
|
|
|
|
20 |
–2,120228 |
3,156975 |
0,795746 |
4,574711 |
3,73767 |
4,934474 |
|
|
|
|
|
|
|
Побудуємо три економетричні моделі такого виду:
ln ui2 i j Xij ; де i ln u2 .
1)ln u12 35,313 8,002 X1 ;
2)ln u22 10,195 2,609 X 2 ;
3)ln u32 9,534 1,677 X 3 .
4. Розрахуємо t-критерії для перевірки статистичної значущості кожної з оцінок параметрів цих моделей:
а) для першої моделі:
t 3,34; |
t 1 |
3,27; |
б) для другої моделі: |
|
|
t 2,05; |
t 2 |
1,83; |
в) для третьої моделі: |
|
|
23
t 0,248; |
t 3 |
3,70. |
5. Порівнюючи кожний із фактичних t-критеріїв із табличним, робимо висновки, що гетероскедастичність в досліджуваній статистичній інформації може існувати за рахунок кожного з чинни-
ків, бо t j tкрит .
tкрит = 2,10 |
для |
0,05 |
і ступенів свободи n 2 18; |
tкрит = 1,77 |
для |
0,10 |
і ступенів свободи n 2 18. |
Зауважимо, що висновки про наявність гетероскедастичності на основі декількох тестів не завжди збігаються. Це пов’язано з тим, що кожний із них побудований за певних вихідних припущень, тобто має свої переваги та недоліки. На наш погляд, тест Глейзера є найбільш простий і одночасно найбільш точний. Хоч і в даному випадку може бути присутня помилка специфікації рівняння регресії залишків від пояснювальної змінної.
7.4. Визначення матриці S
Щоб оцінити параметри моделі на основі інформації, коли існує гетероскедастичність і коли дисперсії залишків визначаються
M (uu ) u2 S, потрібно визначити матрицю S.
Спинимося на визначенні матриці S.
Оскільки явище гетероскедастичності пов’язане лише з тим, що змінюються дисперсії залишків, а коваріація між ними відсутня, то матриця S має бути діагональною і додатно визначеною, а саме:
24
|
1 |
0 |
0 ... |
0 |
|
|
|||
|
|
|
|
||||||
1 |
|
||||||||
|
1 |
|
|
|
|
|
|
||
0 |
0 ... |
0 |
|
|
|||||
|
2 |
|
|||||||
S |
|
|
1 |
|
|
. |
(7.11) |
||
|
|
|
|
|
|||||
|
0 |
0 |
|
|
... |
0 |
|
|
|
|
3 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
... ... |
... ... |
... |
|
||||||
|
0 |
0 |
0 ... |
1 |
|
|
|||
|
|
|
|
||||||
|
|
||||||||
|
|
|
|
|
|
|
n |
|
Аби пояснити, чому саме такий вигляд має ця матриця, потрібно ще раз наголосити: за наявності гетероскедастичності для певних вихідних даних одна (або кілька) пояснювальних змінних можуть різко змінюватися від одного спостереження до іншого, тоді як залежна змінна має такі самі коливання, як і для попередніх спостережень.
Це означає, що дисперсія залишків, яка змінюватиметься від одного спостереження до іншого (чи для групи спостережень), може бути пропорційною до значення пояснювальної змінної Xj (або до її квадрата), яка зумовлює гетероскедастичність, або пропорційною до квадрата залишків.
Звідси в матриці S значення i можна обчислити, користуючись такими гіпотезами:
а) M (uu ) u2 xij , тобто дисперсія залишків пропорційна до зміни пояснювальної змінної xij ;
б) M (uu ) u2 xij2 , тобто зміна дисперсії пропорційна до зміни
квадрата пояснювальної змінної ( x2 ); |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
ij |
|
2 |
|
ui |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
в) M (uu ) u |
|
|
, тобто дисперсія залишків пропорційна до |
|||||||||
зміни квадрата |
|
розрахункових залишків за модулем. |
||||||||||
Для першої гіпотези: yi |
|
1 |
; |
|||||||||
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
xij |
||
для другої гіпотези: |
yi |
|
|
1 |
; |
|
||||||
|
|
|||||||||||
|
|
|
|
|
|
|
|
xij2 |
25