пр 2 Стешенко, (екон)
.doc
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ БІОРЕСУРСІВ І ПРИРОДОКОРИСТУВАННЯ УКРАЇНИ
Кафедра економетрики
Практичне заняття 2
Виконала студентка 3 курсу 1 групи
Економічного факультету
Стешенко Тетяна Олександрівна
Київ - 2023
Практична робота 2.
Завдання:
Дослідити серед наведених чинників наявність мультиколінеарності. Нехай:
У – середньомісячна зарплата, грн.; на неї впливають наступні незалежні змінні
Х1 – продуктивність праці, люд.-дні;
Х2 – фондомісткість, млн. грн.;
Х3 – коефіцієнт плинності робочої сили, %.
Для зручності всі необхідні розрахунки оформлюватимемо у таблиці
Варіант 1.
Вихідні дані та розрахунок стандартизованих змінних та коефіцієнтів кореляції
|
Таблиця 1 |
||||||||||||||||
№ п/п |
Х1 |
Х2 |
Х3 |
Х1-Х1(сер) |
Х2-Х2(сер) |
Х3-Х3(сер) |
(Х1-Х1(сер))2 |
(Х2-Х2(сер))2 |
(Х3-Х3(сер))2 |
Х1* |
Х2* |
Х3* |
Х1**Х2* |
Х1**Х3* |
Х2**Х3* |
|
|
1 |
7,50 |
11,80 |
9,70 |
-14,44 |
-15,73 |
-10,38 |
208,51 |
247,43 |
107,74 |
-0,03051 |
-0,01299 |
-0,01546 |
0,00040 |
0,00047 |
0,00020 |
|
|
2 |
9,40 |
10,80 |
9,40 |
-12,54 |
-16,73 |
-10,68 |
157,25 |
279,89 |
114,06 |
-0,01637 |
-0,01774 |
-0,02062 |
0,00029 |
0,00034 |
0,00037 |
|
|
3 |
11,40 |
11,90 |
9,10 |
-10,54 |
-15,63 |
-10,98 |
111,09 |
244,30 |
120,56 |
-0,00149 |
-0,01251 |
-0,02577 |
0,00002 |
0,00004 |
0,00032 |
|
|
4 |
15,40 |
12,80 |
7,90 |
-6,54 |
-14,73 |
-12,18 |
42,77 |
216,97 |
148,35 |
0,02827 |
-0,00823 |
-0,04639 |
-0,00023 |
-0,00131 |
0,00038 |
|
|
5 |
12,30 |
12,40 |
8,40 |
-9,64 |
-15,13 |
-11,68 |
92,93 |
228,92 |
136,42 |
0,00521 |
-0,01013 |
-0,03780 |
-0,00005 |
-0,00020 |
0,00038 |
|
|
6 |
6,80 |
13,10 |
10,10 |
-15,14 |
-14,43 |
-9,98 |
229,22 |
208,22 |
99,60 |
-0,03571 |
-0,00680 |
-0,00859 |
0,00024 |
0,00031 |
0,00006 |
|
|
7 |
7,90 |
15,40 |
9,70 |
-14,04 |
-12,13 |
-10,38 |
197,12 |
147,14 |
107,74 |
-0,02753 |
0,00414 |
-0,01546 |
-0,00011 |
0,00043 |
-0,00006 |
|
|
8 |
10,40 |
13,90 |
10,60 |
-11,54 |
-13,63 |
-9,48 |
133,17 |
185,78 |
89,87 |
-0,00893 |
-0,00300 |
0,00000 |
0,00003 |
0,00000 |
0,00000 |
|
|
9 |
11,60 |
14,50 |
11,40 |
-10,34 |
-13,03 |
-8,68 |
106,92 |
169,78 |
75,34 |
0,00000 |
-0,00014 |
0,01375 |
0,00000 |
0,00000 |
0,00000 |
|
|
10 |
9,80 |
14,70 |
10,10 |
-12,14 |
-12,83 |
-9,98 |
147,38 |
164,61 |
99,60 |
-0,01339 |
0,00081 |
-0,00859 |
-0,00001 |
0,00012 |
-0,00001 |
|
|
11 |
11,40 |
15,10 |
11,70 |
-10,54 |
-12,43 |
-8,38 |
111,09 |
154,50 |
70,22 |
-0,00149 |
0,00271 |
0,01890 |
0,00000 |
-0,00003 |
0,00005 |
|
|
12 |
10,60 |
11,40 |
9,90 |
-11,34 |
-16,13 |
-10,18 |
128,60 |
260,18 |
103,63 |
-0,00744 |
-0,01489 |
-0,01203 |
0,00011 |
0,00009 |
0,00018 |
|
|
13 |
11,80 |
15,90 |
10,80 |
-10,14 |
-11,63 |
-9,28 |
102,82 |
135,26 |
86,12 |
0,00149 |
0,00652 |
0,00344 |
0,00001 |
0,00001 |
0,00002 |
|
|
14 |
12,70 |
16,20 |
11,50 |
-9,24 |
-11,33 |
-8,58 |
85,38 |
128,37 |
73,62 |
0,00818 |
0,00794 |
0,01546 |
0,00007 |
0,00013 |
0,00012 |
|
|
15 |
13,70 |
16,80 |
11,50 |
-8,24 |
-10,73 |
-8,58 |
67,90 |
115,13 |
73,62 |
0,01563 |
0,01080 |
0,01546 |
0,00017 |
0,00024 |
0,00017 |
|
|
16 |
14,30 |
17,50 |
12,40 |
-7,64 |
-10,03 |
-7,68 |
58,37 |
100,60 |
58,98 |
0,02009 |
0,01413 |
0,03093 |
0,00028 |
0,00062 |
0,00044 |
|
|
17 |
14,90 |
18,90 |
12,90 |
-7,04 |
-8,63 |
-7,18 |
49,56 |
74,48 |
51,55 |
0,02455 |
0,02079 |
0,03952 |
0,00051 |
0,00097 |
0,00082 |
|
|
18 |
16,50 |
18,40 |
13,70 |
-5,44 |
-9,13 |
-6,38 |
29,59 |
83,36 |
40,70 |
0,03646 |
0,01841 |
0,05326 |
0,00067 |
0,00194 |
0,00098 |
|
|
Всього |
208,40 |
261,50 |
190,80 |
0,00 |
0,00 |
0,00 |
2059,67 |
3144,92 |
1657,75 |
-0,00298 |
-0,00019 |
0,00000 |
0,00238 |
0,00415 |
0,00442 |
|
|
В сер. |
21,94 |
27,53 |
20,08 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Крок 1.
Обчислимо елементи стандартизованих векторів . Для допоміжних розрахунків вводимо у таблиці колонки (5-10). Для нашого прикладу п=18; m=3
Необхідно згадати формули сер. арифм. та дисп
;
¯x1= Ƹx1/n = 208,40/18 = 11,6 люд.днів;
¯х2 = Ƹx2/n = 261,50/18 = 14,53 млн.грн;
¯ х3 = Ƹx3/n = 190,80/18 = 10,6%
Q²x1 = 2059,67/18 = 114,43;
Q²x2 = 3144,92/18 = 174,72;
Q²x3 = 1657,75/18 = 92,1;
Тоді знаменник для стандартизації кожної незалежної змінної буде такий:
Для : √Q²x1*n =√114,43*18 = 134,4
Для : √Q²x2*n =√174,72*18 = 210,24
Для : √Q²x3*n =√92,1*18 = 58,2
Усі розраховані елементи записуємо до таблиці та отримуємо матрицю стандартизованих змінних (колонки 11,12,13), і обчислюємо елементи кореляційної матриці r (колонки 14,15,16).
Крок 2.
Знаходимо кореляційну матрицю, що має розмір :
r= |1 0,00238 0,00415|
|0,00238 1 0,00442|
|0,00415 0,00442 1|
Кожен елемент даної матриці є коефіцієнтом кореляції, що характеризує тісноту зв’язку однієї незалежної змінної з іншою. Оскільки діагональні елементи характеризують тісноту зв’язку кожної незалежної з цією самою змінною, то вони дорівнюють 1-ці. Зауважимо, що при знаходженні добутку матриць і за рахунок зміщеності коефіцієнтів парної кореляції числові значення діагональних елементів можуть наближатись до одиниці. Якщо це так, то вони заміняються одиницями, а інші значення матриці r збільшуються на величину, що визначається як різниця між одиницею і діагональним елементом.
Інші елементи матриці r дорівнюють:
r x1x2 = 0,00238; r x1x3 = 0,00415; r x2x3 = 0,00442, тобто вони є парними коефіцієнтами кореляції між пояснювальними змінними. Користуючись цими коефіцієнтами, можна зробити висновок, що між змінними існує зв'язок. Але чи можна стверджувати, що цей зв'язок є виявленням мультиколінеарності, а тому негативно впливатиме на оцінку економетричної моделі?
Відповідь отримаємо знайшовши статистичні критерії оцінки мультиколінеарності за алгоритмом Фаррара-Глобера.
Крок 3.
Обчислимо детермінант кореляційної матриці r і критерій :
r= (-1)¹+¹*1*|1 0,00442|+(-1)¹+²*0,00238*|0,00238 0,00442|+(-1)¹+³*0,00415*|0,00238 1|
|0,00442 1| |0,00415 1 | |0,00415 0,00442|
= 0,999
x² = -[n-1-1/6(2m+5)]In|r| = -[18-1-1/6(2*3+5)]In0,999 = -91/6*In0,999 = 0,0151
При ступені свободи і рівні значущості , критерій . Оскільки , робимо висновок, що в масиві незалежних змінних мультиколінеарності не існує.
Крок 4.
Знайдемо матрицю, обернену до матриці r: , де І так звана «приєднана» матриця до матриці r. Вона складається з алгебраїчних доповнень до елементів матриці r, а саме:
Алгебраїчні доповнення і-го та j-го елемента – це мінор , який береться зі знаком . Тобто, , де (мінор) - визначник матриці (n-1)-порядку, в якій викреслено і-й рядок та j-стовпчик.
r= |1 0,00238 0,00415|
|0,00238 1 0,00442|
|0,00415 0,00442 1|
r11 = (-1)¹+¹*|1 0,00442|= 1 -0,999 = 0,1
|0,00442 1|
r12 = (-1)¹+²*|0,00238 0,00442| = 0,00236
|0,00415 1 |
r13 =(-1)¹+³*|0,00238 1 | = -0,00414
|0,00415 0,00442|
r12 = 0,00236; r21 = 0,00236 r13 = -0,00414; r31 = -0,00414 r22 = 0,999; r23 = -0,00414; r32 = -0,00414 r33 = 0,989.
Отже,
r¯¹= (X*×X*)¯¹= 1/|r|*I = 1/0,999*|0,1 0,00236 -0,00414| |0,1 0,00236 -0,00414|
|0,00236 0,999 -0,00414| = |0,00236 1 -0,00414|=C
|-0,00414 -0,00414 0,989| |-0,00414 -0,00414 0,989|
Крок 5.
Використовуючи діагональні елементи матриці С, обчислимо F-критерії:
F1 = (c11-1)((n-m)/(m-1)) = (0,1 - 1)15/2 = -6,75
F2 = (c22-1)((n-m)/(m-1) = (1-1)15/2 = 0
F3 = (c33-1)((n-m)/(m-1) = (0,989-1)15/2 = -0,0825
Для рівня значущості і (n-m)=15 ступенів вільності і (m-1)=2 критичне (табличне) значення F-критерію F=19,35
Оскільки робимо висновок, що не мультиколінеарна з та .
Оскільки робимо висновок, що не мультиколінеарна з та .
Оскільки робимо висновок, що не мультиколінеарна з та .
Для визначення наявності попарної мультиколінеарності розглянемо t-критерій Стюдента.
Крок 6.
Обчислимо частинні коефіцієнти кореляції, використавши елементи матриці С:
r12.3 = -C12/√C11*C22 = -0,00236/√0,1*1= -0,00746
r13.2 = -C13/√C11*C33 = 0,00414/√0,1*0,989= 0,0132
r23.1= -C23/√C22*C33 = 0,00414/√1*0,989 = 0,00416, тоді порівняємо
r x1x2 = 0,00238>-0,00746 = r12.3 r x1x3 = |0,00415|<|0,0132| = r13.2 r x2x3 = |0,00442|>|0,00416| = r23.1
Частинні коефіцієнти кореляції характеризують тісноту зв’язку між двома змінними за умови, що третя не впливає на цей зв'язок.
Крок 7.
Визначимо t-критерій на основі частинних коефіцієнтів кореляції.
t12 = (|r12.3|√n-m)/(√1-r²12.3) = (-0,00746-√15)/(√1-(-0,00746)²) = -3,88
t13 = (|r13.2|√n-m)/(√1-r²13.2) = (0,0132-√15)/(√1-0,0132²) = -3,86
t13 = (|r23.1|√n-m)/(√1-r²23.1) = (0,00416-√15)/(√1-0,00416²) = -3,86
Табличне значення t-критерію при рівні значущості і (n-m)=15 ступенів вільності є t0,05(15)=1,75. Оскільки , то пара змінних та не мультиколінеарна. Оскільки , то пара змінних та не мультиколінеарна.
Оскільки , то пара змінних та не мультиколінеарна.
Короткі висновки:
Якщо порушується одна з чотирьох умов, необхідних для застосування 1МНК, а саме коли між пояснювальними змінними існує лінійна залежність, то це явище називається мультиколінеарністю.
Мультиколінеарність негативно впливає на кількісні характеристики економетричної моделі або взагалі робить неможливою її побудову.
Найповніше дослідити мультиколінеарність можна з допомогою алгоритму Фаррара—Глобера.
Цей алгоритм має три види статистичних критеріїв, за якими перевіряється мультиколінеарність усього масиву незалежних змінних (2 – «хі»-квадрат); кожної незалежної змінної з усіма іншими (F-критерій); кожної пари незалежних змінних (t-критерій).