эконометрика для очно-заочного 4 курс 2024-2025 год / Эконометрика 2023
.pdf
не значимы, т.е. равны нулю в генеральной совокупности. Значит, в модели №2 могут быть исключены факторы 4 , 6 .
Вывод: модели №1 и №2 значимы в целом (по тесту Фишера). Но в модели №2 большее количество незначимых оценок параметров (по тесту Стьюдента). В модели №1, ввиду большего количества значимых оценок параметров, предположительно, допущено меньше ошибок спецификации.
Б) По имеющимся в таблице 73 данным можно рассчитать ряд характеристики качества аппроксимации (формулы (49)-(61)). Сделаем расчеты для модели №1.
|
|
|
̅ |
∑ |
|
|
|
|
0 |
|
0; |
|
|
|
|
|
∑ |
|
| |
| |
22,46 |
|
0,432; |
|
|
|
|||||
|
|
|
|
|
52 |
|
52 |
|
|
|
|
|
|
|
52 |
|
52 |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
∑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 3 |
|
18,18 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
0,38 0,62; |
|
exp |
∙ |
|
||||||||||||||||
|
|
|
|
|
|
|
52 |
3 |
|
1 |
|
52 |
52 |
|
|||||||||||||||||
|
|
|
|
52 ln |
|
0,35; |
52 |
|
|
∙ |
18,18 |
|
0,44; |
, ∙ 100% |
23,8% |
|
|||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
, |
|
|
6 ln ln 52 |
52 |
46,4; ̅ |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
̅ |
24,49 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
∑ |
|
|
|
|
|
|
|
|
|
|
|
̃ |
52 |
|
∙ 100% |
47% |
|
|
|
|
|
|
|
||||||||
|
|
расчета коэффициента детерминации необходимо найти |
|||||||||||||||||||||||||||||
Для |
|
. |
|
Для этого |
|
воспользуемся соотношением |
|
|
|
|
|||||||||||||||||||||
|
∑ |
|
|
|
|
, |
|
18,18 64,83 |
|
|
. |
|
|
|
|
|
|||||||||||||||
∑ |
|
|
|
|
|
|
83,01 |
|
|
|
|
|
. |
|
|||||||||||||||||
|
|
|
1 |
|
|
|
1 |
, |
|
|
|
0,7809; |
|
|
1 |
0,2191 |
|
0,772 |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Аналогично рассчитываются характеристики для модели №2. Все результаты сведены в таблицу 75.
Таблица 75– Показателикачествааппроксимациимоделейдлярешениязадача1. Контрольная работа 2 по темам 4,5.
№ |
|
|
|
|
|
|
|
|
|
|
|
модели |
|
|
Характеристики |
качества |
|
|
,% |
|
|
||
|
|
|
аппроксимации |
|
|
||||||
1 |
0 |
0,432 |
0,62 |
0,35 |
0,44 |
-46,4 |
23,8 |
47 |
78,09 |
77,2 |
|
|
̅ |
|
|
|
|
HQC |
,%̅ |
|
|
,% |
,% |
2 |
0 |
0,392 |
0,57 |
0,34 |
0,38 |
-53,72 |
17,6 |
39 |
80,97 |
79,78 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Ответ: коэффициенты детерминации модели №2 больше, а остальные характеристики меньше, чем у модели №1. Это означает, что модель №2 лучше
101
аппроксимирует ряд рекомендованных розничных цен, чем модель №1. Поэтому, несмотря на наличие большего числа незначимых оценок параметров, модель №2 лучше объясняет изменения цены на автомобили ( ).
Задача 2. Контрольная работа 2 по теме 5. |
|
|
|
|
|
|
|
|
|
|
|
||||||
Имеются три временных ряда: индекс потребительских цен (ИПЦ) |
|
, цена |
|||||||||||||||
на нефть марки Brent |
|
|
, номинальный обменный курс рубля к |
доллару |
|
||||||||||||
|
|
|
|
|
|
|
|
2 |
|||||||||
(таблица 76). |
Используя спецификацию модели |
|
|
|
|
|
|
|
|
||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
спрогнозируйте в каком диапазоне с вероятностью 99% |
окажется значение ИПЦ |
||||||||||||||||
|
1 |
2 |
|
||||||||||||||
в ноябре и декабре 2017 года. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Таблица 76 – Исходные данные к задаче 2. Контрольная работа 2 по теме 5 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Цена на |
Курс, |
ИПЦ, в |
|
|
|
|
Цена на |
Курс, |
|
ИПЦ, в |
||||||
Период |
нефть, |
|
руб. за |
% к пред. |
|
Период |
|
|
нефть, |
руб. за |
% к пред. |
||||||
|
$/баррель |
|
1$ |
мес. |
|
|
|
$/баррель |
1$ |
|
|
мес. |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Январь 2016 |
45,93 |
|
32,49 |
102,37 |
Январь 2017 |
|
|
|
71,18 |
29,84 |
|
|
101,64 |
||||
Февраль 2016 |
45,84 |
|
35,81 |
101,65 |
Февраль 2017 |
|
|
|
78,03 |
30,16 |
|
|
100,86 |
||||
Март 2016 |
48,68 |
|
34,66 |
101,31 |
Март 2017 |
|
|
|
82,17 |
29,56 |
|
|
100,63 |
||||
Апрель 2016 |
50,64 |
|
33,58 |
100,69 |
Апрель 2017 |
|
|
|
87,35 |
29,19 |
|
|
100,29 |
||||
Май 2016 |
65,80 |
|
31,99 |
100,57 |
Май 2017 |
|
|
|
74,60 |
30,43 |
|
|
100,50 |
||||
Июнь 2016 |
69,42 |
|
31,06 |
100,60 |
Июнь 2017 |
|
|
|
74,66 |
31,17 |
|
|
100,39 |
||||
Июль 2016 |
71,52 |
|
31,51 |
100,63 |
Июль 2017 |
|
|
|
78,26 |
30,68 |
|
|
100,36 |
||||
Август 2016 |
69,32 |
|
31,65 |
100,00 |
Август 2017 |
|
|
|
74,42 |
30,35 |
|
|
100,55 |
||||
Сентябрь 2016 |
68,92 |
|
30,86 |
99,97 |
Сентябрь 2017 |
|
|
|
82,11 |
30,81 |
|
|
100,84 |
||||
Октябрь 2016 |
75,09 |
|
29,46 |
100,00 |
Октябрь 2017 |
|
|
|
83,75 |
30,32 |
|
|
100,50 |
||||
Ноябрь 2016 |
78,36 |
|
28,90 |
100,29 |
Ноябрь 2017 |
|
|
|
84,4 |
30,99 |
|
|
|
??? |
|||
Декабрь 2016 |
77,93 |
|
29,96 |
100,41 |
Декабрь 2017 |
|
|
|
89,7 |
30,86 |
|
|
|
??? |
|||
Рассчитаем вектор оценок коэффициентов регрессии.
|
|
45,96 |
45,84 |
|
… |
83,75 |
45,96 |
32,49 |
112911,6 |
47950,74 |
||||
|
|
|
45,84 |
35,81 |
||||||||||
32,46 |
35,81 |
|
… |
30,32 |
|
… |
… |
47950,74 |
21357,02 |
|||||
|
|
|
83,75 |
30,32 |
||||||||||
|
|
|
|
|
|
|
|
0,00019 |
0,00043 |
|
||||
|
|
|
|
|
|
|
0,00043 |
0,001007 |
|
|||||
|
|
|
|
|
45,96 |
|
45,84 |
… |
83,75 |
102,37 |
156365,3 |
|||
|
|
|
|
|
|
101,65 |
||||||||
|
|
|
32,46 |
|
35,81 |
… |
30,32 |
… |
68923,24 |
|||||
|
|
|
|
|
|
100,5 |
||||||||
|
|
|
|
|
|
0,00019 |
0,00043 |
156365,3 |
0,31 |
|||||
|
|
|
0,00043 |
0,001007 |
68923,24 |
2,54 |
||||||||
102
Итоговый вид модели для прогнозирования ИПЦ:
|
|
|
|
0,31 1 |
2,54 2 |
|
|
|
||
Рассчитаем стандартную ошибку модели. Произведем дополнительные |
||||||||||
расчеты в таблице 77. |
31,44 |
12,37 |
|
6,13 |
101,56 |
|
|
|||
|
∑ |
|
|
5,078, |
||||||
|
|
|
22 |
2 |
|
|
20 |
|
||
|
|
|
|
|
5,078 |
2,25. |
|
|
|
|
Рассчитаем стандартную ошибку прогноза |
для ноября и декабря 2017 г. |
|||||||||
Вектор значений независимых переменных для ноября |
2017: |
|
84,4 . |
|||||||
|
|
|
|
|
|
|
|
|
30,99 |
|
Таблица 77 – Расчет стандартной ошибки модели для задачи 2. Контрольная работа 2 по теме 5
1 |
1 |
|
2 |
|
|
|
|
13 |
1 |
2 |
|
|
|
|
|
|
||||
t |
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
|
45,93 |
32,49 |
|
102,37 |
|
96,76 |
31,44 |
|
71,18 |
29,84 |
101,64 |
|
97,86 |
14,29 |
||||||
2 |
45,84 |
35,81 |
|
101,65 |
|
105,17 |
12,37 |
14 |
78,03 |
30,16 |
100,86 |
100,80 |
0,00 |
|||||||
3 |
48,68 |
34,66 |
|
101,31 |
|
103,13 |
3,30 |
15 |
82,17 |
29,56 |
100,63 |
100,56 |
0,01 |
|||||||
4 |
50,64 |
33,58 |
|
100,69 |
|
100,99 |
0,09 |
16 |
87,35 |
29,19 |
100,29 |
101,22 |
0,87 |
|||||||
5 |
65,80 |
31,99 |
|
100,57 |
|
101,65 |
1,17 |
17 |
74,60 |
30,43 |
100,50 |
100,42 |
0,01 |
|||||||
6 |
69,42 |
31,06 |
|
100,60 |
|
100,41 |
0,04 |
18 |
74,66 |
31,17 |
100,39 |
102,32 |
3,71 |
|||||||
7 |
71,52 |
31,51 |
|
100,63 |
|
102,21 |
2,49 |
19 |
78,26 |
30,68 |
100,36 |
102,19 |
3,34 |
|||||||
8 |
69,32 |
31,65 |
|
100,00 |
|
101,88 |
3,54 |
20 |
74,42 |
30,35 |
100,55 |
100,16 |
0,15 |
|||||||
9 |
68,92 |
30,86 |
|
99,97 |
|
99,75 |
0,05 |
21 |
82,11 |
30,81 |
100,84 |
103,71 |
8,25 |
|||||||
10 |
75,09 |
29,46 |
|
100,00 |
|
98,11 |
3,59 |
22 |
83,75 |
30,32 |
100,50 |
102,98 |
6,13 |
|||||||
11 |
78,36 |
28,90 |
|
100,29 |
|
97,70 |
6,72 |
23 |
84,40 |
30,99 |
|
|
104,88 |
|
||||||
12 |
77,93 |
29,96 |
|
100,41 |
|
100,26 |
0,02 |
24 |
89,70 |
30,86 |
|
|
106,19 |
|
||||||
|
Ошибка прогноза на ноябрь (формула 66): |
|
|
|
|
|
|
|
||||||||||||
точ |
крит |
; |
1 |
∙ |
точ |
|
точ |
|
крит |
; |
|
1 ∙ |
(66) |
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|||
Где s – стандартная ошибка модели; – вектор-столбец значений переменныхвнаблюдениидляпрогноза(первыйэлементвектораравенединице,
103
если в модели есть константа); Z – расширенная матрица независимых переменных.
|
5,078 |
1 |
84,4 |
30,99 |
0,00019 |
|
0,00043 |
84,4 |
5,29 |
0,00043 |
0,001007 |
30,99 |
|||||||
Аналогично, ошибка прогноза на декабрь |
|
5,38. |
|
|
|||||
Рассчитаем |
интервальный прогноз с 99% |
-й надежностью. Табличное |
|||||||
значение t-статистики Стъюдента для уровня значимости 1% и 20 степеней
|
инт |
|
0,01; 20 |
2,85 |
|
|
|
|
свободы |
КРИТ |
|
КРИТ |
∙ |
. |
2,85 ∙ 5,29 |
104,88 |
15,08 |
|
инт |
104,88 |
||||||
|
|
|
КРИТ |
∙ |
106,19 |
2,85 ∙ 5,38 |
106,19 |
15,33 |
Ответ: значение ИПЦ в ноябре 2107 года будет находиться в диапазоне от
89,8% до 119,96%, а в декабре 2017 года от 90,86% до 121,52%.
ТЕМА 6. ОШИБКИ СПЕЦИФИКАЦИИ ФАКТОРНОЙ РЕГРЕССИОННОЙ МОДЕЛИ
6.1 Мультиколлинеарность
Задача для решения в классе (пример 6.1).
Врезультатеанализавзаимосвязей четырех показателей получена матрица коэффициентов линейной корреляции Пирсона (таблица 78).
Таблица 78 – корреляционная матрица переменных для примера 6.1
|
X1 |
X2 |
X3 |
X4 |
|
X1 |
|
0,99 |
0,41 |
0,90 |
Коэффициент корреляции |
|
|
-37 |
-37 |
-37 |
(Число наблюдений) |
|
|
0,0000 |
0,0112 |
0,0000 |
P-значение |
X2 |
0,99 |
|
0,43 |
0,90 |
Коэффициент корреляции |
|
-37 |
|
-37 |
-37 |
(Число наблюдений) |
|
0 |
|
0,0085 |
0,0000 |
P-значение |
X3 |
0,41 |
0,43 |
|
0,41 |
Коэффициент корреляции |
|
-37 |
-37 |
|
-37 |
(Число наблюдений) |
|
0,0112 |
0,0085 |
|
0,0114 |
P-значение |
X4 |
0,90 |
0,90 |
0,41 |
|
Коэффициент корреляции |
|
-37 |
-37 |
-37 |
|
(Число наблюдений) |
|
0,0000 |
0,0000 |
0,0114 |
|
P-значение |
104
А) Пусть все перечисленные переменные являются факторами (независимыми переменными). Каковы выводы о наличии мультиколлинеарности в такой модели?
Б) Какую из четырех переменных нельзя рекомендовать в качестве зависимой переменной регрессионной модели и почему?
В) На основе корреляционной матрицы дайте спецификации конкурирующих регрессионных моделей, выбрав любую зависимую переменную (с учетом вывода по второму вопросу задачи).
Примем следующие обозначения переменных:
X1 – продажи;
X2 – издержки;
X3 – затраты на НИОКР;
X4 – затраты капитала.
Решение А) Мультиколлинеарность факторов возникает, когда в модели линейной
множественнойрегрессииприсутствуетзначимаякорреляциямеждуфакторами. Коэффициент корреляции значим, если соответствующее P-значение меньше уровня значимости (0,05). В данной корреляционной матрице все P-значения меньше 0,05, следовательно, все коэффициенты корреляции значимы. При этом присутствуют значения коэффициентов, по модулю близкие к единице (X1 и X2: 0,99, X1 и X4: 0,90), что говорит об очень сильной связи факторов. Следовательно, мультиколлинеарность есть.
Ответ на пункт А: в модели множественной регрессии присутствует мультиколлинеарность факторов.
Б) В модели множественной регрессии корреляция Y с независимыми переменными должна быть значимой, а корреляция независимых переменных (факторов) между собой должна быть незначимой, либо быть ниже, чем корреляция факторов с Y. В приведенной корреляционной матрице все коэффициенты корреляции значимы. При этом корреляция переменной X3 со всеми остальными относительно слабая (от 0,41 до 0,43). Корреляция остальных переменных между собой либо сопоставима с этими цифрами (X2 и X3: 0,43), либо гораздо сильнее (X1 и X2: 0,99, X1 и X4: 0,90).
Ответ на пункт Б: Переменную X3 (затраты на НИОКР) нельзя рекомендовать в качестве зависимой переменной множественной регрессии.
105
В) Было установлено, что X3 не может быть зависимой переменной. Пусть зависимая переменная – X1 (продажи). Тогда составим спецификации моделей, руководствуясь принципами:
1.Связь X1 с факторами должна быть значимой.
2.Связь факторов друг с другом должна быть незначимой, либо не
должна быть сильнее связи X1 с каждым из факторов.
Из факторов модели (X2, X3, X4) все связаны друг с другом. При этом между X2 и X4 связь очень сильная, коэффициент корреляции по модулю близок к 1. Включать эту пару факторов в одну модель нельзя. Составим список всех возможных моделей множественной регрессии и вычеркнем из него модели с
мультиколлинеарностью: |
|
|
|
|
|
|
факторы X2, X4 в одной модели |
|||||||
1. |
1 |
|
|
2 |
3 |
4 |
||||||||
2. |
1 |
|
|
|
|
|
|
|
|
|
|
|
по модулю |
|
|
|
2 |
3 |
корреляция между X2 и X3 |
||||||||||
|
|
|
|
|
|
|
|
|
||||||
|
больше, чем корреляция между X1 и X3 (0,43>0,41) |
|
||||||||||||
3. |
1 |
|
|
|
|
2 |
|
|
4 |
факторы X2, X4 в одной модели |
||||
4. |
|
|
|
корреляция между X3 и X4 |
по модулю |
|||||||||
|
1 |
|
|
3 |
4 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сравнима с корреляцией между X1 и X3 (0,41=0,41)
Таким образом, модели множественной регрессии для зависимой переменной X1 на этих данных строить нельзя.
В парной линейной регрессии R-квадрат равен квадрату парного линейного коэффициента корреляции. Запишем все возможные модели парной регрессии для X1 с указанием R-квадрата:
1. |
1 |
|
, R2 = 98,8% |
2. |
2, R2 = 17,0% |
||
3. |
1 |
|
3, R2 = 81,7% |
Строить1 |
парную |
4 |
|
всего две модели. |
регрессию с R2 = 17% не имеет смысла, поэтому остается |
||
|
|||
Ответ на пункт В: |
|||
2. |
1 |
|
2 |
1. |
1 |
|
4 |
|
|||
6.2 Гетероскедастичность
Задача для решения в классе (пример 6.2)
Имеется статистика по двадцати странам по соответствию уровней государственных затрат на здравоохранение объему производимого за год валового выпуска (в млрд ден. ед.).
106
В Excel была построена линейная регрессионная модель, описывающая зависимость государственных расходов от валового внутреннего продукта:
0,05 ∙ .
Исходные данные и представлены в таблице 79 , отчет Excel представлен на рисунке 32. Необходимо проверить остаточную компоненту е на предмет присутствия гетероскедастичности.
Таблица 79 – Исходные данные и вспомогательные расчеты к примеру 6.2
i |
Y |
G |
0,51 |
|
0,29 |
|
1 |
10,13 |
0,22 |
||||
2 |
18,88 |
1,23 |
0,94 |
|
0,29 |
|
3 |
20,94 |
1,81 |
1,05 |
|
0,76 |
|
4 |
22,16 |
1,02 |
1,11 |
|
0,09 |
|
5 |
24,67 |
1,07 |
1,23 |
|
0,16 |
|
6 |
27,57 |
0,67 |
1,38 |
|
0,71 |
|
7 |
40,15 |
0,75 |
2,01 |
|
1,26 |
|
8 |
51,62 |
2,8 |
2,58 |
|
0,22 |
|
9 |
57,71 |
4,9 |
2,89 |
|
2,01 |
|
10 |
66,32 |
4,45 |
3,32 |
|
1,13 |
|
11 |
76,88 |
4,26 |
3,85 |
|
0,41 |
|
12 |
101,65 |
5,31 |
5,09 |
|
0,22 |
|
13 |
115,97 |
6,4 |
5,80 |
|
0,60 |
|
14 |
124,15 |
11,22 |
6,21 |
|
5,01 |
|
15 |
140,98 |
8,66 |
7,06 |
|
1,60 |
|
16 |
169,38 |
13,41 |
8,48 |
|
4,93 |
|
17 |
186,33 |
5,46 |
9,33 |
|
3,87 |
|
18 |
211,78 |
4,79 |
10,60 |
|
5,81 |
|
19 |
261,41 |
18,9 |
13,08 |
|
5,82 |
|
20 |
395,51 |
15,95 |
19,80 |
|
3,85 |
|
|
|
|
|
||||
Примечание: i – номер страны, – ВВП, млрд ден. ед., |
– государственные |
||||||
расходы на здравоохранение, млрд ден. ед., |
|
– модельные значения: |
|
0,05 ∙ |
|||
|
|
|
|
|
|||
Решение
6.2.1 Реализация проверки остаточной компоненты на гетероскедастичность с помощью коэффициента ранговой корреляции Спирмена
Для расчета коэффициента ранговой корреляции Спирмена, необходимо подготовить данные, получив новые переменные:
– ранг независимой переменной для i-й страны;
– рангабсолютногозначенияостаточнойкомпонентые дляi-йстраны. Ранги должны быть стандартизированы, то есть должно выполняться условие
107
|
|
|
1 |
, |
|
2 |
Упорядочивание значений осуществляется по убыванию. Одинаковым значениям рангов присваивается среднее их рангов.
Вспомогательные расчеты сведем в таблицу 80.
Рис. 32 – Отчет Excel Регрессия для решения задачи 6.1
Проверим, является ранжировка стандартизированной, сравнив левую и правую части равенства.
210
, |
1 |
|
|
|
|
2 |
Таким210 |
|
|
|
образом, ранжировка является стандиртизированной. |
Таблица 80 – Вспомогательные расчеты для решения примера 6.2
i |
|
|
| |
е |
| |
|
|
|
|
|
|||||||
1 |
10,13 |
20 |
|
|
0,29 |
15,5 |
4,5 |
20,25 |
2 |
18,88 |
19 |
|
|
0,29 |
15,5 |
3,5 |
12,25 |
3 |
20,94 |
18 |
|
|
0,76 |
11 |
7 |
49,00 |
4 |
22,16 |
17 |
|
|
0,09 |
20 |
-3 |
9,00 |
|
|
|
|
|
|
108 |
|
|
|
|
|
i |
5 |
|
|
|
|
|
|
|
|
| |
е |
| |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
24,67 |
|
|
16 |
|
|
|
|
0,16 |
|
|
|
|
19 |
|
|
-3 |
|
|
9,00 |
|
||||||||
|
|
|
|
6 |
|
27,57 |
|
|
15 |
|
|
|
|
0,71 |
|
|
|
|
12 |
|
|
|
|
|
3 |
|
|
9,00 |
|
|||||
|
|
|
|
7 |
|
40,15 |
|
|
14 |
|
|
|
|
1,26 |
|
|
|
|
9 |
|
|
|
|
|
5 |
|
|
25,00 |
|
|||||
|
|
|
|
8 |
|
51,62 |
|
|
13 |
|
|
|
|
0,22 |
|
|
|
|
18 |
|
|
-5 |
|
|
25,00 |
|
||||||||
|
|
|
|
9 |
|
57,71 |
|
|
12 |
|
|
|
|
2,01 |
|
|
|
|
7 |
|
|
|
|
|
5 |
|
|
25,00 |
|
|||||
|
|
|
|
10 |
|
66,32 |
|
|
11 |
|
|
|
|
1,13 |
|
|
|
|
10 |
|
|
|
|
|
1 |
|
|
1,00 |
|
|||||
|
|
|
|
11 |
|
76,88 |
|
|
10 |
|
|
|
|
0,42 |
|
|
|
|
14 |
|
|
-4 |
|
|
16,00 |
|
||||||||
|
|
|
|
12 |
|
101,65 |
|
|
9 |
|
|
|
|
0,23 |
|
|
|
|
17 |
|
|
-8 |
|
|
64,00 |
|
||||||||
|
|
|
|
13 |
|
115,97 |
|
|
8 |
|
|
|
|
0,60 |
|
|
|
|
13 |
|
|
-5 |
|
|
25,00 |
|
||||||||
|
|
|
|
14 |
|
124,15 |
|
|
7 |
|
|
|
|
5,01 |
|
|
|
|
3 |
|
|
|
|
|
4 |
|
|
16,00 |
|
|||||
|
|
|
|
15 |
|
140,98 |
|
|
6 |
|
|
|
|
1,61 |
|
|
|
|
8 |
|
|
-2 |
|
|
4,00 |
|
||||||||
|
|
|
|
16 |
|
169,38 |
|
|
5 |
|
|
|
|
4,94 |
|
|
|
|
4 |
|
|
|
|
|
1 |
|
|
1,00 |
|
|||||
|
|
|
|
17 |
|
186,33 |
|
|
4 |
|
|
|
|
3,86 |
|
|
|
|
5 |
|
|
-1 |
|
|
1,00 |
|
||||||||
|
|
|
|
18 |
|
211,78 |
|
|
3 |
|
|
|
|
5,80 |
|
|
|
|
2 |
|
|
|
|
|
1 |
|
|
1,00 |
|
|||||
|
|
|
|
19 |
|
261,41 |
|
|
2 |
|
|
|
|
5,83 |
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
1,00 |
|
|||||
|
|
|
|
20 |
|
395,51 |
|
|
1 |
|
|
|
|
3,83 |
|
|
|
|
6 |
|
|
-5 |
|
|
25,00 |
|
||||||||
∑ |
|
|
|
|
|
|
|
|
|
210 |
|
|
|
|
|
|
|
|
|
|
|
210 |
|
|
|
|
|
|
|
|
338,5 |
|
||
|
|
|
|
Рассчитаем значение корреляции между зависимой переменной и |
||||||||||||||||||||||||||||||
остатками по формулам (67): |
1 |
|
|
6 ∙ ∑ |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
,| |
е |
| |
|
|
|
|
|
|
|
|
|
|
(67) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
∙ |
|
1 |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
Проверим |
,|е | |
|
1 |
|
|
|
|
6 ∙ 338,5 |
1 |
|
1 |
|
|
0,255 |
0,745 |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
20 ∙ |
|
20 |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
полученный коэффициент на значимость. |
|
|
|
|||||||||||||||||||||||||
,| |
|
| |
H0: |
в генеральной совокупности |
коэффициент |
Спирмена незначим, |
||||||||||||||||||||||||||||
е |
0 |
, т.е. гетероскедастичности остатков нет. |
|
|
|
|
|
|||||||||||||||||||||||||||
|
,| |
|
| |
H1: в генеральной |
|
|
совокупности |
коэффициент |
Спирмена значим, |
|||||||||||||||||||||||||
е |
0 |
, т.е. гетероскедастичность в остатках есть. |
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
Проверку гипотезц осуществляем с помощью t-критерия Стъюдента |
||||||||||||||||||||||||||||||
|
|
|
|
Расчетное значение критерия примет значение: |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
РАСЧ |
|
,| |
е |
| ∙ |
|
|
|
√ |
|
2е |
|
|
|
0,745 ∙ |
|
√ |
18 |
|
|
4,74 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
,| |
| |
|
|
|
1 |
|
0,555 |
||||||||||||||
|
|
|
|
Табличное значение критерия определяется по таблицам распределения |
||||||||||||||||||||||||||||||
Стъюдента |
КРИТ |
; |
|
2 |
|
|
|
0,05; 18 |
|
|
2,10 |
. Поскольку расчетное значение |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
109 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
больше критического |
| РАСЧ| |
КРИТ |
, принимается конкурирующая гипотеза о |
|||||||||
неслучайности |
|
|
абсолютными |
величинами остатков |
||||||||
|
корреляции |
между |
||||||||||
е и значениями объясняющей переменной |
|
. |
Гипотезу об |
отсутствии |
||||||||
гетероскедастичности остатков отвергаем с |
вероятностью не менее 95%. |
|
||||||||||
|
|
|
|
|
|
|||||||
6.2.2. Реализация |
проверки |
|
остаточной |
компоненты |
на |
|||||||
гетероскедастичность с помощью теста Голфельда-Квандта |
|
|
||||||||||
Воспользуемся |
тестом |
Голдфелда-Квандта |
для |
проверки |
гипотезы |
об |
||||||
отсутствии гетероскедастичности в остатках. Результаты расчетов сведем в таблицу 81.
Предварительная подготовка данных в тесте Голфельда-Квандта состоит в получении упорядоченных значений независимой переменной по возрастанию и разбиении данных на три подвыборки. Если независимых переменных несколько, выбираем наиболее значимую переменную. Обозначим , , – число наблюдений в первой, второй и третьей подвыборках соответственно.
По первой и третьей подвыборкам строится регрессионная модель, имеющая одинаковую спецификацию (значения второй подвыборки
исключаются). |
|
|
|
7 |
|
|
|
6 |
. Построим 20 |
разобьём подвыборки следующим образом : |
, |
||||
При |
|
|
|
|
|||
регрессионные модели по первой и третьей |
подвыборкам, а |
||||||
|
|
|
|
|
|
|
|
значения наблюдений второй подвыборки исключим. Результаты идентификации моделей на подвыборках представлены в таблицах дисперсионного анализа Excel на рисунке 33.
Рис. 33 Таблицы дисперсионного анализа для реализации теста ГолфельдаКванда (пример 6.2)
Сформулируем основную и альтернативную гипотезы:
110
