Таблица 3.8 Матрица корреляции между регрессорами
|
CB |
B_P |
C_R |
I_IMOEX |
INVEST |
M_S |
||
|
CB |
1.000000 |
-0.144243 |
-0.121418 |
0.062049 |
0.341204 |
0.232032 |
|
|
B_P |
-0.144243 |
1.000000 |
-0.045471 |
0.008605 |
0.033104 |
0.108493 |
|
|
C_R |
-0.121418 |
-0.045471 |
1.000000 |
-0.131323 |
-0.407909 |
-0.147593 |
|
|
I_IMOEX |
0.062049 |
0.008605 |
-0.131323 |
1.000000 |
0.148375 |
-0.087931 |
|
|
INVEST |
0.341204 |
0.033104 |
-0.407909 |
0.148375 |
1.000000 |
0.085569 |
|
|
M_S |
0.232032 |
0.108493 |
-0.147593 |
-0.087931 |
0.085569 |
1.000000 |
Построение матрицы корреляций также является необходимым шагом при эконометрическом моделировании, поскольку оно позволяет избежать проблем, связанных с мультиколлинеарностью. Матрица корреляций всех показателей приведена в Таблице 3.8. Максимальный по модулю коэффициент парной корреляции между показателями составляет 0.407, что указывает на отсутствие сильной зависимости между регрессорами и делает возможным включение всех объясняющих переменных в модель. Далее перейдем к построению гистограмм для отслеживания выбросов (рис. 3.2).
Рис. 3.2. Гистограммы по всем переменным
Несмотря на то, что гистограммы указывают на наличие выбросов в переменных инвестиции в основной капитал, индекс ММВБ и объем обращающихся облигаций нефинансовых организаций, их удаление нецелесообразно, так как оно может привести к ограниченности выборки, снижению вариативности данных и искажению коэффициентов модели.
Таким образом, анализ данных можно считать завершенным, после чего можно приступать к построению модели. Оценим регрессию с помощью метода наименьших квадратов. В ходе анализа всех возможных моделей, было выявлено, что данная регрессия обладает наилучшими показателями (табл.3.9).
Стоит отметить, что данная регрессионная модель включает в себя переменную денежная масса, которая преобразована в стационарный временной ряд с помощью метода "детрендизации". Остальные переменные, включая зависимую переменную, были рассчитаны с помощью метода взятия разностей. Кроме того, мы можем наблюдать, что такие переменные, как индекс ММВБ и средневзвешенная ставка по кредитам включены в исследуемую модель с использованием временного лага. Это означает, что значения прошлого периода индекса ММВБ и средневзвешенной ставки по кредитам влияют на текущее значение корпоративных облигаций.
Таблица 3.9 Модель зависимости объема корпоративных облигаций нефинансовых организаций от исследуемых факторов
|
Dependent Variable: CB |
|||||
|
Method: Least Squares |
|||||
|
Sample (adjusted): 2011M03 2017M12 |
|||||
|
Included observations: 82 after adjustments |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
B_P |
-3.238773 |
1.753058 |
-1.847499 |
0.0686 |
|
|
C_R(-1) |
-19.74866 |
18.42862 |
-1.071630 |
0.2873 |
|
|
I_IMOEX(-1) |
0.144389 |
0.067839 |
2.128413 |
0.0365 |
|
|
INVEST |
0.065916 |
0.019253 |
3.423704 |
0.0010 |
|
|
M_S |
0.028251 |
0.011672 |
2.420368 |
0.0179 |
|
|
C |
42.93313 |
9.697179 |
4.427384 |
0.0000 |
|
|
R-squared |
0.258514 |
Mean dependent var |
45.42195 |
||
|
Adj. R-square |
0.209732 |
S.D. dependent var |
98.27759 |
||
|
S.E. of regression |
87.36587 |
Akaike info criterion |
11.84844 |
||
|
Sum Sq. resid |
580092.5 |
Schwarz criterion |
12.02454 |
||
|
Log likelihood |
-479.7861 |
F-statistic |
5.299363 |
||
|
Durbin-Watson stat |
1.870002 |
Prob(F-statistic) |
0.000314 |
Начнем с проверки исследуемой регрессии на значимость. Очевидно, что при уровне значимости б = 5%, Prob.(Fst) = 0.0003 < 0.05, значит данная регрессионная модель является значимой. Тест на значимость регрессии в целом показал, что все регрессоры оказывают влияние на объем обращающихся корпоративных облигаций нефинансового сектора.
Что касается значимости отдельных коэффициентов, то на уровне значимости 5 % такие переменные, как денежная масса, цена на нефть, индекс ММВБ, инвестиции в основной капитал, а также константа являются значимыми для объема корпоративных облигаций : при б = 0.05, б > Prob.(t-st), значит H0 - отклоняется. Регрессор - средневзвешенная ставка по кредитам не является значимым на любом возможном уровне значимости.
Кроме того, величина коэффициента детерминации данной регрессионной модели составляет 0,2585, и это указывает на то, что практически 25,9 % изменения зависимой переменной объясняется включенными в модель факторами.
Исходная регрессия обладает наилучшими показателями среди других, в связи с чем было принято решение выбрать ее для дальнейшего анализа.
На следующем этапе стоит проверить выбранную для исследования модель на наличие гетероскедастичности. Для начала стоит сделать предположение о наличии гетероскедастичности в данных на основе графика зависимости между остатками модели и регрессорами (рис. 3.3).
Рис. 3.3 Диаграммы рассеивания. Зависимость между остатками модели и ее регрессорами.
На выше представленных графиках мы можем наблюдать неоднородные данные, что может являться признаком гетероскедастичности в модели. Для того чтобы протестировать данную регрессионную модель на наличие гетероскедастичности, необходимо воспользоваться тестом Уайта, который представлен в таблице 3.10.
Таблица 3.10 Тест Уайта на наличие гетероскедастичности
|
White Heteroskedasticity Test: |
|||||
|
F-statistic |
8.964273 |
Probability |
0.000000 |
||
|
Obs*R-squared |
45.75809 |
Probability |
0.000002 |
||
|
Test Equation: |
|||||
|
Dependent Variable: RESID^2 |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/23/18 Time: 20:55 |
|||||
|
Sample: 2011M03 2017M12 |
|||||
|
Included observations: 82 |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
C |
463.3471 |
2593.925 |
0.178628 |
0.8587 |
|
|
B_P |
29.20111 |
355.1276 |
0.082227 |
0.9347 |
|
|
B_P^2 |
186.9152 |
41.66938 |
4.485673 |
0.0000 |
|
|
C_R(-1) |
628.3334 |
3502.328 |
0.179405 |
0.8581 |
|
|
C_R(-1)^2 |
-3021.731 |
3488.136 |
-0.866288 |
0.3893 |
|
|
I_IMOEX(-1) |
13.85578 |
13.89890 |
0.996898 |
0.3222 |
|
|
I_IMOEX(-1)^2 |
-0.022980 |
0.038873 |
-0.591167 |
0.5563 |
|
|
INVEST |
24.69961 |
5.529842 |
4.466603 |
0.0000 |
|
|
INVEST^2 |
0.015437 |
0.004317 |
3.576008 |
0.0006 |
|
|
M_S |
2.347331 |
2.400184 |
0.977979 |
0.3314 |
|
|
M_S^2 |
-0.003085 |
0.001866 |
-1.653506 |
0.1026 |
|
|
R-squared |
0.558025 |
Mean dependent var |
7074.298 |
||
|
Adj. R-squared |
0.495776 |
S.D. dependent var |
21563.27 |
||
|
S.E. of regression |
15311.81 |
Akaike info criterion |
22.23489 |
||
|
Sum Sq. resid |
1.66E+10 |
Schwarz criterion |
22.55774 |
||
|
Log likelihood |
-900.6305 |
F-statistic |
8.964273 |
||
|
Durbin-Watson stat |
2.302349 |
Prob(F-statistic) |
0.000000 |
||
H0 - модель гомоскедастична, Ha- модель гетероскедастична
Как мы видим, Prob.(F-st.) = 0.0000, 0.0000 < 0.05, следовательно, гипотеза H0 о наличии гомоскедастичности отклоняется, из чего можно сделать окончательный вывод, что модель является гетероскедастичной.
Следующей этап включает в себя тестирование исследуемой регрессионной модели на наличие автокорреляции.
Для начала, можно отследить наличие автокорреляции в данных на основе коррелограмы для остатков регрессии (рис. 3.4).
Рис. 3.4 Коррелограма остатков модели
Согласно представленному рисунку, мы можем предположить об отсутствии автокорреляции в данных, о чем свидетельствуют значения автокорреляционной функции, находящиеся в пределах пунктира. Однако, стоит отметить, что одно из значений автокорреляционной функции находятся на границе допустимой области, что может указывать на наличие автокорреляции 11 порядка.
Следовательно, стоит воспользоваться тестом Бреуша Годфри, чтобы проверить выдвинутые предположения ( табл. 3.11).
Таблица № 3.11 Тест Бреуша Годфри
|
Breusch-Godfrey Serial Correlation LM Test: |
|||||
|
F-statistic |
0.951320 |
Probability |
0.498695 |
||
|
Obs*R-squared |
11.37078 |
Probability |
0.412743 |
||
|
Test Equation: |
|||||
|
Dependent Variable: RESID |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/19/18 Time: 11:25 |
|||||
|
Presample missing value lagged residuals set to zero. |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
C_R(-1) |
26.45087 |
22.57700 |
1.171585 |
0.2456 |
|
|
B_P |
0.066184 |
1.854714 |
0.035684 |
0.9716 |
|
|
I_IMOEX(-1) |
0.005681 |
0.070806 |
0.080231 |
0.9363 |
|
|
INVEST |
-0.008672 |
0.021017 |
-0.412649 |
0.6812 |
|
|
M_S |
0.010372 |
0.013505 |
0.768018 |
0.4453 |
|
|
C |
3.659996 |
9.841267 |
0.371903 |
0.7112 |
|
|
RESID(-1) |
-0.025441 |
0.134689 |
-0.188887 |
0.8508 |
|
|
RESID(-2) |
-0.095203 |
0.134646 |
-0.707063 |
0.4821 |
|
|
... |
... |
... |
... |
... |
|
|
RESID(-11) |
0.182323 |
0.130456 |
1.397589 |
0.1670 |
H0: с=0 (автокорреляции нет)
Ha : с ? 0
Как мы видим, Prob.(Fst) > 0.05, Prob.(Chi-Square) > 0.05, значит гипотеза Н0 об отсутствии автокорреляции подтверждается, из чего можно заключить, что в модели отсутствует автокорреляция одиннадцатого порядка.
Кроме того, согласно данной регрессии мы можем наблюдать, что в модели отсутствует зависимость текущих значений остатков модели от остатков в предыдущих периодах, так как переменные RESID01(-1) и RESID01(-2), RESID01(-11) не являются значимыми (0.850 > 0.05, 0.0482 > 0.05 и 0.167 > 0.05 ).
Таким образом, достаточно будет применить поправки Уайта( табл. 3.12)., которые помогут избавиться от неблагоприятных последствий гетероскедастичности и скорректировать стандартные ошибки регрессии.
Таблица 3.12 Применение поправок White
|
Dependent Variable: CB |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/19/18 Time: 12:00 |
|||||
|
Sample (adjusted): 2011M03 2017M12 |
|||||
|
Included observations: 82 after adjustments |
|||||
|
White Heteroskedasticity-Consistent Standard Errors & Covariance |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
B_P |
-3.238773 |
3.506200 |
-0.923727 |
0.3586 |
|
|
C_R(-1) |
-19.74866 |
12.89743 |
-1.531209 |
0.1299 |
|
|
I_IMOEX(-1) |
0.144389 |
0.066055 |
2.185904 |
0.0319 |
|
|
INVEST |
0.065916 |
0.028652 |
2.300564 |
0.0242 |
|
|
M_S |
0.028251 |
0.011268 |
2.507151 |
0.0143 |
|
|
C |
42.93313 |
8.903950 |
4.821808 |
0.0000 |
|
|
R-squared |
0.258514 |
Mean dependent var |
45.42195 |
||
|
Adj. R-squared |
0.209732 |
S.D. dependent var |
98.27759 |
||
|
S.E. of regression |
87.36587 |
Akaike info criterion |
11.84844 |
||
|
Sum Sq. resid |
580092.5 |
Schwarz criterion |
12.02454 |
||
|
Log likelihood |
-479.7861 |
F-statistic |
5.299363 |
||
|
Durbin-Watson stat |
1.870002 |
Prob(F-statistic) |
0.000314 |
||