Navch._posibnuk_Ivaschyk
.pdfПідставимо отримані значення у вихідну формулу:
|
|
|
ry |
= ryx |
= |
|
|
|
|
|
|
|
|
10 141,84 − 41,1 30,8 |
|
|
|
|
|
|
|
= |
|
|
|
|||||||||||||
|
|
|
|
( |
|
|
|
|
|
|
|
|
)( |
|
|
|
|
|
|
|
|
|
|
) |
|
|
|
|
||||||||||
|
|
|
1 |
|
1 |
|
|
|
|
|
|
|
|
|
− 41,12 |
|
|
|
|
|
−30,82 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
10 181,51 |
10 113,6 |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
= |
|
|
|
|
|
|
1418,4 −1265,88 |
|
|
|
|
= |
152,52 |
|
|
= |
|||||||||||||||||||
|
|
|
|
(1815,1−1689,21)(1136 −948,64) |
|
125,89 187,36 |
||||||||||||||||||||||||||||||||
|
|
|
= |
|
152,52 |
|
|
|
= |
152,52 |
= 0,993. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
23586,75 |
153,58 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
ry |
|
= ryx = |
|
|
|
|
|
|
|
|
|
10 152,77 − 47,1 30,8 |
|
|
|
|
|
= |
|
|
|
|
||||||||||||||||
|
( |
|
|
|
|
|
|
|
|
|
|
)( |
|
|
|
|
|
|
|
) |
|
|
|
|
|
|||||||||||||
|
|
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
47,12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
10 226,19 − |
|
10 113,6 −30,82 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
= |
|
|
|
|
|
1527,4 −1450,68 |
|
|
|
= |
|
|
76,72 |
|
= |
|
||||||||||||||||||||||
|
(2261,9 − 2218,41)(1136 −948,64) |
|
43,49 187,36 |
|
||||||||||||||||||||||||||||||||||
= |
|
76,72 |
|
|
= |
|
|
76,72 |
= 0,85. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
8148,29 |
90,27 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
r12 = rx x = |
|
|
|
|
|
|
|
|
|
|
10 200,2 − 41,1 47,1 |
|
|
|
|
|
|
= |
|
|
|
|
|
|||||||||||||||
|
( |
|
|
|
|
|
|
|
|
|
|
|
|
)( |
|
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
|||||||||
|
|
|
1 |
2 |
|
|
|
|
181,51− 41,12 |
|
|
|
|
− |
47,12 |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
10 |
|
10 226,19 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
= |
|
|
|
|
|
|
|
|
2002 −1935,81 |
|
|
|
|
= |
|
66,19 |
|
=. |
||||||||||||||||||||
|
|
(1815,1−1689,21)(2261,9 − 2218,41) |
|
125,89 43,49 |
||||||||||||||||||||||||||||||||||
= |
|
|
66,19 |
|
|
= |
|
66,19 |
= 0,89. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
5474,956 |
73,99 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Отримані коефіцієнти парної кореляції свідчать про наявність тісного зв’язку між вибраними факторами. Оскільки r12 має досить
високе значення, то є підозра про наявність явища мультиколінеарності.
16.4.3. Частинна кореляція
Коефіцієнти частинної кореляції так само представляють лінійні зв’язки ознак, але при цьому береться до уваги чистий зв’язок пари ознак за умови, що зв’язки всіх інших ознак з ознаками із зазначеної пари не діють. Отже, частинною кореляцією між ознаками xk i xj називається кореляційна залежність між цими ознаками при фіксованих значеннях інших ознак.
491
Для знаходження коефіцієнтів частинної кореляції використаємо альтернативні формули. Нехай число незалежних змінних m = 2 . Тоді мають місце такі формули:
ry1.2 = |
ry1 − ry2r12 |
, ry2.1 |
= |
ry2 − ry1r12 |
, |
(16.52) |
(1− ry22 )(1− r122 ) |
(1− ry21 )(1− r122 ) |
де ry1.2 – коефіцієнт частинної кореляції між змінними y та x1 при виключенні впливу x2; ry2.1 – коефіцієнт частинної кореляції між змінними y та x2 при виключенні впливу x1; ry1 , ry 2 , r12 – відповідно
парні коефіцієнти кореляції.
Аналогічно для випадку m = 3 маємо:
ry1.23 = |
|
ry1.3 − ry2.3r12.3 |
; ry2.13 |
= |
ry2.3 − ry1.3r21.3 |
; |
|
|
(1− ry22.3 )(1− r122.3 ) |
(1− ry21.3 )(1− r212 .3 ) |
|||||
|
|
ry3.2 − ry1.2r31.2 |
|
|
|
|
(16.53) |
ry3.12 = |
|
, |
|
|
|
||
(1− ry21.2 )(1− r312 .2 ) |
|
|
|
де ry1.23 , ry 2.13 , ry3.12 – коефіцієнти частинної кореляції третього порядку.
Узагальнюючи наведені вище формули для будь-якого числа пояснювальних змінних, отримаємо:
r |
= |
ry1.3…m − ry2.3…mr12.3…m |
. |
(16.54) |
|
(1− ry22.3…m )(1− r122.3…m ) |
|||||
y1.2…m |
|
|
|
Рекурентне співвідношення (16.54) показує, що обчислення коефіцієнта частинної кореляції порядку m зводиться до визначення таких же коефіцієнтів, але порядку (m −1). Тому при використанні
(16.54) спочатку необхідно знайти парні коефіцієнти кореляції, а потім перейти до обчислення коефіцієнтів вищих порядків.
Великим інтересом є знаходження коефіцієнтів частинної кореляції здопомогоюматричноїалгебри. Представимоматрицю(16.50) так:
r00 |
r01 |
r02 |
r0m |
|
|
|
r |
r |
r |
r |
|
|
|
10 |
11 |
12 |
1m |
|
|
|
[Rn ]= r20 |
r21 |
r22 |
r2m |
|
, |
(16.55) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
r |
r |
r |
|
|
|
m0 |
m1 |
m2 |
mm |
|
|
де індекс рівний нулю служить для відображення парних коефіцієнтів кореляції залежної змінної та пояснювальних змінних.
492
Тоді буде мати місце загальна формула для коефіцієнта частинної кореляції:
|
|
|
|
r0 j .12…j−1 j+1…m = − |
R0 j |
, |
(16.56) |
||||
|
|
|
|
R00 Rjj |
|||||||
|
|
|
|
|
|
|
|
|
|
||
де Rij означають алгебраїчні |
доповнення |
до елемента |
|||||||||
rkj (k = |
|
|
|
)матриці [Rn |
]. |
|
|
|
|
|
|
0, m, j = |
0, m |
|
|
|
|
|
|||||
Розглянемо випадок, коли результативний показник формується |
|||||||||||
під впливом двох факторів x1 та x2. Отже, маємо: |
|
||||||||||
|
|
|
|
r00 |
r01 |
r02 |
|
|
R01 |
|
|
|
|
|
[Rn ]= r10 |
r11 |
r12 |
; r01.2 = − |
; |
||||
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
R00 R11 |
|
|
|
|
|
|
r21 |
|
|
|
|
|
|
|
|
|
|
r20 |
r22 |
|
|
|
|
r02.1 |
= − |
R02 |
; r12.0 = − |
R12 |
. |
|
R00 R22 |
R11R22 |
|||||
|
|
|
|
Приклад 16.4. Розрахувати частинні коефіцієнти кореляції, взявши дані із прикладу 16.1.
♦Розв’язування.
Коефіцієнти частинної кореляції обчислюємо за формулами
(16.52):
|
|
|
|
|
|
ry1.2 = |
|
|
|
ry1 − ry2r12 |
|
|
|
, ry2.1 = |
|
ry2 − ry1r12 |
|
. |
|
||||||||||||||||||||
|
|
|
|
|
|
|
(1− ry22 )(1− r122 ) |
(1− ry21 )(1− r122 ) |
|
||||||||||||||||||||||||||||||
Значення коефіцієнтів парної кореляції ry1 , ry 2 , r12 візьмемо з |
|||||||||||||||||||||||||||||||||||||||
прикладу 16.3. |
|
|
0,993 −0,85 0,89 |
|
|
|
|
|
|
|
|
0,993 −0,757 |
|
|
|
||||||||||||||||||||||||
|
|
r |
|
|
= |
|
|
|
|
|
= |
|
|
|
|
|
= |
||||||||||||||||||||||
|
|
|
|
|
( |
|
|
|
|
|
|
|
|
)( |
|
|
|
|
|
|
|
|
|
) |
|
|
|
|
|
|
|
||||||||
|
|
y1.2 |
|
|
|
|
− |
|
|
|
|
2 |
|
− |
|
|
|
|
|
2 |
|
|
|
(1−0,7225)(1−0,7921) |
|
||||||||||||||
|
|
|
|
|
|
1 0,85 |
|
|
|
1 |
|
|
0,89 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
= |
|
|
|
|
|
0,236 |
|
|
|
|
|
|
= |
|
|
0,236 |
|
|
= |
0,236 |
= 0,983. |
|
|
|
|||||||||||||
|
|
|
0,2475 0,2079 |
|
|
0,0577 |
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
0,24 |
|
|
|
|
||||||||||||||||||||||||||
r = |
( |
0,85 −0,993 0,89 |
|
|
) |
= |
|
|
|
|
|
|
0,85 −0,884 |
= |
|
||||||||||||||||||||||||
y2.1 |
|
|
|
|
− |
|
|
|
2 |
|
|
|
− |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
(1−0,986)(1−0,7921) |
|
|
|
|||||||||
|
|
1 |
|
|
0,993 |
|
1 |
|
|
0,89 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
= |
|
|
−0,034 |
|
|
|
|
= |
|
|
−0,034 |
|
|
= |
−0,034 = −0,63. |
|
|
|
|||||||||||||||||||||
|
0,014 0,2079 |
|
|
|
0,0029 |
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
0,054 |
|
|
|
|
|
493
16.5. Оцінка якості економетричних моделей
Якість лінійних економетричних моделей оцінюється стандартним для економіко-математичних задач способом: за адекватністю та точністю. Адекватність регресійних моделей може бути визначена на основі статистичного аналізу залишкової послідовності. При цьому їх значення отримують внаслідок підстановки у модель фактичних значень всіх включених до моделі факторів.
Для оцінки точності регресійних моделей використовуються статистичні критерії значущості Фішера та Стьюдента.
Стандартні помилки оцінок. Якість вибраної функції регресії можна оцінити на основі стандартних помилок або дисперсій залишків оцінок параметрів моделі.
Розглянемо спочатку алгоритм оцінки якості моделі з допомогою залишків послідовності. Стандартна помилка залишків називається також стандартною помилкою оцінки регресії у зв’язку з інтерпретацією величини u (збурення) як результату помилки специфікації функцій регресії.
Збурююча змінна u є випадковою з визначеним розподілом ймовірності. Математичне сподівання цієї змінної рівне нулю, а дисперсія – σu2 на основі передумов застосування МНК. Таким чином,
σu2 – це дисперсія збурення у генеральній сукупності. Проте нам не відомі значення збурення. Про нього можна судити лише на основі залишків e. Знайдена за цими залишками дисперсія σe2 буде оцінкою
дисперсії збурюючої змінної. Тоді незміщену оцінку дисперсії збурення знайдено з формули:
|
|
n |
|
|
|
|
|
|
|
|
|
2 |
|
∑ei2 |
|
|
1 |
|
|
′ |
|
|
|
|
i=1 |
|
|
|
|
|
|
||||
σe |
= |
|
|
= |
|
|
|
e |
e. |
(16.57) |
|
n −(m +1) |
n − m |
−1 |
|||||||||
|
|
|
|
|
|
||||||
У знаменнику |
формули |
міститься |
число ступенів |
вільності |
n −(m +1), де n – обсяг вибірки, m – число пояснювальних змінних.
Такий вираз числа ступенів вільності зумовлений тим, що число
залишків повинно задовольняти m+1 |
умовам. |
|
|||
Для безпосереднього проведення розрахунків можна |
|||||
використати такі формули: |
|
|
|
||
n |
n |
n |
n |
n |
|
∑ei2 |
= ∑yi2 − a0 |
∑yi − a1 |
∑xi1 yi −…− am ∑xim yi |
(16.58) |
|
i=1 |
i=1 |
i=1 |
i=1 |
i=1 |
|
494
або в матричній формі:
n
∑ei2 = e′e =Y′Y − A′X ′Y .
i=1
Вираження сум у правій частині (16.58) містяться в робочій таблиці МНК, а оцінки параметрів вже знайдено. Беручи до уваги поняття коефіцієнта детермінації, стає зрозумілим фізичний зміст дисперсії (стандартного відхилення) залишків – це та частка загальної дисперсії σ2y , яка не може бути пояснена залежністю змінної від
факторів x j (j =1,m).
Приклад 16.5. Для економетричної моделі (прикл. 16.1) обчислити незміщену оцінку дисперсії залишків.
♦ Розв’язування.
Незміщену оцінку дисперсії залишків розрахуємо за формулою
(16.38).
|
n |
|
|
|
|
|
|
|
∑ei2 = e′ e =Y′ Y − A′ X ′ Y . |
|
|
|
|||
|
i=1 |
|
|
|
|
|
|
Обчислимо добуток Y ′ Y : |
|
|
|
|
|
|
|
|
|
|
|
|
1,2 |
|
|
|
|
|
|
|
1,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1,9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2,2 |
|
|
|
|
|
|
|
2,8 |
|
=113,6 . |
Y′ Y = (1,2 1,5 1,9 2,2 2,8 3,1 3,4 4,3 4,8 5,4)× |
3,1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3,4 |
|
|
|
|
|
|
|
4,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4,8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5,4 |
|
|
|
− 0,36). |
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
A = (− 0,97 1,4 |
|
|
|
|
|
|
|
Значення добутку (X ′ Y ) візьмемо з прикладу 16.1: |
|
|
|||||
|
|
|
30,8 |
|
|
|
|
|
X ′ Y = |
|
141,84 |
|
|
|
|
|
|
. |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
152,77 |
|
|
|
|
495
Тоді
|
|
|
|
|
30,8 |
|
|
′ |
X |
′ |
Y = (− 0,97 1,4 |
|
141,84 |
|
=112,175. |
A |
|
− 0,37)× |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
152,77 |
|
|
Отже, незміщена оцінка дисперсії залишків буде:
10
∑ei2 =113,6 −112,175 =11425, ≈1,4 .♦
i=1
Тепер перейдемо до розгляду алгоритму оцінки якості моделі з допомогою стандартних помилок або дисперсій оцінок параметрів моделі. При цьому вважаємо, що вже знайдено значення оцінок параметрів моделі. Дані оцінки є випадковими величинами, що мають певний розподіл імовірностей. Можливі значення оцінок розсіюються навколо дійсного значення параметра В. Знайдемо міру розсіювання оцінок параметрів моделі, використавши для цього математичний апарат матричної алгебри. Окреслений апарат допомагає знайти не тільки дисперсії оціночних параметрів А, але й розрахувати значення коваріацій між двома параметрами аk і aj k ≠ j,k = 0,m, j = 0,m . Ці величини служать характеристиками випадкових змінних aj і утворюють дисперсійно-коваріаційну матрицю оцінок параметрів моделі.
За означенням коваріаційна матриця для оператора А буде:
cov(А)= M (A − B) (A − B)′ = σu2 |
(X ′ X )−1. (16.59) |
|
|
|
|
Для доведення цього твердження виконаємо такі перетворення:
A − B = (X ′ X )−1 X ′Y − B = (X ′ X )−1 X ′(XB +u)− B =
= (X ′ X )−1 X ′XB +(X ′ X )−1 X ′ u − B =
(16.60)
= EB +(X ′ X )−1 X ′ u − B = (X ′ X )−1 X ′ u.
Підставимо отримане значення (16.60) у вираз (16.59). Одержимо:
496
|
′ |
−1 |
′ |
u |
|
′ |
−1 |
′ |
u |
′ |
= |
|
cov(A)=M (X X ) |
|
X |
(X X ) |
|
X |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
=M (X′X )−1 X′ u u′X (X′ X )−1 =M (X′ X )−1 X′ M (uu′)X (X′ X )−1 =
=(X′ X )−1 X′Xσu2 (X′ X )−1 =σu2 (X′ X )−1 X′ X (X′X )−1 =σu2 (X′X )−1 .
Позначимо через σа дисперсійно-коваріаційну матрицю:
|
′ |
|
|
σa2 |
… σa a |
|
2 |
−1 |
|
1 |
1 |
m |
|
σa = σu |
(X X ) |
|
= |
|
|
. |
|
|
|
σa a |
… σa2 |
|
|
|
|
|
|
m 1 |
m |
(16.61)
(16.62)
Ця матриця є симетричною і на головній діагоналі містяться
дисперсії оцінок параметрів моделі aj (j =1,m |
): |
|
|
|
|
|||
|
σa2 j |
= M (aj −bj )2 , |
(16.63) |
|||||
а поза діагоналлю – їхня коваріація: |
|
|||||||
σa j ak |
= M [(aj −bj )(ak |
−bk )], j ≠ k; j = |
|
|
|
|
(16.64) |
|
0,m;k = 0,m. |
||||||||
Оскільки |
σu2 невідоме, |
використаємо в подальшому його |
незміщену оцінку σe2 і отримаємо оцінку для матриці σa :
σˆ 2
a1
SA =σe2 (X ′ X )−1 = σˆama1
... |
σˆa a |
|
|
|
|
1 m |
|
|
|
|
|
|
, |
(16.65) |
... |
σˆ 2 |
|
|
|
|
am |
|
|
|
де σˆaj a k – |
оцінка коваріацій між aj та |
|
ak, σˆa2j – |
оцінка дисперсій |
|||||||||||||||||
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
параметрів моделі a |
, σ2 = |
|
|
|
e e |
. |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
j |
e |
n − m −1 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Враховуючи (16.38), вираз для σe2 |
матиме вид: |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
Y |
′ |
|
|
′ |
′ |
|
|
|
|
|
|
|
|
|
|
|
|
σe2 |
|
= |
|
Y − A X Y |
. |
|
|
|
|
|
|
|
|
(16.66) |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
n − m −1 |
|
|
|
|
|
|
|
|
|
||||
Для проведення кінцевих розрахунків позначимо (k, j) елемент |
|||||||||||||||||||||
матриці (X ′ X )−1 через |
Ckj. Внаслідок |
чого отримаємо |
матрицю |
||||||||||||||||||
C = (X ′ X )−1 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тоді |
дисперсії оцінок |
|
параметрів |
заданої |
функції |
регресії |
|||||||||||||||
знайдемо за формулами: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
var (aj )= σˆ a2 j = σe2Cjj , j = |
|
|
, |
|
|
|
|
(16.67) |
|||||||||||
|
|
0,m |
|
|
|
||||||||||||||||
|
cov(ak aj )= σˆ a2 |
|
|
= σe2Ckj , j ≠ k, j = |
|
|
|
|
(16.68) |
||||||||||||
|
a |
|
0,m,k = 0,m . |
||||||||||||||||||
|
|
|
|
k |
j |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
497
Стандартна помилка оцінки параметра регресії використовується для оцінки якості вибору функції регресії. Позначимо її через Sa j , тоді
|
Sa j = |
σˆ a2 j , j = |
0,m |
. |
(16.69) |
|||||||
На завершення аналізу знайдемо відносний показник |
||||||||||||
розсіювання |
|
|
|
Sa j |
|
|
|
|
|
|
|
|
δa j |
= |
|
|
|
|
|
100%, j = |
|
(16.70) |
|||
|
|
|
|
|
0,m, |
|||||||
|
|
aj |
|
|
|
|||||||
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
який характеризує, наскільки вдало обрахована оцінка відповідного параметра моделі. Якщо він малий, то це вказує на її незміщеність, у протилежному випадку – вона зміщена. Невиконання умови M(e)=0 якраз і породжує зміщеність. Така ситуація виникає внаслідок того, що залишки мають систематичну складову, яка викликана некоректною специфікацією процесу дослідження. Отже, чим вища відносна стандартна помилка оцінки параметра, тим більше оцінена величина відрізняється від сукупності значень залежної змінної, і тим менш надійні оцінки прогнозу, розрахованої з допомогою окресленої моделі.
Опишемо основні чинники, від яких залежить стандартна помилка коефіцієнта регресії. Їх можна розділити на різні складові:
1)розсіювання залишків. Чим більша частка варіації значень змінної Y, непояснена її залежністю від x, тим більша стандартна помилка коефіцієнта регресії. Отже, чим більше фактичні значення змінної Y відхиляються від розрахункових значень регресії, тим менш точною є знайдена оцінка параметра;
2)розсіювання значень пояснювальної змінної x. Чим сильніше зазначене розсіювання, тим менша стандартна помилка коефіцієнта
регресії. Звідси випливає, що при витягнутій хмарці точок на діаграмі розсіювання має надійнішу оцінку функції регресій, ніж при незначному скупченні точок, близько розміщених одна від одної.
3) об’єм вибірки. Чим більший об’єм вибірки, тим менша стандартна помилка коефіцієнта регресії. Тут існує безпосередній зв’язок оцінки параметра моделі з властивістю її асимптотичності незміщеності.
Приклад 16.6. Провести оцінку якості економетричної моделі (приклад 16.1) з допомогою процедури стандартних помилок оцінок параметрів моделі.
498
♦ Рoзв’язування.
Стандартні помилки оцінок параметрів моделі. Обчислимо за формулою (16.69)
|
|
|
|
|
|
|
|
Sa j |
= σˆ a2 j |
, j = |
0,m |
,σˆ a2 j = σe2 Cij , |
||||||
|
Y |
′ |
|
|
|
′ |
′ |
|
|
|
|
|
|
|
|
|
|
|
де σe2 = |
|
Y − A X Y |
, а Сij – |
|
елементи головної діагоналі матриці |
|||||||||||||
|
|
n − m −1 |
|
|||||||||||||||
C=(X’ X)-1. |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Значення Cjj |
візьмемо з прикладу 16.1: |
|||||||||||||||||
|
|
|
|
|
|
|
|
C00 = 8,92; C11 = 0,4; C22 =1,51. |
||||||||||
Значення виразу Y |
′ |
|
|
|
|
′ |
′ |
|
|
|||||||||
|
Y − A X Y ми обчислили в прикл. 16.5: |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
Y |
′ |
|
|
′ |
′ |
||
|
|
|
|
|
|
|
|
|
|
|
|
Y − A X Y =1,4. |
||||||
Тоді σ2 |
= |
|
|
1,4 |
|
|
= |
1,4 |
= 0,2, а дисперсії оцінок a : |
|||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
e |
|
10 |
− 2 −1 |
7 |
|
|
|
|
|
|
j |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
σˆa2 |
= 0, 2 8,92 =1,78; |
Sa = 1,78 =1,33 , |
|||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
σˆ a2 = 0,2 0,4 = 0,08; Sa = 0,08 = 0,28, |
|||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
σˆ a2 = 0,2 1,51 = 0,3; Sa = 0,3 = 0,55 . |
|||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
Далі знайдемо відносний показник розсіювання δaj , тобто порівняємо кожну стандартну помилку Saj з відповідним числовим значенням оцінки параметра.
δa0 |
= |
|
Sa0 |
|
|
|
100% = |
1,33 |
|
|
|
|
|
100% =137,11% ; |
|||||
|
a0 |
|
|
|
|
−0,97 |
|
|
|||||||||||
|
|
|
|
||||||||||||||||
δa |
= |
Sa1 |
|
|
100% = |
0,28 |
100% = 20%; |
||||||||||||
|
a1 |
|
|
|
|
1,4 |
|
|
|||||||||||
|
|
|
|
|
|
|
|||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Sa |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δa2 |
= |
|
2 |
|
|
|
100% = |
0,55 |
|
|
100% =14,86% . |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
a2 |
|
|
|
|
−3,7 |
|
|
|||||||||||
|
|
|
|
|
Отже, стандартні помилки оцінок параметрів відносно рівня оцінок параметрів становлять відповідно 137,11 %, 20,0 % і 14,86 %, а це є підтвердженням зміщеності оцінок.
Значущість економетричної моделі. Для перевірки адекватності множинної регресійної моделі, як і у випадку парної регресії, використовується F-критерій Фішера. У цьому випадку нульова гіпотеза узагальнюється:
H0: a1 = a2 = ... = am = 0.
499
Тоді альтернативною гіпотезою буде Н1 : хоча б одне значення aj відмінне від нуля. У випадку невиконання гіпотези Н0 приймається гіпотеза Н1. Отже, не всі параметри незначною мірою відрізняються від нуля. Це свідчить про те, що включені до моделі фактори пояснюють змінну результативного показника.
Для перевірки гіпотези Н0 використовують F-критерій Фішера, з (m–1) та (n–m–1) ступенями вільності:
|
|
n |
|
y ) |
|
|
n |
|
y ) |
|
|
n |
|
y ) |
|
|
|
∑(yi |
− |
2 |
|
∑(yi |
− |
2 |
|
∑(yi |
− |
2 |
|||
|
|
ˆ |
|
|
|
|
|
ˆ |
|
||||||
F = F |
= |
i=1 |
|
|
|
: |
i=1 |
|
|
|
= |
i=1 |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
||||
m−1,n−m−1 |
|
m |
−1 |
|
|
n − m −1 |
|
|
|
|
|
||||
|
|
|
|
|
|
∑(yi − y )2 |
i=1
(n − m −1), (16.71)
(m −1)
де m – кількість незалежних факторів, які включено до моделі разом із фіктивним, n – загальна кількість спостережень.
Можна показати що має місце альтернативне представлення окресленого показника у матричній формі:
|
|
|
|
|
′ |
′ |
|
|
|
|
F |
|
= |
|
A X Y (n − m −1) |
|
. |
(16.72) |
|
|
|
(Y Y − A X Y )(m −1) |
|||||||
|
m−1,n−m−1 |
|
|
′ |
′ ′ |
|
|
||
Далі для заданого рівня значущості α і ступенів вільності k1=m– |
|||||||||
–1 і k2=n–m–1 |
знаходимо табличне значення критерію Фішера – |
||||||||
Fтабл.(k1,k2,α). |
|
Знайдене |
|
розрахункове |
значення |
критерію |
|||
Fm−1,n−m−1 = Fрозр. . |
Порівнюємо з табличним: |
якщо Fрозр > Fтабл., тоді |
гіпотеза Н0 відхиляється і приймається альтернативна, що свідчить про адекватність побудованої моделі, іншими словами, підтверджується наявність істотного зв’язку між залежною та незалежними змінними побудованої економетричної моделі. У протилежному випадку вона приймається і модель вважається неадекватною. ♦
Приклад 16.7. Перевірити гіпотезу про значущість економетричної моделі (прикл.16.1)
♦ Розв’язування.
Висуваємо нульову гіпотезу H0: a1 = a2 = 0.
Для перевірки гіпотези про значущість економетричної моделі використаємо F-критерій Фішера.
Обчислимо розрахункове значення критерію за формулою
(16.72).
Fрозр = (Y ′Y - A′X ′Y )(m −1).
500