cgiirbis_64 (6)
.pdf
пояснюючих змінних (F-критерій); в) кожної пари пояснюючих змінних (t-критерій).
Усі ці критерії при порівнянні з їх критичними значеннями дають змогу зробити конкретні висновки щодо наявності чи відсутності мультиколінеарності пояснюючих змінних.
Опишемо алгоритм Фаррара—Глобера.
Крок 1. Нормалізація змінних.
Нехай Х1, Х2, Х3, ..., Хm — вектори пояснювальних змінних економетричної моделі. Елементи стандартизованих векторів обчислимо за формулою:
|
X ik |
|
|
|
|
X ik |
|
k |
, |
X ik* |
X |
k |
, або |
X ik* |
X |
||||
|
|
|
|
|
|
||||
|
n 2 X k |
|
X k |
||||||
де п — число спостережень i 1, n ; т — число пояснювальних змінних k 1, m ; X k — середня арифметична k-ї пояснювальної змінної; 2X k — дисперсія k-ї пояснювальної змінної.
Крок 2. Знаходження кореляційної матриці першого порядку
r X * X * ,
де X * — матриця нормалізованих пояснювальних змінних; X * — матриця, транспонована до матриці X * .
Крок 3. Визначення критерію «хі-квадрат»:
2 n 1 1 2m 5 ln r ,6
де |r| — визначник кореляційної матриці r.
Значення цього критерію порівнюється з табличним за таких умов: 12 m m 1 ступенів свободи і рівень значущості . Якщо
факт2 2табл , то в масиві незалежних змінних мультиколінеарність відсутня.
Крок 4. Визначення матриці, оберненої до r:
Cr 1 X * X * 1 .
29
Крок 5. Розрахунок F-критеріїв:
Fk ckk 1 nm m1 ,
де ckk — діагональні елементи матриці С. Фактичні значення критеріїв Fk порівнюються з табличними, коли маємо т – 1 і п – т
ступенів свободи і рівень значущості α. Якщо Fk факт Fтабл , від-
повідна k-та пояснююча змінна мультиколінеарна з іншими. Коефіцієнт детермінації для кожної змінної обчислюється так:
R2 |
1 |
1 |
. |
|
|||
X k |
|
ckk |
|
|
|
||
Крок 6. Знаходження частинних коефіцієнтів кореляції:
rkj |
ckj |
, |
|
ckk * c jj |
|||
|
|
де сkj — елемент матриці С, що міститься в k-му рядку і j-му стовпці ( k 1, m , j l, m , ckk і cjj — діагональніелементи матриціС).
Крок 7. Розрахунок t-критеріїв:
tkj rkj n m .
Фактичні значення критеріїв tkj порівнюються з табличними, коли маємо n–m ступенів свободи і рівень значущості α. Якщо
tkj факт tтабл , між пояснюючими змінними Xk і Xj існує мультиколінеарність.
Дослідження наявності мультиколінеарності згідно з алгоритмом Фаррара—Глобера
Розглянемо застосування алгоритму Фаррара—Глобера для розв’язування конкретної задачі.
Приклад. Оскільки між пояснюючими змінними можлива тісна залежність (мультиколінеарність), метод 1МНК може дати зміщені кількісні характеристики взаємозв’язку. Тому дослідимо наявність мультиколінеарності між пояснюючими змінними, скориставшись алгоритмом Фаррара—Глобера.
30
Якщо виявлено мультиколінеарність, потрібно перетворити певним чином вихідну інформацію, щоб уникнути мультиколінеарності, а далі знову оцінити параметри моделі методом 1МНК.
Для дослідження мультиколінеарності наведемо статистичну сукупність спостережень (табл. 1).
|
|
|
|
Таблиця 1 |
|
|
|
|
|
Питомі інвестиції Х3, |
|
Місяць |
Прибуток на |
Фондовіддача |
Продуктивність |
|
|
місяць Y, грн. |
Х1, грн. |
праці Х2, грн. |
грн. |
|
|
|
|
||||
1-й |
40 |
12 |
5 |
15 |
|
|
|
|
|
|
|
2-й |
45 |
17 |
7 |
18 |
|
|
|
|
|
|
|
3-й |
40 |
13 |
6 |
16 |
|
|
|
|
|
|
|
4-й |
43 |
14 |
7 |
17 |
|
|
|
|
|
|
|
5-й |
48 |
16 |
6 |
20 |
|
|
|
|
|
|
|
6-й |
39 |
15 |
5 |
15 |
|
|
|
|
|
|
|
7-й |
42 |
14 |
6 |
16 |
|
|
|
|
|
|
|
8-й |
45 |
17 |
9 |
18 |
|
|
|
|
|
|
|
9-й |
38 |
12 |
5 |
19 |
|
|
|
|
|
|
|
10-й |
48 |
18 |
10 |
20 |
|
|
|
|
|
|
|
11-й |
50 |
20 |
11 |
22 |
|
|
|
|
|
|
|
12-й |
48 |
17 |
10 |
21 |
|
|
|
|
|
|
|
13-й |
49 |
18 |
12 |
21 |
|
|
|
|
|
|
|
14-й |
45 |
19 |
8 |
20 |
|
|
|
|
|
|
|
15-й |
49 |
20 |
9 |
22 |
|
|
|
|
|
|
|
16-й |
52 |
22 |
14 |
23 |
|
|
|
|
|
|
|
17-й |
54 |
24 |
15 |
24 |
|
|
|
|
|
|
|
18-й |
51 |
21 |
13 |
20 |
|
|
|
|
|
|
|
19-й |
55 |
25 |
16 |
24 |
|
|
|
|
|
|
|
20-й |
56 |
27 |
18 |
25 |
|
|
|
|
|
|
|
Необхідно:
1. Обчислити середнє значення та стандартні відхилення.
31
2.Нормалізувати пояснюючі змінні.
3.Знайти кореляційнуматрицюr.
4.Визначити визначник матриці r (det r).
5.Обчислити критерій 2 .
6.Визначити F-критерії.
7.Визначити частинні коефіцієнти кореляції r.
8.Визначити t-критерії.
9.Зробити висновки щодо мультиколінеарності.
Розв’язання
Для розв’язання використаємо стандартні функції Excel. 1.Обчислимо середні значення та стандартні відхилення пояс-
нюючих змінних Х1, Х2, Х3.
Для цього в мастер функції знайдемо категорію «статистичні» і функції «СРЗНАЧ» та «СТАНДВІДХИЛ».
Х1 |
Х2 |
Х3 |
|
|
|
18,05 |
9,6 |
19,8 |
|
|
|
4,260899 |
3,965642 |
3,036619 |
|
|
|
2.Нормалізуємо пояснюючі змінні.
Серед статистичних функцій знайдемо функцію «НОРМАЛІЗАЦІЯ» та нормалізуємо Х1, Х2, Х3:
32
|
* |
* |
* |
|
|
X1 |
X 2 |
X 3 |
|
|
1,41989 |
1,15996 |
1,58071 |
|
|
0,24643 |
0,65563 |
0,59276 |
|
|
1,1852 |
0,9078 |
1,25139 |
|
|
0,9505 |
0,65563 |
0,92208 |
|
|
|
|||
|
0,48112 |
0,9078 |
0,065963 |
|
|
0,71581 |
1,15996 |
1,58071 |
|
|
0,9505 |
0,9078 |
1,25139 |
|
|
0,24643 |
0,1513 |
0,59276 |
|
|
1,41989 |
1,15996 |
0,26345 |
|
X * |
0,01173 |
0,100866 |
0,065863 |
. |
|
0,45765 |
0,353032 |
0,72449 |
|
|
|
|||
|
0,24643 |
0,100866 |
0,395176 |
|
|
0,01173 |
0,605198 |
0,395176 |
|
|
0,222958 |
0,40347 |
0,065863 |
|
|
0,45765 |
0,1513 |
0,72449 |
|
|
0,927034 |
1,10953 |
1,053804 |
|
|
1,396419 |
1,361696 |
1,383117 |
|
|
0,692342 |
0,857364 |
0,065863 |
|
|
1,631111 |
1,613862 |
1,383117 |
|
|
|
|||
|
2,100496 |
2,118194 |
1,712431 |
|
Транспонуємо матрицю Х* (нормалізовану) в матрицю X *
|
|
1,4199 |
0,2464 |
1,1852 |
0,9505 |
0,4811 |
0,7158 |
0,9505 |
|
|
|
|
1,16 |
0,6556 |
0,9078 |
0,6556 |
0,9078 |
1,16 |
0,9078 |
X * |
|
||||||||
|
|
|
1,5807 |
0,5928 |
1,2514 |
0,9221 |
0,06586 |
1,5807 |
1,2514 |
|
|
|
|
|
|
|
|
Продовження матриці X * |
|
|
|
|
|
||
0,2464 |
1,4199 |
0,0117 |
0,45765 |
0,2464 |
0,0117 |
0,22296 |
0,45765 |
0,1513 |
1,16 |
0,10087 |
0,35303 |
0,10087 |
0,6052 |
0,4035 |
0,1513 |
0,5928 |
0,2635 |
0,06586 |
0,72449 |
0,39518 |
0,39518 |
0,06586 |
0,72449 |
|
|
|
|
|
|
|
|
|
|
Завкінчення матриці X * |
|
|
|
|
|
|
|
|
|
0,92703 |
|
1,39642 |
0,69234 |
1,63111 |
2,1005 |
|
|||
1,10953 |
|
1,3617 |
0,85736 |
1,61386 |
2,11819 |
. |
|||
1,0538 |
|
1,38312 |
0,06586 |
1,38312 |
1,71243 |
|
|
||
|
|
|
|||||||
Перемножимо матриці X * |
та Х*: |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
19 |
|
17,8964552 |
16,9413894 |
|
|||
|
|
|
|
|
19 |
16,6415575 |
|
||
X * X * |
17,8964552 |
|
. |
||||||
|
16,9413894 |
16,6415575 |
|
19 |
|
|
|||
3.Знайдемо кореляційну матрицю r.
33
Для знаходження кореляційної матриці r необхідно кожний
|
|
|
|
1 |
|
|
|
|
|
n 1 (у нашому випадку |
|||
елемент матриці X * X * помножити на |
|
|||||
n – 1 = 19, оскільки п = 20): |
|
|
|
|
|
|
|
1 |
0,941918693 |
0,891652074 |
|
||
r |
0,941918693 |
1 |
0,875871449 |
. |
||
|
0,891652074 |
0,875871449 |
1 |
|
||
|
|
|||||
4.Знайдемо визначник матриці r (det r).
Для знаходження det r необхідно серед математичних функцій 6Excel знайти функцію «МОПРЕД». Скориставшись нею, дістанемо:
det r |
0,02182033 |
|
|
Оскільки det r наближається до 0, то в масиві пояснюючих змінних може існувати мультиколінеарність.
Прологарифмуємо визначник матриці r:
ln det r |
–3,824913185 |
|
|
5.Обчислимо критерій 2 |
за формулою: |
||||||||||
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
|
|
n 1 |
6 |
2m 5 ln det r . |
||||||
Маємо: |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
64,38603861 |
|
|
Знайдене |
значення |
2 |
порівнюємо з табличним значенням |
||||||||
kp |
2 7,80 , |
коли маємо |
1 m m 1 3 ступенів свободи та при |
||||||||
|
|
|
|
|
|
2 |
|
|
|
|
|
рівні значущості = 0,05.
34
Оскільки факт2 2крит , то в масиві пояснюючих змінних (про-
дуктивність праці, питомі інвестиції та фондовіддача) існує мультиколінеарність.
6.Обчислимо F-критерії.
Для визначення F-критеріїв необхідно знайти матрицю С, яка є оберненою до матриці r:
|
10,67120467 |
7,375970017 |
3,05460023 |
|
C |
7,375970017 |
9,392918454 |
1,65019013 |
|
|
3,054600232 |
1,65019013 |
5,168995053 |
|
|
|
|||
Безпосередньо F-критерії обчислюються за формулою:
FCkk 1 n m ,
m 1
де сkk — діагональний елемент матриці С.
F1 |
51,57975822 |
|
|
F2 |
44,76223175 |
|
|
F3 |
22,23464028 |
|
|
Обчислені F-критерії порівнюються з табличним значенням (F = 5,29), коли є n – m (16) ступенів свободи та при рівні значущос-
ті = 0,05.
У розглядуваному випадку F1 > Fкрит, F2 > Fкрит, F3 > Fкрит. Це означає, що кожна з пояснюючих змінних мультиколінеарна з іншими.
7.Визначимо частинні коефіцієнти кореляції r.
Частинні коефіцієнти кореляції показують тісноту зв’язку між двома пояснюючими змінними за умови, що інші змінні не впливають на цей зв’язок і обчислюються за формулою:
rkj.s |
ckj |
. |
|
ckk c jj |
|||
|
|
||
r12.3 |
0,736736 |
||
|
|
||
r13.2 |
0,411287 |
||
|
|
||
r23.1 |
0,236827 |
||
|
|
|
|
35
Отже, спираючись на здобуті нами значення окремих (частинних) коефіцієнтів кореляції, можна сказати, що зв’язок між фондовіддачею та продуктивністю праці є тісним, якщо не враховувати вплив питомих інвестицій; зв’язок між фондовіддачею та питомими інвестиціями є слабким, якщо не брати до уваги вплив продуктивності праці. Зв’язок між продуктивністю праці та питомими інвестиціями також є слабким, якщо не враховувати фондовіддачу.
8.Визначимо t-критерії.
Ці критерії застосовуються для визначення мультиколінеарності двох пояснюючих змінних і обчислюються за формулою:
tkj |
rkj.s |
n m |
|||||
|
1 |
r |
2 |
|
|
||
|
|||||||
|
|
|
|
||||
|
|
|
kj.s |
|
|||
|
|
|
|
||||
t12 |
|
|
6,445364 |
||||
|
|
|
|
||||
t13 |
|
|
1,980092 |
||||
|
|
|
|
||||
t23 |
|
|
1,003595 |
||||
|
|
|
|
|
|
|
|
Обчислені t-критерії порівнюються з табличним значенням (t = 1,746), коли маємо n – m (16) ступенів свободи та при рівні значущості = 0,05.
Оскільки t12 > tкрит, то продуктивність праці та фондовіддача є відповідно мультиколінеарними між собою;
t13 > tкрит, тому відповідно продуктивність праці та питомі інвестиції є мультиколінеарними між собою;
t23 < tкрит, тому продуктивність праці та питомі інвестиції не є мультиколінеарними між собою.
Дослідження показали, що мультиколінеарність існує. Отже, для того, щоб можна було застосувати метод 1МНК для оцінювання параметрів моделі за цією інформацією, необхідно звільнитися від мультиколінеарності.
Покажемо, як побудувати економетричну модель методом, 1МНК на основі нормалізованих даних.
36
|
|
|
* |
* |
* |
|
|
|
|
|
|
|
X1 |
X 2 |
X 3 |
|
|
|
|
|
|
|
1,41988804 |
1,15996353 |
1,58070558 |
|
|
1,30559111 |
|
|
|
|
0,24642685 |
0,65563156 |
0,59276459 |
|
|
0,3526049 |
|
|
|
|
|
|
|
||||
|
|
|
1,1851958 |
0,90779755 |
1,25139192 |
|
|
1,30559111 |
|
|
|
|
0,95050356 |
0,65563156 |
0,92207826 |
|
|
0,73379939 |
|
|
|
|
0,48111909 |
0,90779755 |
0,06586273 |
|
|
0,21918683 |
|
|
|
|
0,71581133 |
1,15996353 |
1,58070558 |
|
|
1,49618836 |
|
|
|
|
0,95050356 |
0,90779755 |
1,25139192 |
|
0,92439663 |
||
|
|
|
0,24642685 |
0,15129959 |
0,59276459 |
|
|
0,3526059 |
|
|
|
|
1,41988804 |
1,15996353 |
0,26345093 |
|
1,6876856 |
|
|
X |
* |
|
0,01173461 |
0,100866394 |
0,06586273 |
|
|
0,21918683 |
|
|
|
0,457649864 |
0,353032379 |
0,72449006 |
. |
Y * |
0,60038132 |
. |
|
|
|
|
|
|
|
||||
|
|
|
0,24642685 |
0,100866394 |
0,3951764 |
|
|
0,21918683 |
|
|
|
|
0,01173461 |
0,605198364 |
0,3951764 |
|
|
0,40978407 |
|
|
|
|
0,222957626 |
0,40346558 |
0,06586273 |
|
|
0,3526049 |
|
|
|
|
0,457649864 |
0,151299559 |
0,72449006 |
|
|
0,40978407 |
|
|
|
|
0,927034339 |
1,109530335 |
1,05380372 |
|
|
|
|
|
|
|
1,396418814 |
1,36169632 |
1,38311738 |
|
|
0,9815758 |
|
|
|
|
|
|
1,36277029 |
|
|||
|
|
|
0,692342101 |
0,85736435 |
0,06586273 |
|
|
0,79097856 |
|
|
|
|
1,631111052 |
1,613862305 |
1,38311738 |
|
|
1,55336753 |
|
|
|
|
2,100495528 |
2,118194275 |
1,71243105 |
|
|
1,74396477 |
|
|
|
|
1,631111052 |
1,613862305 |
1,38311738 |
|
|
|
|
|
|
|
2,100495528 |
2,118194275 |
1,71243105 |
|
|
|
|
|
|
|
|
|
|
|
|||
Застосуємо функцію «лінійн» для знаходження всіх оцінок параметрів нової моделі, стандартні похибки цих оцінок, коефіцієнт детермінації, дисперсію залишків і т. ін.
b3 |
b2 |
b1 |
|
0,3183466 |
0,31026984 |
0,39186698 |
0 |
|
|
|
|
0,1529347 |
0,2061594 |
0,21974024 |
#Н/Д |
|
|
|
|
0,9269234 |
0,2932106 |
#Н/Д |
#Н/Д |
|
|
|
|
71,877548 |
17 |
#Н/Д |
#Н/Д |
|
|
|
|
18,538468 |
1,46153177 |
#Н/Д |
#Н/Д |
|
|
|
|
Обчислюємо €* — залежну змінну, підставивши пояснюючі
Y
змінні у відповідну формулу:
€* |
* |
* |
* |
Y |
b1 X1 |
b2 X 2 |
b3 X 3 . |
37
1,41952111
0,48869381
1,1444776
0,86943409
0,44922971
1,14361669
1,05250946
0,33221481
1,00017763
€0,04766457 . Y * 0,51951208
0,0605323
0,308979440,01684629
0,36303308
1,043002731,41001437
0,55828748
1,58022201
2,02547315
Щоб обчислити оцінки параметрів моделі в абсолютному виразі, скористаємося формулами:
€ |
|
|
Y |
|
|
€ |
|
|
|
|
X 2 ; |
€ |
|
|
|
|
X 3 |
; |
|
a1 |
|
b1 |
|
X1 ; a2 |
b2 Y |
|
a3 |
|
b3 |
Y |
|
||||||||
|
|
|
|
a€0 |
Ycp a€1 X cp1 a€2 X cp2 |
a€3 X cp3. |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
a€0 |
|
|
23,30887 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
a€1 |
|
|
0,482526 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
a€2 |
|
|
0,410496 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
a€3 |
|
|
0,550039 |
|
|
|
|
|
|
|
||
Отже, економетрична модель має вигляд:
|
|
€ |
0,410X 2 |
0,550X 3 , |
|
|
€ |
Y 23,31 0,483X1 |
|||
де |
— дохід; Х1 — фондовіддача; Х2 — продуктивність праці; |
||||
Y |
|||||
Х3 |
— питомі інвестиції. |
|
|
||
Завдання для самостійної роботи
Нехай на витрати обігу впливають: обсяг вантажообігу, грн., запаси з вантажообігу, грн., та трудомісткість одиниці вантажообігу, людино-год. Щоб побудувати економетричну модель цієї зале-
38
