А. В. Полбин
Эконометрическая оценка структурной макроэкономической модели российской экономики
В статье анализируется применимость динамических стохастических моделей общего равновесия с широким набором номинальных и реальных «жесткостей» для изучения делового цикла российской экономики. Оценивается структурная модель российской экономики и производится разложение циклической компоненты основных макроэкономических переменных по экзогенным шокам.
ключевые слова: динамические стохастические модели общего равновесия; малая открытая экономика; деловые циклы; байесовский подход в эконометрике. JEL classification: C11; E32; E40; E47; F41.
1. введение
В последнее десятилетие был достигнут значительный прогресс как в спецификации, так и в оценивании динамических стохастических моделей общего равновесия (DSGE). Целью построения DSGE моделей, как правило, является изучение делового цикла. Но, в отличие от векторных авторегрессий, DSGE модели прочно опираются на экономическую теорию и предлагают формальный экономико-математический аппарат для изучения факторов делового цикла и анализа экономической политики, с данной точки зрения они меньше подвержены критике Лукаса (Lucas, 1976). С другой стороны, поскольку DSGE модели прочно основаны на теории, они часто являются слишком стилизованными, чтобы можно было привести их в непосредственное соответствие с данными (см., например, (Ireland, 2004)).
В качестве ключевых работ относительно спецификации и оценки DSGE моделей можно выделить (Christiano et al., 2001, 2005; Smets, Wouters, 2003, 2007). В работах (Christiano et al., 2001, 2005) было показано, что DSGE модель с широким набором реальных и номинальных жесткостей, таких как жесткость цен и номинальных заработных плат1, привычки в потреблении, издержки на установку капитала и издержки интенсивности загрузки капитальных мощностей, достаточно хорошо способна воспроизводить динамические отклики на шоки денежно-кредитной политики.
В работах (Smets, Wouters, 2003, 2007) авторы строят DSGE модели для Европы и США, спецификации которых в целом согласуются с работами (Christiano et al., 2001, 2005), и оценивают их с помощью байесовских эконометрических методов. Согласно результатам (Smets, Wouters, 2003, 2007), оцененные DSGE модели обеспечивают хорошую согласованность
1 Это означает, что данные ценовые показатели не являются гибкими в краткосрочном периоде и приспосабливаются к изменениям спроса только с течением некоторого времени.
моделей с эмпирическими данными и конкурируют в прогнозных свойствах с сугубо эко-нометрическими моделями.
В настоящей статье делается попытка применить данный подход для анализа делового цикла российской экономики. Несомненно, при этом возникает множество ограничений. Достаточно спорным является вопрос о применимости DSGE моделей для изучения делового цикла в развивающихся экономиках. В работе (Aguiar, Gopinath, 2007) авторы приходят к выводу, что основным источником флуктуаций развивающихся экономик являются шоки трендового роста, а не временные колебания вокруг детерминированного тренда.
Вопрос о применимости DSGE моделей является важным и для российской экономики. В частности, возможно, больший интерес заключается в изучении именно трендовой, а не циклической компоненты динамики отечественной экономики. Здесь также возникает проблема разделения динамики макроэкономических переменных на трендовую и циклическую составляющие. Вызывает сомнение вопрос о том, насколько хорошо технические методы фильтрации могут произвести такую декомпозицию для имеющихся данных.
С другой стороны, источники флуктуаций развивающихся экономик могут быть совершенно другими по сравнению с развитыми экономиками. Например, в работе (Mendoza, 1995) приводятся доводы в пользу условий торговли в качестве превалирующей движущей силы делового цикла в развивающихся экономиках. В работе (Uribe, Yue, 2006) — в пользу шоков премии за риск к отечественным активам в качестве важного источника делового цикла, при этом наличие процентного спрэда к отечественным активам, величина которого зависит от внутренних факторов, может значительно усиливать деловой цикл в развивающихся экономиках.
Предлагаемая в настоящей работе модель является достаточно стилизованной и не призвана полностью объяснить действительность. В работах (Smets, Wouters, 2003, 2007) рассматривался обширный набор структурных шоков экономики2, и соответствующие экзогенные переменные являлись ненаблюдаемыми (латентными). Можно утверждать, что провести аналогичный качественный эконометрический анализ для российской экономики невозможно из-за короткой длины временных рядов, что непосредственно обуславливает проблему идентификации структурных шоков на имеющихся данных.
В настоящей работе ставится лишь вопрос о том, насколько ряд экзогенных шоков способен объяснить деловой цикл отечественной экономики. В частности, эконометрические оценки влияния мировых цен на нефть на российскую экономику (см., например, Казакова, Синельников-Мурылев, 2009; Мельников, 2010; Ito, 2009) говорят о достаточно сильной зависимости отечественной экономики от данного фактора внешнеэкономических условий. Особый интерес представляет вопрос, насколько изменения цен на нефть способны объяснить деловой цикл в российской экономике в рамках теоретической модели общего равновесия.
Макроэкономические эффекты влияния на российскую экономику изменений цен на нефть в контексте DSGE моделей изучались в рамках имитационного анализа, например, в работах (Sosunov, Zamulin, 2006а, 2006б). В настоящей статье рассматривается модель с большим набором механизмов, обеспечивающих динамику макроэкономических переменных, для ее оценки используются байесовские эконометрические методы.
2 В работе (Smets, Wouters, 2003) рассматривалось десять структурных шоков и семь наблюдаемых переменных, а в (Smets, Wouters, 2007) — семь шоков и семь переменных.
В качестве структурных шоков рассматриваются: 1) шок экспортных доходов, опреде- § ляющим фактором которых являются мировые цены на нефть; 2) шок обменного курса; § 3) шок совокупной факторной производительности (СФП) и 4) шок премии за риск к отече- ^ ственным активам. При этом в эконометрическом анализе экспортные доходы и обменный курс используются в качестве наблюдаемых переменных, что в некотором смысле упрощает проблему идентификации. Шок премии за риск к национальным активам трактуется как некоторый обобщенный фактор несовершенства финансового рынка и шоков финансового сектора.
Работа построена следующим образом. В разделе 2 приводится описание теоретической модели, в разделе 3 описываются используемые данные и результаты эконометрического анализа. В заключении приводятся основные выводы и результаты работы.
2. теоретическая модель
В данном разделе приводится описание теоретической модели российской экономики. Предлагаемая модель, основанная на работах (Christiano et а1., 2001, 2005; Smets, Wouters, 2003, 2007), представляет собой малую открытую экономику с пятью отечественными экономическими агентами: 1) домохозяйства; 2) центральный банк; 3) фирмы, производящие товар конечного потребления; 4) фирмы, производящие отечественный промежуточный товар; 5) импортирующие фирмы.
Разделение всех фирм на три типа позволяет моделировать жесткость цен как на отечественные промежуточные товары, так и на импортные товары, а также анализировать эффект неполного переноса изменений валютного курса в цены импортных товаров на отечественном рынке. Фирмы, производящие товары конечного потребления, используют в качестве факторов производства промежуточные импортные и отечественные товары и производят продукт, который в дальнейшем идет на потребление домашних хозяйств и инвестиции. Такое предположение является сильным упрощением и не учитывает различное соотношение импортных и отечественных товаров в инвестициях и потреблении.
На данном этапе построения модели происходит абстрагирование от фискального блока и отнесение фискального сектора к экономической деятельности домохозяйств. В дальнейшем предполагается расширение модели до детализированного описания экономической деятельности фискального блока.
Совокупный экспорт моделируется как экзогенная переменная, представляющая собой чистый трансферт из-за границы. Такая предпосылка может быть мотивирована преобладанием нефти, газа и других сырьевых товаров в экспорте российской экономики.
2.1. Домохозяйства
Предполагается, что в экономике существует континуум домохозяйств т е [0,1]. Каждое домохозяйство максимизирует свое ожидаемое благосостояние иг (т), являющееся дисконтированной суммой полезностей от потребления Сг (т) и количества отработанных часов Ц (т). Полезность домохозяйства положительно зависит от потребления и отрицательно —
от количества отработанных часов. Предполагается следующая сепарабельная функция благосостояния:
где Ег — оператор условного математического ожидания на основе всей доступной информации на момент времени I, ос — обратная величина к межвременной эластичности замещения, оь — обратная величина к эластичности предложения труда по реальной заработной плате, т — нормировочная константа, 3 — субъективный коэффициент дисконтирования, Н — переменная, характеризующая привычки домохозяйств в потреблении, отражающие инерционность динамики потребления домохозяйств во времени.
В соответствии с (Smets, Wouters, 2003, 2007) считается, что переменная привычек домохозяйств в потреблении Нг пропорциональна агрегированному потреблению в прошлом периоде времени:
В модели предполагается, что домохозяйства владеют капиталом и сдают его в аренду фирмам. При этом объем капитала, который арендуется фирмами, может быть изменен как за счет соответствующих инвестиционных решений, так и за счет интенсивности загрузки физического капитала. В свою очередь, интенсивность загрузки физического капитала связана с некоторыми реальными издержками в терминах конечного выпуска.
В модель также вводится жесткость номинальных заработных плат. Предполагается, что домохозяйства наделены дифференцированным трудом и имеют некоторую монопольную власть на рынке труда. Таким образом, каждое домохозяйство устанавливает цену на свой труд при эластичной функции спроса на труд.
При моделировании жесткости зарплат используется подход к ценообразованию по (Calvo, 1983) и предполагается, что домохозяйство в каждый период времени может изменить свою заработную плату с некоторой экзогенной вероятностью (1 — £W). Привлекательность ценообразования по (Calvo, 1983) заключается в том, что в данном подходе естественным образом возникает понятие продолжительности ценового контракта и параметр XW имеет наглядную экономическую интерпретацию. Но в случае моделирования жесткости заработных плат подход (Calvo, 1983) сильно усложняет задачу, т. к. среди домохозяйств возникает гетерогенность. Они начинают сталкиваться как с агрегированной неопределенностью, так и с индивидуальной неопределенностью ввиду вариации в индивидуальных трудовых доходах.
Чтобы обойти эту проблему, в модель вводятся отечественные контингентные ценные бумаги, позволяющие застраховать домохозяйство от индивидуальной неопределенности. Такая предпосылка является стандартной при построении современных DSGE моделей (см., например, (Erceg et al., 2000; Christiano et al., 2001, 2005; Smets, Wouters, 2003, 2007)). Она обеспечивает, что в равновесии домохозяйства являются однородными относительно потребления и владения активами. Неоднородность наблюдается только для отработанных часов и заработной платы.
Домохозяйства максимизируют свое благосостояние при следующих динамических бюджетных ограничениях:
(1)
H. = hC,
Bt (r) StB't (r) Wt(т) ^ Rk v , ч I
4- t t V ' --tV ' T t\-I--V Ъ
Ct (т) + Л (т) + ^ + ^^ = (т) + (т) +
Ri = Rt exP
Чв; ^
\ -'t-'t /
+ e R
где Г* — агрегированный конечный выпуск, который будет определен ниже.
Спецификация данной зависимости соответствует работе (Adolfson et а1., 2007). Параметр рв определяет чувствительность эндогенной составляющей премии за риск к изменениям в сбережениях (долге) национальной экономики. В модель также вводится экзогенная составляющая премии за риск е^, которая описывается следующим стационарным AR(1) процессом:
е Я =РК е ^ + П^,
где рЯ — коэффициент автокорреляции, и~ N(0, оЯ) — шок премии за риск к национальным активам.
Вторым ограничением оптимизационной задачи домохозяйств является уравнение динамики накопления физического капитала:
р р р 1 р 1 1-1 б
+Щ) + зДф^ + рр1+^+АТ _ ^ ткт I
р р р р
где Сг (т) — потребление, (т) — инвестиции, Вг (т) — номинальная стоимость отечественных облигаций, В*(т) — номинальная стоимость облигаций в иностранной валюте, приобретенных на внешнем финансовом рынке, р — индекс агрегированных цен товаров конечного потребления, 5* — номинальный обменный курс, Щ (т) — номинальная заработная плата, Я — рентная цена капитала, и — интенсивность загрузки капитальных мощностей, К (т) — физический капитал, Я — валовая доходность по отечественным облигациям, ЯЯ* — валовая доходность по вложениям в иностранные облигации, Ртг (т) — прибыль фирм, полученная в виде дивидендов, X* — экспортный доход, А (т) — отечественные контингентные ценные бумаги, G (и* (т)) — издержки на загрузку единицы физического капитала с интенсивностью иг (т).
Экспорт рассматривается как экзогенный трансферт дохода из-за границы, а динамика экспорта моделируется в виде следующего стохастического процесса:
1п(Х,) = (1 - рх) 1п(X) + Рх 1п(X) + иХ,
где рх — коэффициент автокорреляции, X — долгосрочный уровень экспорта, ^ ~ N(0, оX) — шок экспортного дохода.
В модели также предполагается, что валовая доходность Я* по вложениям в иностранные облигации определяется с некоторой премией за риск (к безрисковой доходности) Я* на внешнем финансовом рынке, зависящей от агрегированного уровня сбережений (долга) отечественной экономики в долях ВВП. Предполагается следующая функциональная форма:
Kt (t) = (1-d)K- (t) +
1-Ф
{ т t \ W
I (t)
It-1 (t)
\ t-
It (t),
))
где д — норма выбытия капитала, Ф() — издержки на установку нового капитала. Предполагается, что Ф (1) = 0, Ф'(1) = 0 и Ф" (1) = р.
В соответствии с работами (Christiano et а1., 2001, 2005) предполагается, что издержки на установку нового капитала зависят от темпа изменения инвестиций. Такая спецификация обеспечивает инерционную динамику инвестиций в ответ на шоки экономики.
Определим также переменную загруженного капитала К( (т), который арендуется фирмами у домохозяйств и в дальнейшем используется в производственном процессе:
К; (т) = и (т)К— (т).
Как отмечалось выше, предпосылка о существовании отечественных контингентных ценных бумаг обеспечивает, что в равновесии домохозяйства являются однородными относительно потребления и владения активами. Условия оптимальности задачи домохозяйства относительно выбора потребления, объема покупки отечественных облигаций и внешних облигаций приводит к следующим уравнениям (индекс т опущен):
At ={Ct - НС— )С,
Л t =bEt
Л.
R
Л t =bEt
Л
R St+i pt+i St
где Л, — множитель Лагранжа при бюджетном ограничении, являющийся теневой ценой потребления.
Условия оптимальности по физическому капиталу, инвестициям и интенсивности загрузки капитала приводят к следующим динамическим уравнениям (индекс т опущен):
Q = Et
b^(Qt+1 (1 — d) + uMrth -G(ut+i))
\2
Qt
4 x
\It~1) it-1
L t t+1 A
\ t /
+ 1 = Qt
\ t /
1-Ф
Ч,л
It-1 ))
гк = С ' (и,),
где Qt — теневая реальная стоимость единицы капитала, г/ = Я /р — реальная рентная цена.
Предполагается, что С(1) = 0 и С'(1) = 1/в -1 + 5. Введем параметр 'фи = С'(1)/С"(1), представляющий собой эластичность интенсивности загрузки физического капитала по рентной цене.
Перейдем к задаче выбора номинальной заработной платы. Предполагается следующий вид спроса на труд домохозяйства г:
L (т) =
г W (т) v(1+A)/Ar
W
v t
L
L =
W =
^ (т))
[£(w (т) У
i/(i+Ar)
ёт
1/Aw
ёт
1 + AW
-Ar
Wt (t) =
P_
P
\ t—2 /
W-i (t):
W (t)
P
Et ж v
(P / p—i )n
p / p
t+s ' t+s—l
L+ (t)V+ (t)
1 + A„
+ Et +s (t) L+s (t) = 0,
J I
CQ «i
где (1+ ) / — эластичность спроса на труд отдельного домохозяйства, Ц — агрегированный спрос на труд и — агрегированный уровень зарплаты, которые определяются с помощью следующих функций агрегирования:
Предполагаем, что домохозяйство в каждый период времени с экзогенной заданной вероятностью (1 ) получает сигнал, при реализации которого оно устанавливает зарплату исходя из максимизации своего благосостояния (1).
В случае неполучения данного сигнала в настоящей модели, как и в работах (Smets, Wouters, 2003, 2007), предполагается частичная индексация номинальных заработных плат к инфляции агрегированного уровня цен:
где уш — степень индексации заработных плат. Если уш = 0, то не происходит индексации номинальных заработных плат, а если уш = 1, то номинальные заработные платы полностью индексируются к инфляции прошлого периода.
Условие оптимальности для домохозяйств, максимизирующих в данный период времени зарплату Щ (г), принимает вид:
где Vt (т) =
(C (т) - Ht у
1 1 + aL
(Lt (т)У L, VtC (т) = (Ct (т) - Ht) стс — предельная полез-
ность потребления, УЩ (г) = — (Ц (г)) Ц — предельная «пагубность» работы. Индекс г при потреблении можно опустить ввиду гомогенности по данной переменной.
Поскольку все домохозяйства, оптимизирующие зарплаты в данный период времени, устанавливают одну и ту же зарплату, общий индекс номинальных зарплат определяется в следующем виде: , \-1/ящ
(W )—1/v =Xr
W_
P_
P
Vt— 2 /
+ (1 —Xw )(Wt (t))
2.2. Производство товара конечного потребления
В модели предполагается, что фирмы, производящие товар конечного потребления, действуют на рынке совершенной конкуренции и в качестве факторов производства используют дифференцированные промежуточные отечественные товары и импортные товары. Производственный процесс в данном секторе описывается с помощью следующей функции Коб-ба -Дугласа: /„л \1-®
_(Im£M
' тт (1-т)т '
где — агрегированный конечный выпуск, 1т{ — агрегированный импорт, У^ — агрегированный выпуск отечественных промежуточных товаров, т — эластичность агрегированного конечного выпуска по импорту. В работе (Карев, 2009) также предполагалась функция Кобба -Дугласа при формализации выбора между отечественными и импортными товарами.
В свою очередь, 1т{ и У^ определяются с помощью следующих функций агрегирования:
Imt =
Yd =
1
io (m (0)
1/(1+Xpm)
Г1 / , 41/(1+Xpd) i0( j d
(2)
(3)
где 1т{ (7) — товар 7-й фирмы-импортера, Г/ (/) — товар]-й фирмы, производящей отечественные промежуточные товары. В уравнениях (2) и (3) интеграл берется по континууму фирм-импортеров и производителей отечественных промежуточных товаров соответственно. Параметры Xрт и Xра определяют величину маржи в соответствующих секторах.
Задача минимизации расходов фирмы, производящей конечный продукт, приводит к следующим функциям спроса на отдельные промежуточные товары:
\-(1 + ^рт )¡^рг.
Imt (i) =
Pm (i)
pm
V t
Ytd (j) =
/ pd (j)\-(m pd)'я pd pd
p
a—^Y
m
p
(1-w) p^Y
где Р. (7) — цена товара 7-й фирмы-импортера, Р/ (j) — цена товара]-й фирмы, производящей отечественные промежуточные товары, Рт — агрегированный уровень цен импортных товаров на внутреннем рынке, р — агрегированный уровень цен отечественных промежуточных товаров, Р — агрегированный уровень цен товаров конечного потребления.
В свою очередь, индексы цен определяются следующими выражениями:
P"
pd =
=[/0(pm (i))-[10 (Pd (j)
p=(pm ) (pd)
X
2.3. Импортеры Ц
\о
I
Щ
В настоящей модели предполагается, что существует отдельный сектор импортирующих фирм, которые покупают однородный импортный товар на внешнем рынке, дифференцируют его, после чего продают на внутреннем рынке. Считается, что импортирующие фирмы не несут каких-либо дополнительных издержек на производство своего промежуточного продукта, кроме затрат на покупку данного товара на внешнем рынке.
При моделировании ценообразования на импортные товары на внутреннем рынке используется подход (Calvo, 1983) и предполагается, что импортирующие фирмы в каждый период времени с экзогенной заданной вероятностью (l получают сигнал, при реализации которого они устанавливают цену исходя из максимизации своей целевой функции. В случае неполучения данного сигнала в настоящей модели, как и в работах (Smets, Wouters, 2003, 2007), предполагается частичная индексация цен к инфляции агрегированного уровня цен на импортные товары в прошлом периоде. Такая предпосылка обеспечивает включение первого лага инфляции в кривой Филлипса и вводится для учета инерционности цен.
В случае когда фирмы получили соответствующий сигнал о возможности оптимизировать цену, они устанавливают цену на свой товар, принимая во внимание, что в будущем они смогут изменить (оптимизировать) свою цену только с определенной вероятностью. Целевая функция оптимизирующей цену фирмы принимает следующий вид:
да Л г
e xP^ [(((° s - ^ P s) s (i).
t t + s
Выражение в квадратных скобках представляет собой прибыль фирмы в момент времени t + s, величина St+sPt+* является предельными издержками фирм-импортеров, которые соответствуют затратам на приобретение импортного товара на внешнем рынке и равны произведению номинального обменного курса на цену импортного товара в иностранной валюте. Здесь также предполагается, что фирмы действуют в интересах домохозяйств, и,
соответственно, выражение —представляет собой стохастический дисконт-фактор,
Л' р+*
отражающий текущую ценность для домохозяйства дополнительного рубля в момент времени t + s. Фирмы принимают решение по цене рт (/), которая будет установлена в текущем периоде, с учетом того, что они не смогут оптимизировать свою цену в момент времени t + s с вероятностью [Х,тр ) . Если фирмы не оптимизируют цену в момент времени t + s, то они индексируют ее на прошлую инфляцию цен импортных товаров на внутреннем рынке со степенью индексации утр .
Условие оптимальности задачи импортирующей фирмы, оптимизирующей свою цену в данный момент времени, приводит к следующему уравнению:
да Л
Et(i) "s
Л
Pm (i) ( pm ) ni,
--(1 + Km )-
pm p / pm v Vm> p
t t + s ' t t + s
= 0.
Поскольку все фирмы, оптимизирующие цены в данный период времени, устанавливают одну и ту же цену, общий индекс импортных цен определяется в следующем виде:
s=0
(Ptm)
-1/я „,
_ hn.
x p
pn
{ pm Vp
pm
\ t—2 /
\-1/я p
+
(i-xm )(pppm ( j))~
2.4. Производство отечественного промежуточного продукта
Предполагается, что множеству отечественных товаров соответствует единичный континуум фирм, действующих на рынке монополистической конкуренции. Технология каждой фирмы описывается с помощью производственной функции Кобба -Дугласа:
Yd (j)=< (( (j)(j))
где г" — экзогенный уровень СФП, К( (j) — количество загруженного капитала, арендуемого j-й фирмой у домохозяйств, Ц (j) — количество единиц использованного дифференцируемого труда, а — эластичность выпуска по капиталу.
Предполагается, что экзогенный уровень СФП описывается следующим стохастическим процессом:
^) = Ра ^«-1)+<, где ра — коэффициент автокорреляции, и" ~ N(0, о2а) — шок СФП.
Задача минимизации издержек фирмы на производство единицы продукции приводит к следующему оптимальному соотношению между трудом и загруженным капиталом:
КЛЛ Щ а
L (j) Rf 1-a
Таким образом, оптимальное соотношение труда и капитала определяется отношением цен данных факторов, которые для всех фирм одинаковы. Следовательно, фирмы арендуют труд и капитал в одинаковых пропорциях вне зависимости от объема производства, и предельные издержки определяются следующим выражением:
1 (W )1-a(f г
MCt (j) = —v ' 1 ' . tJ < (1-a) aa
При моделировании ценообразования на отечественные промежуточные товары используется подход (Calvo, 1983) и предполагается, что отечественные фирмы в каждый период
времени с экзогенной заданной вероятностью (1- £dp ) получают сигнал, при реализации
которого фирмы устанавливают цену исходя из максимизации ожидаемой суммы дисконтированных прибылей. При этом, установив цену на товар, фирмы обязуются удовлетворить весь спрос на свою продукцию. Так как фирмы в секторе производства отечественных торгуемых товаров сталкиваются с одинаковыми предельными издержками, их оптимизацион-
ная задача по выбору цены полностью аналогична задаче импортирующих фирм. Целевая функция записывается в следующем виде:
E,
2(Р b)' Лр" [(( ( jK+s - MC1+s ) (j) ].
\Pt+s
J I
CQ «i
Если фирмы не оптимизируют цену в момент времени t + s, то они индексируют ее на прошлую инфляцию отечественных цен на промежуточные товары со степенью индексации ур . Условие оптимальности задачи отечественной фирмы, оптимизирующей свою цену в данный момент времени, приводит к следующему уравнению:
E b)'Yi' (j)
л,
p d (j) (p^/ P- )g p
p / pd
-1 -1 i
- (1 + A pd)
MCh p
= о.
Агрегированный уровень цен на отечественные промежуточные товары записывается
в виде:
(Pd )-1/*pd =xd
Pd
/ pd V~p Pt-1 pd V t-2 /
\-1/* p.
+
(1-Xp )((d (j))
-1/* p.
X
s=0
'=0
2.5. Центральный банк
Предполагается, что центральный банк следует политике управляемого номинального обменного курса. Для описания в модели деятельности центрального банка предполагаем, что номинальный обменный курс является экзогенной переменной, динамика которой описывается следующим стохастическим процессом:
Р - -р (1 + и ),
где р — темп роста номинального обменного курса, Р — некоторый долгосрочный темп роста обменного курса, и г ~ N(0, о5) — шок обменного курса.
Данное предположение является существенным упрощением ввиду сложности описания политики Центрального банка РФ с помощью инструментальных правил политики, учитывающих эндогенную реакцию денежно-кредитных властей на ситуацию в экономике. Данная спецификация позволяет описывать динамику макроэкономических переменных, но возникают трудности при интерпретации шоков денежно-кредитной политики (ДКП). Например, девальвации во время кризисов 1998 и 2008 гг. были в некотором смысле эндогенными реакциями центрального банка на экономическую ситуацию, а не экзогенными шоками политики.
2.6. Равновесие на рынках
Рынок товаров конечного потребления находится в равновесии, если количество произведенного товара конечного потребления равно спросу со стороны домохозяйств на потребление и инвестиции, а также на издержки, связанные с загрузкой капитала:
Yt = Ct + It + G(ut (t ))K- (t).
Рынок отечественных промежуточных товаров находится в равновесии, если объем производства отечественных промежуточных товаров равен спросу на данные товары. Рынок капитала находится в равновесии, если спрос на загруженный капитал со стороны фирм равен предложению загруженного капитала со стороны домохозяйств. Рынок труда находится в равновесии, если спрос на труд со стороны фирм равен предложению труда со стороны домохозяйств при ставках заработных плат, устанавливаемых домохозяйствами. А равновесие платежного баланса формализуется следующим образом:
(B* - B-1) = X - p* Imt + (R-1 -1) B-1.
Все номинальные показатели настоящей модели были продефлированы на агрегированный уровень цен товаров конечного потребления Pt.
3. Результаты эконометрического анализа
В настоящей статье будем следовать подходу байесовской эконометрики, который использовался при оценивании DSGE моделей, например, в работах (Adolfson et al., 2007; Del Negro, Schorfheide, 2008; Smets, Wouters, 2003, 2007). Динамическая система нелинейных уравнений модели приводилась к линейным разностным уравнениям с рациональными ожиданиями с помощью лог-линеаризации около долгосрочного равновесия. Для решения линейных разностных уравнений с рациональными ожиданиями использовался алгоритм, предложенный в работе (Sims, 2002). Для построения функции правдоподобия использовался фильтр Калмана. Эконометрическое оценивание проводилось в специализированном программном пакете Dynare (Adjemian et al., 2011).
3.1. Описание используемых данных
При оценке модели использовались данные для России с 1-го квартала 1995 г. по 4-й квартал 2012 г. (72 наблюдения) по 8-ми макроэкономическим переменным:
1) реальный выпуск отечественных промежуточных товаров;
2) реальное потребление;
3) реальные инвестиции;
4) реальная заработная плата;
5) реальный импорт;
6) относительные цены импортных товаров к внутренним товарам;
7) темп роста номинального эффективного курса рубля;
8) долларовый экспорт.
Первые шесть переменных являются эндогенными переменными, последние две — экзогенными. Для построения первых шести переменных использовались данные Росстата, для номинального эффективного курса рубля использовались данные МВФ (IMF, International Financial Statistics), для экспорта в долларовом выражении использовалась статистика Банка России.
В последние годы Росстат перешел от публикации статистики национальных счетов в по- § стоянных ценах 2003 года к статистике в ценах 2008 года. Для построения согласованных § временных рядов реальных показателей рассматриваемые в настоящей работе ряды с 1-го ^ квартала 1995 г. по 4-й квартал 2003 г. (в постоянных ценах 2003 года) продлевались на период с 1-го квартала 2004 г. по 4-й квартал 2012 г. с использованием индексов физических объемов рассматриваемых реальных показателей в ценах 2008 года в процентах к соответствующему кварталу предыдущего года. Тем самым для «сцепления» рядов была выбрана середина периода.
В качестве реального потребления рассматривалась сумма расходов на конечное потребление (в постоянных ценах 2003 г.) домашних хозяйств, государственного управления и некоммерческих организаций, обслуживающих домашние хозяйства. В качестве реальных инвестиций рассматривалась сумма валового накопления основного капитала и изменения запасов материальных оборотных средств в постоянных ценах 2003 года. В качестве реального импорта рассматривался импорт в постоянных ценах 2003 года. Реальный выпуск отечественных промежуточных товаров определялся как сумма реального потребления и реальных инвестиций за вычетом реального импорта.
Совершенно аналогично можно получить временные ряды рассматриваемых переменных в текущих ценах и построить соответствующие дефляторы цен. Таким образом, в качестве относительных цен импортных товаров к внутренним товарам использовалось отношение дефлятора цен импорта к дефлятору цен выпуска отечественных промежуточных товаров. В качестве реальной заработной платы бралось отношение средней номинальной заработной платы к дефлятору потребления. Согласно методологии статистических данных, средняя номинальная заработная плата определяется как заработная плата на одного работника и, таким образом, может являться не очень хорошим «эталоном» для соответствующей модельной переменной, которая определяется как ставка зарплаты за отработанный час. В качестве экспорта использовался совокупный экспорт Российской Федерации в долларовом выражении.
Из всех рядов удалена сезонная составляющая с помощью Х-12-АММА. Далее с помощью фильтра Ходрика-Прескотта была выделена трендовая составляющая. Вычисления проведены в программе EViews. Предварительная фильтрация данных является стандартной при работе с DSGE моделями, но, как отмечалось во Введении, использование данного подхода может являться весьма спорным при оценке на российских данных. В качестве циклической компоненты рассматривалось отклонение в долях сезонно сглаженного ряда от тренда. Полученные циклические компоненты в дальнейшем трактовались как отклонения от долгосрочного равновесия.
3.2. Оценка параметров модели
В теоретической модели предполагается четыре структурных шока и, соответственно, четыре экзогенные переменные. При этом две экзогенные переменные в эконометрическом анализе используются в качестве наблюдаемых. Чтобы решить проблему сингулярности, в модель вводятся шесть ошибок измерения для всех наблюдаемых эндогенных переменных. Предполагается, что ошибки измерения имеют нормальное распределение, не корре-лированы во времени, между собой и со структурными шоками.
Предварительно перед проведением эконометрического анализа был зафиксирован ряд параметров, отвечающих за долгосрочное равновесие. Оценка же проводилась для параметров, обуславливающих динамику экономики.
Предполагается, что Ь = 0.99. Это соответствует реальной процентной ставке, равной 1% в квартал. Норма амортизации б калибруется на уровне 2.5% в квартал. Выбор таких значений является стандартным в литературе (см., например, (Prescott, 1986; Smets, Wouters, 2003)). Доля импорта т в конечном выпуске бралась на уровне 0.26, что соответствует средней доле импорта в потреблении и инвестициях в текущих ценах за рассматриваемый период времени. Доля капитала в производственной функции — 0.4. Параметр , как и в (Smets, Wouters, 2003), полагался на уровне 0.5, а Хра калибровался так, чтобы обеспечить долю потребления в конечном выпуске в долгосрочном равновесии на уровне 0.78, что является средней долей на рассматриваемом промежутке времени. Параметр Хрт калибровался на таком же уровне.
Остальные параметры оценивались с помощью метода Байеса. Спецификация априорных распределений параметров проводилась, в основном, согласно работам (Smets, Wouters, 2003, 2007). Для коэффициентов автокорреляции рх, ра и ря экзогенных переменных предполагается бета-распределение со средним 0.85 и стандартным отклонением 0.1 (Smets, Wouters, 2003). Бета-распределение ограничивает диапазон возможных значений коэффициента автокорреляции на интервале от 0 до 1. При этом заданные значения параметров априорного распределения закладывают высокую долю вероятности того, что рассматриваемые экзогенные переменные обладают высокой инерционностью.
Для стандартных отклонений структурных шоков и ошибок измерения предполагается обратное гамма-распределение с двумя степенями свободы. При этом для структурных шо-ков априорное среднее равно 0.05, а для ошибок измерения — 0.025.
Спецификация априорных распределений параметров предпочтений домохозяйства также выполнялась на основе (Smets, Wouters, 2003), при этом предполагалось, что обратная величина межвременной эластичности замещения ос имеет априорное нормальное распределение со средним 1 (что соответствует логарифмической спецификации предпочтений относительно потребления) и стандартным отклонением 0.375. Обратная величина к эластичности предложения труда по заработной плате оЦ имеет априорное нормальное распределение со средним 2 и стандартным отклонением 0.75. Параметр привычек потребления ограничен на интервале от 0 до 1, и для него задается априорное бета-распределение со средним 0.7 и стандартным отклонением 0.1.
При оценке параметра эластичности интенсивности загрузки физического капитала по рентной цене ри использовался подход, изложенный в (Smets, Wouters, 2007), и проводилась следующая нормализация:
р = ^.
РР и
Оценка будет проводиться для параметра Рр и, который предполагается ограниченным на интервале от 0 до 1. Если Рри стремится к нулю, то эластичность интенсивности загрузки физического капитала по рентной цене стремится к бесконечности. Если Рр и стремится к единице, то эластичность стремится к нулю. Для параметра Рри задается априорное бета-распределение со средним 0.5 и стандартным отклонением 0.15 (Smets, Wouters, 2007).
Для параметра издержек на установку капитала задается априорное нормальное распре- § деление со средним 4 и стандартным отклонением 1.5 (Smets, Wouters, 2007). При специфи- § кации параметров ценообразования задается априорное бета-распределение со средним 0.5 ^ и стандартным отклонением 0.15 для параметров индексации цен и зарплат, параметров вероятности изменить цены и зарплаты (Smets, Wouters, 2007). При данной спецификации априорного среднего средняя продолжительность ценового контракта составляет два квартала. В качестве априорного распределения для чувствительности эндогенной составляющей премии за риск по уровню сбережений (долга) национальной экономики задается обратное гамма-распределение со средним 0.2 и с двумя степенями свободы, что согласуется с эконо-метрическими оценками DSGE модели для Европы из работы (Adolfson et а1., 2007).
Оценка апостериорного распределения проводилась с помощью методов Монте-Карло по схеме марковских цепей (МСМС). При оценке использовались две цепи с длиной итераций в каждой 1 млн. Для анализа качества сходимости марковских цепей к эргодическо-му распределению использовались тесты, предложенные в работе (Brooks, Ge1man, 1998), которые показали достаточно хорошую сходимость алгоритма. В таблице 1 в столбцах под заголовком «Базовая спецификация» представлены моды, апостериорные средние и 90%-ные доверительные интервалы для каждого оцениваемого структурного параметра (оценки для стандартных ошибок измерения не приведены).
Опишем результаты оценивания ряда параметров.
Как следует из проведенного эконометрического оценивания, апостериорные распределения параметров ос и аь не сильно отличаются от заданных априорных распределений, что свидетельствует о слабой идентификации данных параметров на имеющихся данных.
Распределение параметра к привычек в потреблении демонстрирует сильный сдвиг относительно априорного распределения в сторону более высоких значений. Апостериорное среднее к равно 0.94, в то время как оценки для Европы (Smets, Wouters, 2003) и США (Smets, Wouters, 2007) составляли 0.59 и 0.71 соответственно. Таким образом, результаты оценивания специфицированной модели говорят о том, что в российской экономике более высокие привычки в потреблении и динамика потребления обладают большей инерционностью. Данные выводы могут быть ограничены, например, тем, что в модели не разделяется потребление государства и домохозяйств.
Апостериорное среднее для параметра издержек на установку капитала < равно 4.65, апостериорное среднее для ри равно 0.825, что соответствует эластичности интенсивности загрузки капитала по рентной цене, равной 0.21. Величина эластичности интенсивности загрузки капитала по рентной цене достаточно низка.
Также результаты оценивания свидетельствуют в пользу более высокой жесткости номинальных зарплат по сравнению жесткостью цен. Средняя продолжительность контракта номинальной заработной платы 1/ (1 - ^) в точке апостериорного среднего составляет 3.6 квартала. Средняя продолжительность ценового контракта для импортных цен и цен на отечественные промежуточные товары составляет 3.1 и 1.4 квартала соответственно. Таким образом, у импортирующих фирм меньше вероятность получения сигнала, при котором они могут изменить свою цену, по сравнению с отечественными фирмами. Но, с другой стороны, отечественные фирмы в большей мере индексируют свои цены, что говорит о большей инерционности инфляции отечественных цен. Если следовать интерпретации ^аП et а1., 2001) коэффициентов при первом лаге инфляции и при ожидаемой инфляции в кривой Филипса, то фирмы, производящие отечественные промежуточные
§ 3
и »
н в
и ^
« в
S ^ у ^
Й с
о ^ д §
° I
S S
3
и
и &
13
§
а
я ^
Й в ft ^
° к
К °
о в а ^ & S
о
С <
и &
л u о К ft о
5
6
<
о4
о4
ft о
й
«
о
о4
о4
а
U ft О
й
«
о
о
В &
« §
н
о
и ft О
и &
см
<N 00 00 ОП VO сл 1 1П CD О 00 00 in in on £ VO \D 1 1 in 00 m о 00 VO \D 00 00 c^ сл CD 00 m CD CK CD 00 r^
' 1 го CD VO CD CD CD CD CD CD CD CD CD CD CD CD
vo \D 00 <N in in сл 1 on 00 00 CK VO in \D <N on vo <N on 1 1 <N CD <N CK in in 00 in CD 00 CD <N CD CD
о ' 1 CD <N CD CD CD CD CD CD CD CD CD CD CD CD
ск in on ск <N 00 on сл 1 CK 00 CK F^ in О ^f CK \D \D on ^f 1 1 m in <N m <N m сл 00 00 CD 00 <N CD \D CD m
' 1 <N CD CD CD CD CD CD CD CD CD CD CD CD CD
00 с^ on <N <N сл 1 \D <N on \D \D 1 1 VO О F^ <N in <N сл VO 00 CD VO <N CD 00 in CD
' 1 <N CD CD CD CD CD CD CD CD CD CD CD CD CD
00 о 00 00 on ОП VO сл 00 <N сл vo <N <N ^r 00 in in ^f VO on CK 00 \D <N <N 00 r^ in 00 00 CK in 00 00 c^ сл CK CD m m CD m c^ CD <N
СП CD CD VO CD CD CD CD CD CD CD CD CD CD CD CD CD CD
О m <N <N 00 ,_, <N VO on CK ,_, VO VO
00 m <N <N m о cn \D \D 00 <N О
0^ ■—i сл 4D \D en '—i CD CD 00 CD CD CD CD
CD ' 1 CD CD r^ CD CD CD CD CD CD CD CD CD CD CD CD CD CD
m CK in vo m VO m ^r О on 00 in 00 <N
in m <N <N c^ CK on 00 <N \D
en c^ 00 \D c^ \D l/n on '—i li^ сл CD CD CD '—i
' 1 <N CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD
m CK CK in 00 О CK 00 CK on VO on CK О
00 00 <N О \D in <N о in 00 <N in
'—1 r^ oo CI l/^ on i/^ c^ C^ CD '—i
r^ CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD
in
in in in in in in in in
CI г- '—i '—1 in '—i '—1 '—i '—1 '—1 '—i '—i '—i '—i ж
CD CD CD CD CD CD CD CD CD CD <N CD CD CD <N <N <N <N
in in in in in in in
li^ li^ l/n 1/П l/n li^ li^ 00 00 00 CD CD CD CD
i <N CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD CD
ft о
К к
й н
U
М
й н и М
ft о К
Й й й й й й
H H H H H H
u u u u u 1)
M M M M M W
й й й й й
^ ^ ^ ^ S
S S S S S
Й Й Й Й Й
и и и и и
u U U U U
о о о о о
К К К К К
н н н н н
й ft й й й й ft й ft й ft й ft
ю U М U М U М ю ю ю ю
О О О О О
&
oj^ ojx «x ^ ;C ;C ^
а а а Ь
u &
о Й ft
is
о о
В
о
н
с «
О К л
5
6
и £
товары, в большей мере формируют назадсмотрящие ожидания по сравнению с импорти- § рующими фирмами. §
Апостериорные оценки параметров автокорреляции свидетельствуют в пользу достаточ- ^ но высокой инерционности экзогенных процессов. При этом наибольшее значение автокорреляции (0.93) соответствует экзогенному процессу для премии за риск к отечественным активам. Для данного процесса также высоко значение стандартного отклонения. Как отмечалось выше, в настоящей модели шоки премии за риск к отечественным активам рассматриваются как некоторый обобщенный фактор финансовых шоков, потенциально влияющих на российскую экономику. Высокие значения параметра автокорреляции и стандартного отклонения могут быть обусловлены недостаточно детализированным описанием трансмиссионного механизма влияния шоков премии за риск на экономику, отсутствием финансовых несовершенств в модели.
3.3. Сравнение альтернативных спецификаций модели
Приведенная в настоящей статье модель насчитывает достаточно большое количество параметров, номинальных и реальных «жесткостей», определяющих траектории рассматриваемых переменных. Естественным образом возникает вопрос, насколько необходим данный набор механизмов для описания динамики российских макроэкономических переменных, и можно ли снизить размерность модели без значительных потерь в ее описательных свойствах.
Для ответа на данный вопрос было проведено сравнение исходной модели с рядом других спецификаций, в которых оценивалась модель без некоторых жесткостей, т. е. в качестве альтернативных спецификаций рассматривались модели, вложенные в исходную. При этом применялся стандартный подход при байесовском оценивании DSGE моделей (см., например, (Adolfson et al., 2007; Rabanal, Rubio-Ramírez, 2005; Smets, Wouters, 2007)) и в качестве критерия использовалась маржинальная плотность распределения данных условно на спецификации модели M:
p(Тт | M) = JL(XT | ©м,M)p(©м )d©м .
Маржинальная плотность оценивалась в рамках приближения Лапласа около моды апостериорного распределения. Отказ от привычек в потреблении либо от издержек на установку капитала значительно снижает маржинальную плотность распределения данных по сравнению с базовой спецификацией модели, что говорит о важности данных «жесткостей» при описании динамики рассматриваемых в настоящей работе переменных. Отказ же от спецификации с издержками на загрузку капитала приводит к существенному увеличению маржинальной плотности, поэтому нет необходимости в данной «жесткости» при описании динамики рассматриваемых переменных.
Что касается параметров номинальных жесткостей, то отказ от индексации заработных плат снижает маржинальную плотность, отказ от индексации цен — увеличивает, отказ от индексации цен и зарплат — увеличивает в меньшей мере, чем только отказ от индексации заработных плат, отказ от «жесткости» цен и зарплат — снижает, отказ от «жесткости» зарплат — снижает, отказ от «жесткости» цен — снижает. Таким образом, данные эксперимен-
ты свидетельствуют в пользу важности «жесткостей» цен и зарплат и индексации зарплат при описании данных, в индексации же цен, вероятно, нет необходимости.
В качестве наилучшей спецификации модели была выбрана спецификация с постоянной загрузкой капитала и без индексации цен. Здесь также была протестирована рассматриваемая спецификация модели по сравнению с альтернативными спецификациями, в которых модель дополнительно упрощалась. В результате был сделан вывод, что данная модель превосходит конкурирующие согласно критерию маржинальной плотности распределения данных. Для данной спецификации, так же как и для исходной версии модели, была проведена оценка апостериорного распределения с помощью методов Монте-Карло по схеме марковских цепей. Соответствующие оценки параметров приведены в табл. 1 в столбцах под заголовком «Оптимальная спецификация». При дальнейшем анализе будут использоваться результаты оценки модели без индексации цен и при постоянной загрузке капитала.
3.4. Качество описания модели данных
Зададимся вопросом о том, насколько модель способна объяснить динамику рассматриваемых макроэкономических рядов российской экономики. Для оценки качества модели будем рассматривать коэффициент детерминации Я . В таблице 2 представлены коэффициенты детерминации для всех используемых в эконометрическом анализе макроэкономических переменных. В качестве модельных переменных, которые ставятся в соответствие с данными, понимаются временные ряды сглаженных значений соответствующих переменных в DSGE модели. Сглаженным значением некоторой переменной в момент времени ^ называется условное математической ожидание данной переменной на основе всей доступной информации Тт (подробнее см., например, (НатШоп, 1994)). Таким образом, разница между модельной переменной и наблюдаемой переменной будет соответствовать сглаженным значениям ошибок измерения.
Таблица 2. Оценки качества аппроксимации
Выпуск Потребление Инвестиции Импорт Относительные цены Зарплата
0.98 0.81 0.87 0.94 0.61 0.74
Как следует из табл. 2, модель в целом обеспечивает достаточно хорошую согласованность с данными. Наибольшее значение коэффициента детерминации (0.98) соответствует динамике выпуска. Для потребления, инвестиций и импорта коэффициент детерминации варьируется от 0.81 до 0.94. Несколько хуже модель объясняет динамику относительных цен и реальных зарплат, коэффициент детерминации для этих переменных равен 0.61 и 0.74 соответственно.
Другой характеристикой качества модели является способность модели воспроизводить выборочные кросс-корреляции между переменными. На рисунке 1 в Приложении приводится сравнение эмпирических и модельных кросс-корреляций. Чтобы учесть апостериорную неопределенность в параметрах и неопределенность, связанную с короткой продолжительностью рассматриваемых в настоящей работе временных рядов, проводились эксперименты Монте-Карло, в рамках которых для каждого вектора параметров из апостериорного рас-
пределения, полученного на этапе МСМС, рассматривались симуляции модели на перио- § де длиной 1072 наблюдения. Данные симуляции проводились в предположении отсутствия § ошибок измерения. Далее в каждом эксперименте отбрасывались первые 1000 наблюдений, ^ в предположении, что начальные условия «забыты» и распределение вектора переменных модели уже стало эргодическим, и на оставшихся 72 точках симулированных рядов рассчитывались выборочные кросс-корреляции.
На рисунке 1 в Приложении представлены медиана и границы 95%-ного доверительного интервала, полученного в рамках экспериментов Монте-Карло распределения кросс-корреляций. В целом модель достаточно хорошо воспроизводит эмпирические данные. По большему набору переменных выборочные кросс-корреляции лежат в пределах доверительного интервала и демонстрируют близкие тенденции в сравнении с медианными кросс-корреляциями модели. Наихудшее соответствие модели и данных проявляется при сопоставлении кросс-корреляций инфляции номинального обменного курса с другими переменными, используемыми в эконометрическом анализе. Так, согласно результатам модели, ортогональные шоки обменного курса оказывают положительное влияние на деловую активность в отечественной экономике, поэтому должна наблюдаться положительная кросс-корреляция с выпуском, потреблением и инвестициями.
Фактически же центральный банк допускал девальвацию рубля во время кризисов 1998 и 2008 гг., и на фоне негативного влияния на отечественную экономику других шоков происходило одновременное падение выпуска, потребления и инвестиций наряду с ростом номинального обменного курса, что приводило к отрицательной эмпирической кросс-корреляции. Как отмечалось в разделе 1, описание денежно-кредитной политики с помощью экзогенного стохастического процесса для номинального обменного курса является весьма сильным упрощением, и обусловлено вероятной несимметричностью реакции центрального банка на шоки. Банк, вероятно, сдерживал укрепление рубля при положительных шоках внешнеэкономических условий и допускал девальвацию при отрицательных шоках.
3.5. Историческая декомпозиция вариации
Перейдем к анализу исторической декомпозиции вариации рассматриваемых в настоящей работе переменных, графики декомпозиции для которых представлены на рис. 2-7 в Приложении. Как следует из оцененного разложения, роль шоков экспорта в объяснении динамики российских макроэкономических переменных с течением времени значительно усилилась. Данные шоки обуславливали значительный рост выпуска, потребления, инвестиций и импорта относительно тренда перед кризисом 2008 года. Реальное укрепление национальной валюты в данный период времени, которое выражается в снижении относительных цен импортных товаров по отношению к отечественным товарам, также в значительной мере связано с шоками экспортных доходов.
С другой стороны, во время кризиса 2008 года произошло снижение экспортных доходов, которое было вызвано как падением цен на нефть, так и снижением выручки от других составляющих российского экспорта. Такое падение в значительной мере обуславливало уменьшение выпуска, инвестиций, потребления, импорта и ослабление национальной валюты. Рассматриваемый шок в меньшей мере отразился на падении реальных заработных плат, что связано с высокой степенью жесткости данного ценового показателя.
Другим важным фактором, обуславливающим динамику российских макроэкономических переменных во время кризиса 2008 года, является шок премии за риск к отечественным активам, а в случае кризиса 1998 года данный шок является определяющим фактором. В обоих кризисах модель идентифицирует рост рисков и соответствующее негативное влияние данного фактора на отечественную экономику. Это является достаточно интерпретируемым результатом, т. к. кризис 2008 года соответствует мировому финансовому кризису, а кризис 1998 года — дефолту по государственным долговым обязательствам российской экономики.
Положительные шоки обменного курса, т. е. девальвации рубля, смягчили падение российской экономики как во время кризиса 1998 года, так и во время кризиса 2008 года. Результаты эконометрического анализа свидетельствуют в пользу того, что шоки совокупной факторной производительности обуславливают малую долю вариации российских макроэкономических переменных при анализе факторов делового цикла. Но данный вывод не противоречит гипотезе о том, что динамика трендовой составляющей макроэкономических переменных могла в большей мере обуславливаться именно ростом производительности.
3.6. Функции импульсного отклика
В заключение описания результатов приведем функции импульсного отклика рассматриваемых в эконометрическом анализе переменных на структурные шоки. На рисунке 8 в Приложении представлены функции импульсного отклика на положительный шок экспорта величиной в одно стандартное отклонение. Данный шок является положительным трансфертом со стороны внешнего сектора и увеличивает богатство отечественных домохо-зяйств, что приводит к увеличению спроса как на отечественные, так и на импортные товары. Положительный шок экспорта в 8.8% приводит к увеличению выпуска отечественного промежуточного продукта с пиком роста порядка 1.5%, увеличению потребления с пиком роста около 0.8%, увеличению инвестиций с пиком роста 6%, увеличению импорта с пиком роста примерно 4%, падению относительных цен импортных товаров к отечественным товаром (т. е. реальному укреплению отечественной валюты) с низшей точкой падения порядка 2.5% и росту реальных зарплат с пиком около 0.5%.
На рисунке 9 в Приложении представлены функции импульсного отклика на положительный шок производительности размером в одно стандартное отклонение. Шок производительности ведет к увеличению кривой производственных возможностей отечественной экономики, и, соответственно, к росту выпуска отечественных промежуточных товаров и потребления домохозяйств. Рост СФП приводит к увеличению предельного продукта капитала и увеличению рентной цены капитала в среднесрочной перспективе, что ведет к росту теневой стоимости единицы капитала и росту инвестиций. Рост СФП также приводит к снижению предельных издержек фирм, производящих отечественный промежуточный продукт, что ведет к снижению их цен (относительные цены импортных товаров к отечественным товарам растут) и, соответственно, замещению импорта отечественными товарами в краткосрочной перспективе (импорт падает). После истечения одного года домохозяйства начинают потреблять накопленные активы, и отклонение импорта от долгосрочного равновесия становится положительным.
На рисунке 10 в Приложении представлены функции импульсного отклика на положительный шок премии за риск в одно стандартное отклонение. Шок премии за риск делает
менее привлекательными вложения в отечественные ценные бумаги по сравнению с зару- § бежными активами. При неизменной процентной ставке на мировом финансовом рынке § в отечественной экономике наблюдается рост процента. В свою очередь, рост процентных ^ ставок снижает компоненты спроса, чувствительные к проценту: потребление домашних хозяйств и инвестиции, что далее приводит к падению выпуска отечественных промежуточных товаров и импорта. Происходит снижение предельных издержек и ослабление реального обменного курса.
На рисунке 11 в Приложении представлены функции импульсного отклика на положительный шок обменного курса в одно стандартное отклонение. Данный шок соответствует удорожанию импортных товаров и замещению импорта отечественными товарами, что приводит к увеличению выпуска отечественных промежуточных товаров. Также шок обменного курса снижает реальные заработные платы в течение достаточно продолжительного времени, поскольку номинальные заработные платы — «жесткие». Соответственно, фирмы увеличивают объем используемого труда вследствие превышения предельного продукта труда над реальной заработной платой, что вносит дополнительный вклад в рост выпуска. Происходит рост предельного продукта капитала и, соответственно, рост инвестиций. Рост количества отработанных часов приводит к росту трудовых доходов и, соответственно, росту потребления, что в свою очередь повышает спрос на импортные товары. Таким образом, наблюдается небольшое падение импорта в течение первых двух кварталов после шока с последующим его увеличением. Эффект импортозамещения ослабляется из-за эффекта неполного переноса обменного курса в цены импортных товаров.
4. заключение
В настоящей работе описана динамическая стохастическая модель общего равновесия и проведена ее оценка с помощью байесовских эконометрических методов для российской экономики. В целом эконометрическое оценивание приводит к достаточно непротиворечивым и интерпретируемым результатам. Был протестирован ряд альтернативных спецификаций, проанализировано, какие «жесткости» в модели необходимы при описании динамики российских макроэкономических переменных. В качестве «наилучшей» спецификации модели была выбрана спецификация с постоянной загрузкой капитала без индексации цен импортных и отечественных промежуточных товаров.
Оценки структурных параметров модели лежат в приемлемом диапазоне и согласуются с результатами исследований по построению и оценке DSGE моделей, приведенными в экономической литературе. Результаты оценивания говорят о большей жесткости номинальных заработных плат по сравнению с жесткостью цен. Средняя продолжительность контракта номинальной заработной платы в «наилучшей» спецификации составляет 3.7 квартала, а средняя продолжительность ценовых контрактов для импортных цен и цен на отечественные промежуточные товары составляет 3.0 и 1.7 квартала соответственно.
Оцененная модель достаточно хорошо согласуется с данными и способна воспроизводить циклические компоненты динамики российских макроэкономических переменных. Модель идентифицирует рост рисков, как во время кризиса 2008 года, так и во время кризиса 1998 года, и свидетельствует в пользу значительного влияния данного фактора на российскую экономику. Это, в свою очередь, может являться аргументацией для дальнейшего
расширения модели путем детализации финансового сектора. Шоки совокупной факторной производительности занимают незначительное место в объяснении делового цикла российской экономики. Шоки экспорта обуславливали значительный рост выпуска, потребления, инвестиций, импорта и отечественных цен по отношению к импортным ценам в отклонениях от тренда перед кризисом 2008 года, и падение соответствующих переменных во время и после кризиса 2008 года.
Список литературы
Казакова М. В., Синельников-Мурылев С. Г. (2009). Конъюнктура мирового рынка энергоносителей и темпы экономического роста России. Экономическая политика, 5, 118-135.
Карев М. Г. (2009). Инфляция, реальный обменный курс и денежная политика в экономике с ограниченной эластичностью потока капитала по процентной ставке. Экономический журнал ВШЭ, 13 (3), 329-359.
Мельников Р. М. (2010). Влияние динамики цен на нефть на макроэкономические показатели российской экономики. Прикладная эконометрика, 17 (1), 20-29.
Adjemian A., Bastani H., Juillard M., Mihoubi F., Perendia P., Ratto M., Villemot S. (2011). Dynare: Reference manual, version 4. Dynare Working Papers, No. 1, CEPREMAP.
Adolfson M., Laseen S., Linde J., Villani M. (2007). Bayesian estimation of an open economy DSGE model with incomplete pass-through. Journal of International Economics, 72 (2), 481-511.
Aguiar M., Gopinath G. (2007). Emerging market business cycles: The cycle is the trend. Journal of Political Economy, 115 (1), 69-102.
Brooks S. P., Gelman A. (1998). General methods for monitoring convergence of iterative simulations. Journal of Computational and Graphical Statistics, 7 (4), 434-455.
Calvo G. A. (1983). Staggered prices in a utility-maximizing framework. Journal of Monetary Economics, 12 (3), 383-398.
Christiano L. J., Eichenbaum M., Evans C. (2001). Nominal rigidities and the dynamic effect of a shock to monetary policy. NBER Working Paper, 8403.
Christiano L. J., Eichenbaum M., Evans C. (2005). Nominal rigidities and the dynamic effect of a shock to monetary policy. Journal of Political Economy, 113 (1), 1-45.
Del Negro M., Schorfheide F. (2008). Forming priors for DSGE models (and how it affects the assessment of nominal rigidities). Journal of Monetary Economics, 55 (7), 1191-1208.
Erceg C. J., Henderson D. W., Levin A. T. (2000). Optimal monetary policy with staggered wage and price contracts. Journal of Monetary Economics, 46 (2), 281-313.
Gali J., Getler M., Lopez-Salido D. (2001). European inflation dynamics. European Economic Review, 45 (7), 1237-1270.
Ireland P. N. (2004). A method for taking models to the data. Journal of Economic Dynamics and Control, 28 (6), 1205-1226.
Ito K. (2009). The Russian economy and the oil price: A co-integrated VAR approach. Transitory Studies Review, 16 (1), 220-227.
Lucas R. E. (1976). Econometric policy evaluation: A critique. Carnegie-Rochester Conference Series on Public Policy, 1, 19-46.
Mendoza E. G. (1995). The terms of trade, the real exchange rate, and economic fluctuations. Interna- s tional Economic Review, 36 (1), 101-137.
Prescott E. C. (1986). Theory ahead of business cycle measurement. Carnegie-Rochester Conference Series on Public Policy, 25, 11-44. ®
Rabanal P., Rubio-Ramirez J. F. (2005). Comparing new Keynesian model of the business cycle: A Bayes-ian approach. Journal of Monetary Policy, 25 (6), 1151-1166.
Sims C. A. (2002). Solving linear rational expectations models. Computational Economics, 20 (1-2), 1-20.
Smets F., Wouters R. (2003). An estimated dynamic stochastic general equilibrium model of the Euro Area. Journal of European Economic Association, 1 (5), 1123-1175.
Smets, F., Wouters, R. (2007). Shocks and frictions in US business cycles: a Bayesian DSGE approach. American Economic Review, 97 (3), 586-606.
Sosunov K., Zamulin O. (2006a). The inflationary consequences of real exchange rate targeting via accumulation of reserves. BOFITDiscussion Papers, 11/2006.
Sosunov K., Zamulin O. (20066). Can oil prices explain the real appreciation of the Russian ruble in 1998-2005? CEFIR Working Papers, w0083.
Uribe M., Yue V Z. (2006). Country spreads and emerging countries: Who drives whom? Journal of International Economics, 69 (1), 6-36.
Приложение
Выпуск^ - h) Потребление^ - h) Инвестиции^ - h) Импорт^ - h)
Цены(1 - h) Зарплата^ - h)
Выпуск^)
Потребление^)
Инвестиции^)
Импорт^)
Цены(0
Зарплата^)
Экспорт^) 0
-1 0.5
Инфляция курса^) о -0.5
Экспорт^ - h) Инфляция курса^ - h)
0.5г
10 0 5 h
Рис. 1. Сравнение кросс-корреляций модели с данными. Сплошная толстая линия — данные, сплошная тонкая линия — медиана на основе симуляций модели, штриховые линии — 95%-ный доверительный интервал на основе симуляций модели
«О 40 t^ 00 OS
о\ ^ ^ о\ о\ о\ ^ ^ о\ о\
: Начальные условия
■ ■ ■ ■ ■ Ошибка
.... Шок обменного курса
___Шок СФП
^ Шок премии за риск
1 Шок экспорта ■ Потребление
<N <N <N <N <N <N <N <N <N <N <N <N <N
Рис. 2. Историческая декомпозиция вариации реального потребления
Начальные условия ■ I ■ I ■ Ошибка * и 1 Шок обменного курса ■ ■ 'Шок СФП
Шок премии за риск I Шок экспорта ■Инвестиции
оооооооооооооооооо
оо с^ ^ С\ ^ ^ ^ С\ ^ ^
<N <N <N <N <N<N<N<N<N <N<N<N<N
Рис. 3. Историческая декомпозиция вариации реальных инвестиций
-25 л
а а СУ СУ СУ СУ СУ СУ а а а а а а а а су СУ
«о 40 с- 00 о <N m -t «о 40 с- 00 о <N
о о О о о о о о о о .—1 .—1 .—1
о <N о <N О <N о <N о <N о <N о <N о <N о <N о <N о <N о <N О <N
Рис. 4. Историческая декомпозиция вариации реального выпуска
Начальные условия ■: ■: ■ Ошибка р ' ' Шок обменного курса ^ ^ - Шок СФП
\ Шок премии за риск ' Шок экспорта ■Импорт
л г
CQ
00
^ ^ С\ ^ ^ ^ С\ ^
<N<N<N<N<N<N<N<N<N<N<N<N
Рис. 5. Историческая декомпозиция вариации реального импорта
Начальные условия ■' ■' ■ Ошибка г ж , Шок обменного курса -|-|-| Шок СФП
Шок премии за риск . Шок экспорта ■ Относительные цены
а а а с> о а а а С? а С? С? а СУ СУ СУ СУ СУ
«о 40 с- 00 сл о .—1 <N m -t «о 40 с- 00 с^ о .—1 <N
0\ 0\ о о о о о о О о о о .—1 .—1 .—1
сл сл 0\ 0\ сл о <N о <N о <N о <N о <N о <N О <N о <N о <N о <N о <N о <N О <N
Рис. 6. Историческая декомпозиция вариации относительных цен импортных к отечественным товарам
^^^ Начальные условия ■. ■. ■ Ошибка
Шок обменного курса тип Шок СФП
"Ъ Шок премии за риск I Шок экспорта ■Реальные зарплаты
оооооооооооооооооо
<N<N<N<N <N<N<N<N
<N <N <N <N
Рис. 7. Историческая декомпозиция вариации реальных зарплат
2 1.5 1
0.5 0
Выпуск
о r^i ю СЛ
Импорт
1.5
0.5
Потребление
Инвестиции
О r^i ю 0\
Относительные цены
0 -1
-2 -3 -4
Зарплата
0.5
-0.5
О r^i ю Оу
Рис. 8. Функции импульсного отклика на шок экспорта (в процентных отклонениях от долгосрочного равновесия). Сплошная линия — апостериорное среднее, штриховые линии — 90%-ный апостериорный доверительный интервал
1
0
1
6
0
2.5 2 1.5 1
0.5 0
Выпуск
Потребление
Инвестиции
0.5
о съ о\
О r^i ю Оу
6 4 2 0 -2
Импорт
Относительные цены
2.5 2 1.5 1
0.5 0
Зарплата
о съ сл
Рис. 9. Функции импульсного отклика на шок производительности (в процентных отклонениях от долгосрочного равновесия). Сплошная линия — апостериорное среднее, штриховые линии — 90%-ный апостериорный доверительный интервал
0
Выпуск
Потребление
Инвестиции
1
0 -1
-2 -3 -4
0
-0.5 -1 -1.5 -2 -2.5
10 5 0
-5 -10 -15
Импорт
Относительные цены
1
0.5 0
-0.5
Зарплата
Рис. 10. Функции импульсного отклика на шок премии за риск (в процентных отклонениях от долгосрочного равновесия). Сплошная линия — апостериорное среднее, штриховые линии — 90%-ный апостериорный доверительный интервал
1
0.5 0
-0.5
Импорт
Относительные цены
о m ск
0.5 0
-0.5 -1 -1.5 -2 -2.5
Зарплата
Рис. 11. Функции импульсного отклика на шок обменного курса (в процентных отклонениях от долгосрочного равновесия). Сплошная линия — апостериорное среднее, штриховые линии — 90%-ный апостериорный доверительный интервал
3
2
0