Решение.
Построим таблицу значений факторных переменных Х2 = ln A, Х3 = ln k и результативной переменной Y= ln q.
T | X2 | Х3 | У | Х | Х | У | Х Х |
0,1798 | 1,1939 | -0,1278 | 0,032 328 | 1,425 397 | 0,016 333 | 0,214 663 | |
0,2038 | 1,2030 | -0,1009 | 0,041 534 | 1,447 209 | 0,010 181 | 0,245 171 | |
0,1807 | 1,1878 | -0,1290 | 0,032 652 | 1,410 869 | 0,016 641 | 0,214 635 | |
0,1914 | 1,1314 | -0,1404 | 0,036 634 | 1,280 066 | 0,019 712 | 0,216 550 | |
0,2608 | 1,0986 | -0,0823 | 0,068 017 | 1,206 922 | 0,006 773 | 0,286 515 | |
0,2994 | 1,0543 | -0,0587 | 0,089 640 | 1,111 548 | 0,003 446 | 0,315 657 | |
0,3570 | 1,0006 | -0,0182 | 0,127 449 | 1,001 200 | 0,000 331 | 0,357 214 | |
0,3471 | 0,9969 | -0,0294 | 0,120 478 | 0,993 810 | 0,000 864 | 0,346 024 | |
0,3681 | 1,0225 | 0,0000 | 0,135 498 | 1,045 506 | 0,000 000 | 0,376 382 | |
0,4148 | 0,9783 | 0,0334 | 0,172 059 | 0,957 071 | 0,001 116 | 0,405 799 | |
ВСЕГО | 2,8029 | 10,8673 | -0,6533 | 0,856 289 | 11,879 598 | 0,075 397 | 2,978 610 |
СРЕДНЯЯ | 0,28029 | 1,08673 | -0,06533 | 0,085 628 9 | 1,187 959 8 | 0,007 539 7 | 0,297 861 |
t | Х У | Х У | У= + + Х + + Х | e=У-У | e2 | (У- )2 |
-0,022 978 | -0,152 580 | -0,128 74 | 0,000 94 | 0,000 000 88 | 0,004 021 | |
-0,020 563 | -0,121 383 | -0,100 59 | 0,000 31 | 0,000 000 10 | 0,001 243 | |
-0,023 310 | -0,153 226 | -0,130 09 | 0,001 09 | 0,000 001 19 | 0,004 194 | |
-0,026 873 | -0,158 849 | -0,140 20 | 0,000 20 | 0,000 000 04 | 0,005 606 | |
-0,021 464 | -0,090 415 | -0,080 97 | 0,001 33 | 0,000 001 77 | 0,000 245 | |
-0,017 575 | -0,061 887 | -0,057 79 | 0,000 91 | 0,000 000 83 | 0,000 057 | |
-0,006 497 | -0,018 211 | -0,018 56 | 0,000 36 | 0,000 000 13 | 0,002 187 | |
-0,010 205 | -0,029 309 | -0,030 14 | 0,000 74 | 0,000 000 55 | 0,001 238 | |
0,000 000 | 0,000 000 | 0,001 09 | 0,001 09 | 0,000 001 19 | 0,004 412 | |
0,013 854 | 0,032 675 | 0,032 65 | 0,000 75 | 0,000 00056 | 0,009 600 | |
ВСЕГО | -0,135 611 | -0,753 185 | 0,000 00724 | 0,032 803 | ||
СРЕДНЯЯ | -0,013 5611 | -0,075 318 5 | 0,000 000 724 | 0,003 280 3 |
1. Найдем коэффициенты уравнения линии регрессии
Y= + Х + Х
Как решение системы уравнений
= + +
= + +
= + +
-0,065 33 = + 0,280 29 + 1,086 73
-0,013 56 = 0,280 29 + 0,085 63 + 0,297 86
-0,075 32 = 1,086 73 + 0,297 86 + 1,187 96
Решая систему методом Гаусса получим (решение следует привести в тетради):
= -0,7614; = 1,0305; = 0,3747.
Уравнение регрессии имеет вид:
У= -0,7614 + 1,0305Х2 + 0,3747Х3.
Далее найдем дисперсии переменных У, Х2 и Х3, а также попарные ковариации этих переменных:
0,085 629 – (0,280 29)2 = 0,007 066;
1,187 960 – (1,086 73)2 = 0,006 978;
0,007 540 – (-0,065 33)2 = 0,003 272;
Cov(Х ;У) = - × = -0,013 561 1 – 0,280 29 ×(- 0,065 33) = 0,004 750;
Cov(Х ;У) = - × = -0,07 531 85 – 1,086 73 × (- 0,065 33) = -0,004 322;
Cov(Х ;Х ) = - × = 0,297 861 – 0,280 29 ×1,086 73 = -0,006 739.
2) Найдем общую дисперсию, объясненную и необъясненную регрессией дисперсии по формулам:
общая дисперсия
= D(y) = 0,003 272;
необъясненная регрессией часть дисперсии
× 0,000 007 24 = 0,000 000 724
объясненная регрессией часть дисперсии
×0,032 803 = 0,003 280 3
3) Найдем коэффициент детерминации
1 - = 0,999 8;
Скорректированный коэффициент корреляции найдем по формуле:
1 - = 0,999 7.
Значение скорректированного коэффициента корреляции близко к единице, значит скорее всего имеется высокая степень линейной зависимости значений У от факторных величин Х и Х .
Однако следует проверить наличие корреляции между факторами Х и Х . Для этого найдем частные коэффициенты корреляции.
4) Сначала найдем коэффициенты корреляции для парных регрессий по формулам:
r (x2; y) = = = 0,977 5;
r (x3; y) = = = 0,819 0.
Коэффициент частной корреляции между Уи Х при исключении влияния Х найдем по формуле:
r (у;x2 ½х3 ) = = = 0,9989.
Коэффициент частной корреляции между Уи Х при исключении влияния Х найдем по формуле:
r (у;x3 ½х2 ) = = = 0,9911.
Коэффициенты частной корреляции близки к единице, значит связь между У и факторами Х и Х тесная.
5) Проверим значимость линейной модели в целом. Для этого на уровне доверия 95% проверим нулевую гипотезу
Н: a2 = a3=0
об отсутствии линейной зависимости между У и факторами Х и Х . Для этого рассмотрим F-статистику
F = = = 17 496,5.
Найдем критическое значение F-статистики:
Fс = F(k-1;n-k) = F(3-1;10-3) = F(2;7)
а) по таблице распределения Фишера при 95% уровне доверия находим
Fc = 4,74
Поскольку F>Fc, то гипотеза Н: a2 = a3=0 с вероятностью 95% отвергается и принимается альтернативная гипотеза о значимости линейной регрессии
6) Найдем доверительные интервалы для параметров регрессии.
Составим матрицу
Q = n× = 10×
Определитель этой матрицы D = 0, 004
Найдем диагональные элементы матрицы Q-1:
q11 = ×( ) = ×(0,101724 – 0,088721) =325,06;
q33 = ×( ) = ×D(Х2) = ×0,006 978 = 174,45;
q22 = ×( ) = ×D(Х3) = ×0,007 066 = 176,65;
Используя t-статистику получим доверительные интервалы для коэффициентов a1, a2, a3
на доверительном уровне 95%:
am Î ( - tc × ; + tc × ),
где m = 1,2,3.
Критическое значение tc = t(n-k; 0,95) = t(10-3; 0,95) = t(7; 0,95) найдем по таблицам распределения Стьюдента
tc = 2,365;
Несмещенную оценку S2 дисперсии ошибок 2 найдем по формуле:
= = 0,000 001 03
Отсюда
tc × = 0,0434; tc × = 0,0320; tc × = 0,0318.
В результате получим доверительные интервалы на уровне доверия 95%:
a1 Î (-0,7614 – 0,0434; -0,7614 + 0,0434), т.е. a1 Î (-0,8048; -0,7180 );
a2 Î (1,0305 – 0,0320; 1,0305 + 0,0320), т.е. a2 Î (0,9985; 1,0625 );
a3 Î (0,3747 – 0,0318; 0,3747 + 0,0318), т.е. a3 Î (0,3429; 0,4065 ).
7) Проверим значимость регрессии. Для этого на уровне значимости 95% проверим гипотезы
Н1: α1= 0; Н2: α2= 0; Н3: α3= 0.
Составим – статистики
; ;
и найдем их значения:
, ,
Сравнивая полученные значения -статистик с критическим получим, что , , , т.е. все три гипотезы Н1, Н2 и Н3 должны быть отвергнуты на уровне доверия 95% и приняты альтернативные гипотезы, т.е. с вероятностью 95% все три коэффициента α1, α2, α3 значимы.
8) Проверим наличие мультиколинеарности. Для этого составим матрицу коэффициентов корреляции и найдем ее определитель ,
где ,
Составленный определитель равен
=1-0,000036=0,999964
На уровне значимости 95% проверим гипотезу о независимости факторных переменных
H:
Для этого рассчитаем величину
Критическое значение найдем по таблицам распределения : Пирсона при уровне значимости 95%
Поскольку , то гипотеза Н принимается на уровне доверия 95%, т.е. можно считать, что мультиколлинеарности нет.