PERIOD^2
0.011310
0.001348
8.387553
0.0000
LOG(EMISSION_VOL)
-0.470637
0.072219
-6.516825
0.0000
RATING_INTER
0.159979
0.033468
4.780001
0.0000
TYPE_BOND
-0.477368
0.215455
-2.215634
0.0271
TYPE_ORGANISER
-0.939571
0.233826
-4.018254
0.0001
CALL
0.942421
0.192854
4.886718
0.0000
C
10.21405
0.508146
20.10064
0.0000
R-squared
0.477761
Mean dependent var
2.367726
Adjusted R-squared
0.472313
S.D. dependent var
3.061158
S.E. of regression
2.223689
Akaike info criterion
4.447924
Sum squared resid
3317.955
Schwarz criterion
4.501186
Log likelihood
-1502.070
F-statistic
87.69313
Durbin-Watson stat
1.641631
Prob(F-statistic)
0.000000
Согласно построенной модели, принятые меры способствовали улучшению качества исследуемой регрессии. Переменные PERIOD^2 и LOG(EMISSION_VOL) оказались значимыми для премии за риск корпоративных облигаций на уровне значимости 1%. Стоит отметить, что коэффициент детерминации данной модели значительно выше, чем в прошлой модели и составляет 47,8%.
После того, как была выбрана наилучшая функциональная взаимосвязь между премией за риск и характеристиками облигационного займа, можно приступить к рассмотрению макроэкономических факторов.
На следующем этапе дополним получившуюся регрессионную модель макроэкономическими факторами (таблица 3.3).
Таблица 3.3
Модель зависимости премии за риск от характеристик облигационного займа и макроэкономических факторов
|
Dependent Variable: SPREAD |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/05/20 Time: 13:14 |
|||||
|
Sample: 1 782 IF SPREAD<29 AND DEFAULT_RATE<0.3 AND |
|||||
|
EMISSION_VOL<175000 AND PERIOD<47.5 |
|||||
|
Included observations: 679 |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
PERIOD |
-0.415512 |
0.047507 |
-8.746284 |
0.0000 |
|
|
PERIOD^2 |
0.009990 |
0.001335 |
7.483529 |
0.0000 |
|
|
LOG(EMISSION_VOL) |
-0.483666 |
0.071259 |
-6.787400 |
0.0000 |
|
|
RATING_INTER |
0.161486 |
0.032643 |
4.947070 |
0.0000 |
|
|
TYPE_BOND |
-0.830097 |
0.222280 |
-3.734469 |
0.0002 |
|
|
TYPE_ORGANISER |
-0.932615 |
0.234005 |
-3.985455 |
0.0001 |
|
|
CALL |
0.753102 |
0.191086 |
3.941177 |
0.0001 |
|
|
BRENT_GROWTH |
0.024470 |
0.011547 |
2.119068 |
0.0345 |
|
|
EXCH_RATE |
-0.155754 |
0.103432 |
-1.505857 |
0.1326 |
|
|
GDP_GROWTH |
-0.083946 |
0.015470 |
-5.426464 |
0.0000 |
|
|
C |
11.03468 |
0.528982 |
20.86024 |
0.0000 |
|
|
R-squared |
0.506130 |
Mean dependent var |
2.367726 |
||
|
Adjusted R-squared |
0.498737 |
S.D. dependent var |
3.061158 |
||
|
S.E. of regression |
2.167297 |
Akaike info criterion |
4.400906 |
||
|
Sum squared resid |
3137.714 |
Schwarz criterion |
4.474142 |
||
|
Log likelihood |
-1483.108 |
F-statistic |
68.45838 |
||
|
Durbin-Watson stat |
1.722750 |
Prob(F-statistic) |
0.000000 |
Как мы видим, включение в модель макро-факторов оказало положительное влияние на исследуемую регрессию. Стандартные ошибки уменьшились, коэффициент детерминации увеличился и достиг 50,6%.
Далее перейдем к рассмотрению отраслевых переменных. На основе последней модели проведем тест Чоу, для того, чтобы понять каким образом следует включать в модель отраслевые переменные.
Данный тест проводится для выявления структурных сдвигов в регрессионной модели, что позволит выяснить, есть ли различие между той или иной отраслью. Проведем тест отдельно для каждой отраслевой переменной. Таким образом, будет проведено 8 тестов. Рассмотрим тест Чоу на примере нефтегазовой отрасли.
Сначала необходимо отсортировать значения бинарной переменной OIL_GAS от 0 к 1. Затем можно приступать к тестированию. Переходная точка от 0 к 1 находится в 653 наблюдении. Данное тестирование проводится на основе скорректированной модели.
Таблица 3.4
Тест Чоу для нефтегазовой отрасли
|
Chow Breakpoint Test: 653 |
||||
|
F-statistic |
1.369869 |
Probability |
0.182463 |
|
|
Log likelihood ratio |
15.39385 |
Probability |
0.165167 |
H0: структурные изменения отсутствуют
Hа: есть структурные изменения в модели
Благодаря данному тесту было выявлено, что Prob (Fst) = 0,18 (0,18 > 0,05). Следовательно, нет оснований отклонять гипотезу H0, из чего следует, что в модели отсутствуют структурные изменения и нет необходимости в построении отдельной регрессии для нефтегазовой отрасли. Аналогичным образом был проведен тест Чоу для остальных отраслевых переменных. Согласно проведенным тестам, было выявлено, что построение отдельных регрессий для каждой отрасли не требуется. Результаты проведенного тестирования для всех исследуемых отраслей отражены в Приложении 3.
Таким образом, мы можем дополнить имеющуюся модель отраслевыми переменными. Среди восьми отраслевых переменных было выбрано пять основных переменных, которым принадлежит наибольшая часть рынка обращающихся корпоративных облигаций нефинансового сектора: нефтегазовая отрасль, строительная отрасль, энергетическая отрасль, отрасль торговли и транспортная отрасль. Вместе с тем, было принято решение не использовать в регрессии нефтегазовую отрасль, поскольку в модели уже присутствует такой макроэкономический фактор, как цена на нефть, который содержит в себе тенденции нефтяного рынка. Остальные переменные были включены в модель. Используемые в модели отрасли охватывают более 50% рынка обращающихся корпоративных облигаций.
Получившаяся регрессионная модель отражена в таблице 3.5.
Таблица 3.5
Модель зависимости премии за риск от характеристик облигационного займа, макроэкономических и отраслевых факторов
|
Dependent Variable: SPREAD |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/05/20 Time: 13:15 |
|||||
|
Sample: 1 782 IF SPREAD<29 AND DEFAULT_RATE<0.3 AND |
|||||
|
EMISSION_VOL<175000 AND PERIOD<47.5 |
|||||
|
Included observations: 679 |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
PERIOD |
-0.441471 |
0.047522 |
-9.289808 |
0.0000 |
|
|
PERIOD^2 |
0.010714 |
0.001340 |
7.994663 |
0.0000 |
|
|
LOG(EMISSION_VOL) |
-0.487893 |
0.070035 |
-6.966451 |
0.0000 |
|
|
RATING_INTER |
0.121116 |
0.032973 |
3.673229 |
0.0003 |
|
|
TYPE_BOND |
-0.946508 |
0.219507 |
-4.311964 |
0.0000 |
|
|
TYPE_ORGANISER |
-0.856934 |
0.230738 |
-3.713881 |
0.0002 |
|
|
CALL |
0.904562 |
0.189213 |
4.780664 |
0.0000 |
|
|
BRENT_GROWTH |
0.021477 |
0.011295 |
1.901481 |
0.0577 |
|
|
EXCH_RATE |
-0.172856 |
0.101258 |
-1.707086 |
0.0883 |
|
|
GDP_GROWTH |
-0.074462 |
0.015228 |
-4.889717 |
0.0000 |
|
|
BUILDING |
1.127356 |
0.247447 |
4.555948 |
0.0000 |
|
|
ENERGETICS |
-0.716467 |
0.298522 |
-2.400047 |
0.0167 |
|
|
TRADE_RETAIL |
0.534650 |
0.314685 |
1.698998 |
0.0898 |
|
|
TRANSPORT |
0.669979 |
0.281557 |
2.379548 |
0.0176 |
|
|
C |
10.76037 |
0.532757 |
20.19752 |
0.0000 |
|
|
R-squared |
0.531664 |
Mean dependent var |
2.367726 |
||
|
Adjusted R-squared |
0.521789 |
S.D. dependent var |
3.061158 |
||
|
S.E. of regression |
2.116876 |
Akaike info criterion |
4.359603 |
||
|
Sum squared resid |
2975.492 |
Schwarz criterion |
4.459470 |
||
|
Log likelihood |
-1465.085 |
F-statistic |
53.84178 |
||
|
Durbin-Watson stat |
1.724360 |
Prob(F-statistic) |
0.000000 |
Представленная регрессия в целом значима и имеет наименьшее значение стандартных ошибок, и это свидетельствует о том, что прогнозное значение данной регрессии в наименьшей степени отлично от истинного значения по сравнению с другими моделями. Кроме того, положительным аспектом является то, что значение суммы квадратов остатков также уменьшилось по сравнению с исходной регрессией.
Более того, данная модель характеризуется наивысшим значением коэффициента детерминации (53%), следовательно, 53% изменений зависимой переменной объясняются включенными в модель факторами.
Следующий этап состоит в проверке выбранной регрессионной модели на наличие гетероскедастичности и автокорреляции.
Построим график зависимости остатков модели от исследуемых регрессоров с целью выдвинуть предположение о наличии гетероскедастичности.
Рис. 3.9 Графики функциональной связи остатков модели и регрессоров
Представленные диаграммы рассеивания указывают на неоднородность данных, что может свидетельствовать о наличии гетероскедастичности в модели. Проверим данное предположение с помощью теста Уайта, который представлен в таблице 3.6.
Таблица 3.6
Тест Уайта (проверка гетероскедастичности)
|
White Heteroskedasticity Test: |
|||||
|
F-statistic |
13.72084 |
Probability |
0.000000 |
||
|
Obs*R-squared |
199.8345 |
Probability |
0.000000 |
||
|
Test Equation: |
|||||
|
Dependent Variable: RESID^2 |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/05/20 Time: 13:34 |
|||||
|
Sample: 1 782 IF SPREAD<29 AND DEFAULT_RATE<0.3 AND |
|||||
|
EMISSION_VOL<175000 AND PERIOD<47.5 |
|||||
|
Included observations: 679 |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
C |
59.20257 |
4.830787 |
12.25526 |
0.0000 |
|
|
PERIOD |
0.841338 |
0.497226 |
1.692064 |
0.0911 |
|
|
PERIOD^2 |
-0.015570 |
0.027234 |
-0.571713 |
0.5677 |
|
|
(PERIOD^2)^2 |
4.72E-06 |
1.39E-05 |
0.340204 |
0.7338 |
|
|
LOG(EMISSION_VOL) |
-12.19316 |
1.310625 |
-9.303322 |
0.0000 |
|
|
(LOG(EMISSION_VOL))^2 |
0.637209 |
0.096256 |
6.619926 |
0.0000 |
|
|
RATING_INTER |
0.359396 |
0.679920 |
0.528585 |
0.5973 |
|
|
RATING_INTER^2 |
0.021115 |
0.089125 |
0.236917 |
0.8128 |
|
|
TYPE_BOND |
-2.603202 |
1.162037 |
-2.240205 |
0.0254 |
|
|
TYPE_ORGANISER |
0.295364 |
1.254838 |
0.235380 |
0.8140 |
|
|
CALL |
-3.430008 |
1.038529 |
-3.302755 |
0.0010 |
|
|
BRENT_GROWTH |
0.032268 |
0.064970 |
0.496664 |
0.6196 |
|
|
BRENT_GROWTH^2 |
-0.000371 |
0.005131 |
-0.072247 |
0.9424 |
|
|
EXCH_RATE |
0.650604 |
0.548865 |
1.185364 |
0.2363 |
|
|
EXCH_RATE^2 |
0.173266 |
0.185816 |
0.932459 |
0.3514 |
|
|
GDP_GROWTH |
-0.016952 |
0.244605 |
-0.069305 |
0.9448 |
|
|
GDP_GROWTH^2 |
-0.005070 |
0.007717 |
-0.656971 |
0.5114 |
|
|
BUILDING |
-2.956664 |
1.337496 |
-2.210596 |
0.0274 |
|
|
ENERGETICS |
0.115731 |
1.634983 |
0.070784 |
0.9436 |
|
|
TRADE_RETAIL |
-1.915593 |
1.700131 |
-1.126733 |
0.2603 |
|
|
TRANSPORT |
-3.065657 |
1.542884 |
-1.986965 |
0.0473 |
|
|
R-squared |
0.294307 |
Mean dependent var |
4.382168 |
||
|
Adjusted R-squared |
0.272857 |
S.D. dependent var |
13.00641 |
||
|
S.E. of regression |
11.09091 |
Akaike info criterion |
7.680569 |
||
|
Sum squared resid |
80939.48 |
Schwarz criterion |
7.820382 |
||
|
Log likelihood |
-2586.553 |
F-statistic |
13.72084 |
||
|
Durbin-Watson stat |
1.876123 |
Prob(F-statistic) |
0.000000 |
Гипотеза: H0 - модель является гомоскедастичной, Ha - модель является гетероскедастичной.
Согласно проведенному тесту, Prob (F-st) = 0 (0 < 0.05), следовательно, нулевая гипотеза о гомоскедастичности модели отклоняется, что подтверждает выдвинутые предположения о наличии гетероскедастичности.
Далее проверим исследуемую модель на автокорреляцию. Построим коррелограму остатков модели для того, чтобы выдвинуть предположение о наличии автокорреляции в данных (рис. 3.10).
Рис. 3.10 Коррелограма остатков регрессионной модели
Результаты построения коррелограммы свидетельствуют о наличии автокорреляции в исследуемой модели, что подтверждает автокорреляционная функция, выходящая за пределы пунктирной линии.
Для того, чтобы протестировать гипотезу о наличии автокорреляции необходимо использовать тест Бреуша Годфри, который представлен в таблице 3.7.
Таблица 3.7
Тестирование Бреуша Годфри
|
Breusch-Godfrey Serial Correlation LM Test: |
|||||
|
F-statistic |
11.09855 |
Probability |
0.000912 |
||
|
Obs*R-squared |
11.17925 |
Probability |
0.000827 |
||
|
Test Equation: |
|||||
|
Dependent Variable: RESID |
|||||
|
Method: Least Squares |
|||||
|
Date: 05/05/20 Time: 13:58 |
|||||
|
Presample and interior missing value lagged residuals set to zero. |
|||||
|
Variable |
Coefficient |
Std. Error |
t-Statistic |
Prob. |
|
|
PERIOD |
0.008694 |
0.047237 |
0.184045 |
0.8540 |
|
|
PERIOD^2 |
-0.000297 |
0.001333 |
-0.222941 |
0.8237 |
|
|
LOG(EMISSION_VOL) |
-0.015736 |
0.069668 |
-0.225870 |
0.8214 |
|
|
RATING_INTER |
0.001366 |
0.032727 |
0.041735 |
0.9667 |
|
|
TYPE_BOND |
0.034594 |
0.218104 |
0.158614 |
0.8740 |
|
|
TYPE_ORGANISER |
0.008562 |
0.229018 |
0.037388 |
0.9702 |
|
|
CALL |
-0.021140 |
0.187897 |
-0.112509 |
0.9105 |
|
|
BRENT_GROWTH |
6.24E-05 |
0.011210 |
0.005569 |
0.9956 |
|
|
EXCH_RATE |
0.009795 |
0.100539 |
0.097427 |
0.9224 |
|
|
GDP_GROWTH |
-0.001123 |
0.015118 |
-0.074271 |
0.9408 |
|
|
BUILDING |
-0.080083 |
0.246760 |
-0.324539 |
0.7456 |
|
|
ENERGETICS |
0.035476 |
0.296469 |
0.119661 |
0.9048 |
|
|
TRADE_RETAIL |
-0.043741 |
0.312595 |
-0.139927 |
0.8888 |
|
|
TRANSPORT |
0.015801 |
0.279480 |
0.056536 |
0.9549 |
|
|
C |
0.088635 |
0.529420 |
0.167419 |
0.8671 |
|
|
RESID(-1) |
0.135604 |
0.040704 |
3.331449 |
0.0009 |
|
|
R-squared |
0.016464 |
Mean dependent var |
1.84E-15 |
||
|
Adjusted R-squared |
-0.005788 |
S.D. dependent var |
2.094906 |
||
|
S.E. of regression |
2.100960 |
Akaike info criterion |
4.345948 |
||
|
Sum squared resid |
2926.503 |
Schwarz criterion |
4.452472 |
||
|
Log likelihood |
-1459.449 |
F-statistic |
0.739904 |
||
|
Durbin-Watson stat |
2.021145 |
Prob(F-statistic) |
0.744457 |