Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ИСПОЛЬЗОВАНИЕ КОИНТЕГРАЦИОННОГО АНАЛИЗА ПРИ ПОСТРОЕНИИ CИСТЕМЫ ОДНОВРЕМЕННЫХ УРАВНЕНИЙ

Суханова Е.И. 1 Ширнаева С.Ю. 1
1 ФГБОУ ВПО «Самарский государственный экономический университет»
В статье предложена методика отбора экзогенных и эндогенных переменных в систему одновременных уравнений с использованием коинтеграционного анализа временных рядов показателей, отражающих стабилизационные процессы экономики России. Предварительно временные ряды показателей были проверены на устойчивость и на причинно-следственную зависимость. Проверялась гипотеза о наличии коинтеграции временных рядов методом Энгла‒Грэнжера, в соответствии с которым определен порядок интегрируемости временных рядов экономических показателей I(1) (на основе выводов, полученных по автокорреляционной, частной автокорреляционной функциям и расширенному тесту Дики‒Фуллера); протестированы остатки парной регрессии между соответствующими переменными с использованием расширенного теста Дики‒Фуллера. Было получено, что рассмотренные пары временных рядов экономических показателей коинтегрированы. По исходным данным построена эконометрическая модель взаимосвязей экономических показателей, отражающих стабилизационные процессы экономики России, в виде системы одновременных уравнений. Полученная модель может быть использована для построения точечных и интервальных прогнозов рассмотренных экономических показателей.
временные ряды
стационарность
коинтеграция
система одновременных уравнений
стабилизационные процессы
1. Краткосрочные экономические показатели Российской Федерации. – М., 2013. – URL: http://www.gks.ru/bgd/regl/b13_03/Main.htm (дата обращения: 19.03.2013).
2. Магнус Я.Р., Катышев П.К., Пересецкий А.А. Эконометрика. Начальный курс. – М.: Дело, 2004. – 575 с.
3. Суханова Е.И., Ширнаева С.Ю. Построение интервальных прогнозов эндогенных переменных одной эконометрической модели // Вестник Самарского государственного экономического университета. – Самара, 2013. – № 2(100). – С. 109–114.
4. Суханова Е.И., Ширнаева С.Ю. Статистические методы исследования макроэкономической стабилизации в России // Научные исследования Самарского государственного экономического университета: итоги и перспективы: сб. науч. тр.: в 2 ч. – Самара: Изд-во Самар. гос. экон. ун-та, 2011. – Ч.1. – С. 83–90.
5. Суханова Е.И., Ширнаева С.Ю. Типологизация показателей, отражающих стабилизационные процессы экономики России, по критериям устойчивости // Вестник Самарского государственного экономического университета. – Самара, 2012. – № 5(91). – С. 103–109.
6. Эконометрика: учебник / под ред. И.И. Елисеевой. – М.: Финансы и статистика, 2002. – 344 с.
7. Engle R.F., Granger C.W.J. Cointegration and error correction: representation, estimation and testing // Econometrica. 1987. Vol. 55, pp. 251–276.
8. MacKinnon J.G. Critical values for cointegration tests // UC San Diego Discussion Paper, 1990, pp. 90–94.

При построении эконометрической модели в виде системы одновременных уравнений базовой задачей является отбор экзогенных и эндогенных переменных в модель. Обоснованный подход к решению данной задачи позволяет получить адекватную и качественную эконометрическую модель и дает возможность использовать ее для анализа и прогнозирования. Деление переменных на эндогенные и экзогенные определяется, прежде всего, содержательной стороной модели [2]. Требования, предъявляемые к набору переменных, вытекают из целей моделирования и сущности изучаемого экономического явления или процесса.

При моделировании такого макроэкономического явления, как стабилизационные процессы экономики России, выбор системы одновременных уравнений (СОУ) в качестве модели позволит провести всесторонний анализ взаимосвязей показателей, отражающих данные процессы. Если при этом переменные (показатели) представлены в виде временных рядов, то исследование этих рядов на наличие долговременной зависимости и определение направления этой зависимости становится неотъемлемой частью эконометрического моделирования стабилизационных процессов.

В ходе исследования был сформирован информационный массив, включающий ежемесячные статистические данные по 35 экономическим показателям России. Объем выборки составил 169 наблюдений за период с января 1999 года по январь 2013 года [1]. Обработка статистических данных проводилась с использованием пакетов прикладных программ Microsoft Excel и Econometric Views.

Сначала было проведено исследование устойчивости временных рядов всех рассмотренных показателей. В качестве критериев устойчивости использовались индекс устойчивости уровней, коэффициент устойчивости и ранговый коэффициент корреляции Спирмена [5].

При отборе эндогенных переменных (проверка устойчивости показала, что временные ряды всех отобранных эндогенных переменных устойчивы как по степени устойчивости уровней, так и по наличию необходимой тенденции изменения) использовалось определение стабилизации [4]. Под стабилизацией экономики понимается процесс приведения экономической системы в стабильное состояние, для которого характерен ряд признаков. Основные признаки стабилизации экономики и экономические показатели, поведение которых отражает соответствующие признаки, представлены в табл. 1.

Таблица 1

Экономические показатели (эндогенные переменные), отражающие стабилизационные процессы экономики России

Признак стабилизации

Экономический показатель

Относительно незначительный, но устойчивый рост производства

Индекс промышленного производства, % к соответствующему периоду предшествующего года (Y1)

Рост занятости, соответствующий естественному росту населения

Общая численность безработных (на конец периода), млн чел. (Y2)

Сбалансированность внешнеторговых операций

Чистый экспорт, млрд долл. США (Y3)

Практическая неизменность цен

Индекс потребительских цен, % к соответствующему периоду предшествующего года (Y4)

Стабильный уровень благосостояния населения

Начисленная среднемесячная заработная плата одного работника номинальная, руб. (Y5)

Сокращение бюджетного дефицита

Инвестиции в основной капитал, млрд руб. (Y6)

К временным рядам всех показателей (сезонная составляющая была удалена) был применен тест Грэнжера на причинно-следственную зависимость [2], позволивший определить направление причинно-следственных связей рассмотренных экономических показателей. Следует отметить, что данный тест является чувствительным к количеству лагов, поэтому он проводился для величины лага от 2 до 6 месяцев. Учитывая результаты исследования устойчивости временных рядов, а также результаты применения теста Грэнжера, в качестве экзогенных переменных были отобраны следующие показатели: коммерческий грузооборот транспорта, млрд т-км (X1); погрузка грузов на железнодорожном транспорте, млн т (X2); объем работ по виду деятельности «Строительство», млрд руб. (X3); официальный курс доллара, руб. за 1 долл. США (X4); оборот розничной торговли, млрд руб. (X5); объем платных услуг населению, млрд руб. (X6); денежные доходы в среднем на душу населения, руб. (X7); просроченная кредиторская задолженность организаций в бюджет, млрд руб. (X8); просроченная дебиторская задолженность организаций, млрд руб. (X9); средние цены производителей на нефть, руб. за 1 т (X10); средние цены производителей на газ горючий, руб. за 1 тыс. м3 (X11).

Характерной особенностью временных рядов большинства показателей, отражающих стабилизационные процессы экономики России, является их нестационарность. Наличие в данных рядах тренда наряду с наличием сезонной составляющей, при применении к ним традиционных методов корреляционного и регрессионного анализа может привести к явлениям ложной корреляции и мнимой регрессии [2]. С целью применения этих методов временные ряды рассмотренных показателей были проверены на коинтеграцию.

Под коинтеграцией будем понимать зависимость в уровнях двух (или более) временных рядов, которая выражается в совпадении или противоположной направленности их тенденций и случайной колеблемости [6]. Нестационарные временные ряды Xtи Yt называются коинтегрированными, если существует линейная комбинация этих рядов, представляющая собой стационарный временной ряд.

Для проверки наличия коинтеграции двух временных рядов Xtи Yt в данной работе использовался метод Энгла–Грэнжера [7], в соответствии с которым исследование проводилось по следующей схеме:

1) определялся порядок интегрированности данных временных рядов;

2) с помощью метода наименьших квадратов оценивалось уравнение регрессии Yt = a + bXt + εt (называемое уравнением коинтеграции);

3) с помощью теста Дики–Фуллера остатки регрессии et данного уравнения исследовались на стационарность; основная гипотеза состояла в том, что Xtи Yt не коинтегрированы;

4) если по результатам теста Дики–Фуллера остатки et оказывались стационарными, то основная гипотеза отвергалась, а временные ряды Xtи Yt являлись коинтегрированными.

При применении метода Энгла–Грэнжера использовались специально разработанные таблицы критических точек [8].

Коинтегрированными могут быть только такие два временных ряда, которые имеют одинаковый порядок интегрированности. Зная порядок интегрированности временных рядов рассматриваемых экономических показателей, можно подобрать преобразование временного ряда, которое обеспечит его стационарность. Это, в свою очередь, делает применение эконометрических методов корректным и позволяет избежать мнимой регрессии.

Определение порядка интегрированности временных рядов переменных рассмотрим на примере экзогенной переменной X1 – коммерческий грузооборот транспорта (млрд т-км). Динамика изменения временного ряда данной переменной (сезонная составляющая удалена) представлена на рис. 1.

pic_77.tif

Рис. 1. Динамика коммерческого грузооборота транспорта (X1) за период с января 1999 года по январь 2013 года

Анализ динамики позволяет предположить, что исходный временной ряд не является стационарным (имеет место возрастающий линейный тренд). Для подтверждения этого предположения были построены автокорреляционная (АКФ) и частная автокорреляционная (ЧАКФ) функции данного временного ряда (рис. 2).

Autocorrelation

Partial Correlation

 

AC

PAC

Q-Stat

Prob

.|*******

.|*******

1

0,968

0,968

160,11

0,000

.|*******

.|. |

2

0,939

0,043

311,80

0,000

.|*******

.|. |

3

0,914

0,042

456,34

0,000

.|******|

.|. |

4

0,889

–0,006

593,87

0,000

.|******|

.|. |

5

0,861

–0,052

723,73

0,000

.|******|

.|. |

6

0,832

–0,044

845,64

0,000

.|******|

.|. |

7

0,805

0,018

960,54

0,000

.|******|

.|. |

8

0,781

0,026

1069,3

0,000

.|***** |

.|. |

9

0,761

0,064

1173,2

0,000

.|***** |

*|. |

10

0,734

–0,097

1270,6

0,000

.|***** |

.|. |

11

0,711

0,034

1362,6

0,000

.|***** |

.|. |

12

0,691

0,012

1449,9

0,000

Рис. 2. Автокорреляционная и частная автокорреляционная функции временного ряда переменной X1 (фрагмент)

Анализ рис. 2 показал, что АКФ показателя X1 (коммерческий грузооборот транспорта) экспоненциально убывает, а ЧАКФ имеет наибольшее сильно выделяющееся значение коэффициента автокорреляции первого порядка и не имеет значимых коэффициентов автокорреляции других порядков. Это означает, что временной ряд данного показателя представляет собой модель авторегрессии первого порядка AR(1), то есть не является стационарным.

Еще одним способом проверки временного ряда на стационарность является расширенный тест Дики–Фуллера (ADF-тест) (используется распределение t-статистики, описанное Дики и Фуллером), который был применен для проверки основной гипотезы о нестационарности временного ряда. Вычисленное наблюдаемое значение t-статистики ADF-теста, а также критические значения представлены в табл. 2.

Из табл. 2 видно, что основная гипотеза не отвергается на рассмотренных уровнях значимости (1, 5 и 10 %). Таким образом, как графики АКФ и ЧАКФ, так и результаты ADF-теста показывают, что временной ряд переменной X1 не является стационарным. Для получения стационарного временного ряда к нему была применена операция взятия первых разностей. График преобразованного временного ряда переменной X1, графики АКФ и ЧАКФ, а также результаты ADF-теста представлены на рис. 3, 4 и в табл. 3.

Таблица 2

Результаты расширенного теста Дики–Фуллера для переменной X1

Значение t-статистики ADF-теста

Уровни значимости

Критические значения

–2,43

1 %

–4,01

5 %

–3,44

10 %

–3,14

pic_78.tif

Рис. 3. Временной ряд первых разностей для переменной (коммерческий грузооборот транспорта) за период с января 1999 года по январь 2013 года

Autocorrelation

Partial Correlation

 

AC

PAC

Q-Stat

Prob

*|. |

*|. |

1

–0,154

–0,154

4,0414

0,044

.|. |

.|. |

2

0,030

0,007

4,1969

0,123

*|. |

*|. |

3

–0,071

–0,067

5,0647

0,167

.|. |

.|. |

4

0,070

0,051

5,9250

0,205

.|* |

.|* |

5

0,104

0,128

7,8089

0,167

*|. |

*|. |

6

–0,123

–0,099

10,473

0,106

*|. |

*|. |

7

–0,079

–0,114

11,584

0,115

*|. |

*|. |

8

–0,171

–0,198

16,750

0,033

.|* |

.|. |

9

0,149

0,073

20,710

0,014

.|. |

.|. |

10

–0,064

–0,030

21,440

0,018

*|. |

*|. |

11

–0,120

–0,132

24,035

0,013

.|. |

.|* |

12

0,061

0,077

24,718

0,016

Рис. 4. Автокорреляционная и частная автокорреляционная функции первых разностей временного ряда переменной X1 (фрагмент)

Так как ни один из коэффициентов автокорреляции у АКФ и ЧАКФ не является значимым, то можно предположить, что временной ряд первых разностей переменной X1 является стационарным. Для проверки данного предположения был проведен расширенный ADF-тест по первым разностям уровней этой переменной (табл. 3)

Таблица 3

Результаты расширенного теста Дики–Фуллера для первых разностей временного ряда переменной X1

Значение t-статистики ADF-теста

Уровни значимости

Критические значения

–14,99

1 %

–3,47

5 %

–2,88

10 %

–2,58

Таким образом, по результатам ADF-теста гипотеза о нестационарности временного ряда (наличии единичного корня) отвергается на всех рассмотренных уровнях значимости, следовательно, ряд первых разностей переменной X1 является стационарным, а исходный ряд – интегрированным первого порядка I(1).

Анализируя результаты исследования свойств временных рядов остальных экзогенных и эндогенных переменных, было получено, что они являются нестационарными с порядком интегрированности I(1). В свою очередь, первые разности данных временных рядов являются стационарными, то есть имеют порядок интегрированности I(0).

Результаты проведенного анализа показали, что временные ряды рассматриваемых показателей могут быть подвергнуты коинтеграционному анализу для выявления долговременных связей между ними. Далее в работе проводился коинтеграционный анализ пар временных рядов экономических показателей, у которых ранее (с помощью теста Грэнжера) было выявлено наличие причинно-следственной зависимости.

Исследуем временные ряды показателей X7 (денежные доходы в среднем на душу населения) и Y5 (начисленная среднемесячная номинальная заработная плата одного работника) на коинтеграцию. С этой целью было построено графическое изображение динамики изменения временных рядов данных показателей (рис. 5).

pic_79.tif

Рис. 5. Динамика денежных доходов в среднем на душу населения (X7) и начисленной среднемесячной номинальной заработной платы одного работника (Y5) за период с января 1999 года по январь 2013 года

Графический анализ динамики временных рядов показателей Y5 и X7 показал, что их тенденции совпадают. Чтобы исключить возможность ложной корреляции (согласно методу Энгла–Грэнжера), остатки модели парной линейной зависимости Y5 и X7 были исследованы на стационарность. С помощью метода наименьших квадратов (МНК) было получено следующее выборочное уравнение парной линейной зависимости:

Eqn58.wmf R2 = 0,995.

Далее выдвигалась основная гипотеза: временные ряды показателей Y5 и X7 не коинтегрированы. Для проверки основной гипотезы были найдены остатки регрессии et полученного уравнения, которые исследовались на стационарность с помощью расширенного теста Дики–Фуллера (табл. 4).

По результатам ADF–теста остатки регрессии et стационарны, следовательно, нулевая гипотеза об отсутствии коинтеграции между временными рядами показателей Y5 и X7 отвергается. Таким образом, временные ряды исследуемых показателей Y5 и X7 коинтегрированы.

Аналогично исследовались пары временных рядов каждой переменной правой части уравнений СОУ с каждой зависимой переменной Eqn59.wmf левой части уравнений СОУ. Были вычислены наблюдаемые значения t-статистики расширенного теста Дики–Фуллера для каждой пары переменных. Результаты расчетов представлены в табл. 5.

Таблица 4

Результаты расширенного теста Дики–Фуллера для остатков регрессии et

Значение t-статистики ADF-теста

Уровни значимости

Критические значения

–9,54

1 %

–2,58

5 %

–1,94

10 %

–1,62

Таблица 5

Наблюдаемые значения t-статистики ADF-теста

Зависимые переменные

Переменные правой части уравнений СОУ

Y1

Y2

–3,19

X3

–11,26

X4

–13,39

X6

–11,96

X8

–14,66

X9

–11,56

X11

–17,79

Y2

Y4

–10,63

Y6

–2,43

X1

–10,64

X3

–1,71

     

Y3

Y1

–16,09

Y6

–4,49

X1

–16,59

X4

–7,51

X10

–12,51

   

Y4

Y3

–3,19

X2

–3,19

X4

–3,19

X10

–3,19

     

Y5

Y1

–11,56

Y6

–5,46

X5

–7,79

X7

–9,54

X10

–12,13

   

Y6

Y1

–4,29

Y4

–4,51

X4

–4,08

X6

–5,16

X10

–5,19

   

На основе полученных результатов был сделан вывод о том, что рассмотренные пары временных рядов отобранных экономических показателей коинтегрированы на 10 %-м уровне значимости (критическое значение теста Дики-Фуллера при 10 %-м уровне значимости равно –1,62). Полученные выводы были использованы при построении эконометрической модели стабилизационных процессов в виде СОУ, в которой в качестве исходных данных использовались непосредственно уровни рассматриваемых рядов. Это позволило учесть информацию, содержащуюся в исходных данных, в полном объеме, не переходя, например, к первым разностям уровней. Полученный результат особенно важен при исследовании стабилизационных процессов, так как в этом случае нужно учитывать наличие необходимой для общества тенденции развития. Моделирование таких процессов без учета тенденции не позволит проанализировать долговременную зависимость между переменными и получить качественный прогноз.

Учитывая результаты коинтеграционного анализа, была получена структурная форма СОУ следующего вида:

Eqn60.wmf

где uti – случайные компоненты; αi, αik, γip – структурные параметры, которые были оценены по выборке (Eqn61.wmf Eqn62.wmf Eqn63.wmf Eqn64.wmf).

После оценивания структурных параметров предложенной модели двухшаговым методом наименьших квадратов (так как проверка условий идентифицируемости показала, что все уравнения СОУ сверхидентифицируемы) была получена следующая система уравнений:

Eqn65.wmf

Построенная эконометрическая модель удовлетворяет требованиям адекватности. Практически все оценки параметров уравнений регрессии статистически значимы на 5 %-м уровне значимости (в модели в скобках указаны расчетные значения t-статистики для соответствующих оценок параметров, здесь же указаны значения коэффициентов детерминации R2). Статистически значимы (по F-критерию) и все уравнения системы.

Полученная система одновременных уравнений может быть использована для анализа структуры взаимосвязей рассмотренных экономических показателей, а также для построения точечных и интервальных прогнозов [3].

Рецензенты:

Корнеева Т.А., д.э.н., профессор кафедры бухгалтерского учета и экономического анализа, ФГБОУ ВПО «Самарский государственный экономический университет», г. Самара;

Репин О.А., д.ф.-м.н., профессор, заведующий кафедрой математической статистики и эконометрики, ФГБОУ ВПО «Самарский государственный экономический университет», г. Самара.

Работа поступила в редакцию 11.07.2013.

Библиографическая ссылка

Суханова Е.И., Ширнаева С.Ю. ИСПОЛЬЗОВАНИЕ КОИНТЕГРАЦИОННОГО АНАЛИЗА ПРИ ПОСТРОЕНИИ CИСТЕМЫ ОДНОВРЕМЕННЫХ УРАВНЕНИЙ // Фундаментальные исследования. – 2013. – № 8-5. – С. 1158-1164;
URL: https://fundamental-research.ru/ru/article/view?id=32102 (дата обращения: 23.04.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674