В ряде исследований американских ученых для моделирования динамики процентных ставок используется модель авторегресии–проинтегрированного скользящего среднего Дж.Бокса–Г.Дженкинса (ARIMA)[70]. В частности, П.Кэмпбелл и Р.Шиллер использовали модель ARIMA для описания колебаний долгосрочных процентных ставок[71], Е.Фама и Р.Блисс – для прогнозирования изменений краткосрочных процентных ставок[72], Н.Галтекин и Р.Рогальски – для прогнозирования доходностей бескупонных облигаций[73]. Мы считаем,что моделью ARIMA можно воспользоваться и при построении сценариев изменения значений главных компонент временной структуры процентных ставок.
Модель ARIMA предназначена для описания и прогнозирования динамики нестационарных временных рядов, характеризующихся нестабильным средним значением уровней ряда. Как правило, в ходе анализа рядов динамики процентных ставок можно выявить несколько трендов, последовательно сменяющих друг друга и определяющих движение процентных ставок в течение некоторого промежутка времени. Высокая степень зависимости между элементами ряда, обусловленная этими трендами, находит выражение в высоких значениях выборочного коэффициента автокорреляции для больших лагов. В то же время переход к первым разностям уровней ряда динамики процентной ставки позволяет привести его к стационарному виду: устраняются тренды, стабилизируется среднее значение, а выборочная автокорреляционная функция приобретает затухающую форму. Это свидетельствует о том, что динамика процентных ставок определяется интегрированным нестационарным случайным процессом, который можно описать при помощи модели ARIMA.
Модель ARIMA(p,d,q) задает процесс изменения значений случайной переменной при помощи небольшого числа параметров: степени интегрирования d, p коэффициентов авторегресии и q коэффициентов скользящего среднего. Степень интегрирования d равна числу шагов расчета разностей между последовательными элементами временного ряда, необходимому для приведения исходного ряда к стационарному виду. Полученный стационарный ряд Xt описывается при помощи модели
, (2.3.1)где m – константа (опускаемая в большинстве моделей как незначимо отличающаяся от нуля), fi – коэффициенты авторегрессии, qj – коэффициенты скользящего среднего, et – независимо распределенные случайные ошибки.
Математический аппарат, разработанный Боксом и Дженкинсом, позволяет идентифицировать число порядков p, d и q, оценивать параметры fi и qj, а также строить условные распределения будущих значений уровней исходного временного ряда. Как отмечают С.А.Айвазян и В.С.Мхитарян, в наиболее распространенных моделях ARIMA(p,d,q) используются комбинации порядков (0,1,1), (0,2,2), (1,1,1), (1,1,0) и (2,1,0)[74].
Автор предлагает формировать сценарии будущих перемещений временной структуры процентных ставок, взяв за основу квантили условных распределений будущих значений главных компонент, построенные при помощи модели ARIMA. Для того, чтобы используемые квантили адекватно отражали спектр возможных состояний временной структуры процентных ставок, они должны представлять различные участки распределения, отражать его симметрию, а среднее значение и стандартное отклонение выборки квантилей должны соответствовать прогнозу модели ARIMA и его стандартной ошибке.
Поскольку главные компоненты независимы друг от друга, сценарии временной структуры строятся исходя из всех возможных комбинаций квантилей условного распределения будущих значений главных компонент. В частности, при использовании двухфакторной модели временной структуры процентных ставок и пяти сценариев будущих значений каждой главной компоненты общее число сценариев перемещения временной структуры процентных ставок оказывается равным 25.
Каждый сценарий состояния временной структуры процентных ставок через заданный промежуток времени определяет множество цен облигаций, которые установятся на рынке при реализации данного сценария. Располагая информацией о структуре портфеля, можно рассчитать значение его будущей стоимости при реализации каждого из сценариев по формуле
, (2.3.2)где FVs(m) – стоимость портфеля через промежуток времени m при реализации сценария перемещения временной структуры процентных ставок s, qj – число облигаций выпуска j, включенных в состав портфеля,
– цена облигации выпуска j через промежуток времени m при реализации сценария перемещения временной структуры процентных ставок s.Используя выборку значений рыночной стоимости портфеля FVs(m) при различных сценариях перемещения временной структуры процентных ставок, можно построить выборку значений доходности портфеля h(m), которая рассчитывается по формуле
, (2.3.3)где PV – рыночная стоимость портфеля в начальный момент времени.
Формулы (2.3.2) и (2.3.3) можно применять лишь при условии, что в течение периода вложений по портфелю не осуществляется никаких выплат. В противном случае будущая стоимость портфеля определяется не только размером начальных вложений в облигации и состоянием временной структуры процентных ставок через промежуток времени m, но также стратегией реинвестирования и состояниями временной структуры процентных ставок в моменты реинвестирования.
Если не использовать дополнительных предположений, построение репрезентативной выборки значений доходности портфеля требует построения траекторий движения временной структуры процентных ставок в течение всего периода вложений, а не сценариев ее состояния на конец периода. Для решения этой задачи необходимо воспользоваться техникой имитационного моделирования, задавая последовательность изменений значений главных компонент на основе оцененной модели ARIMA и значений ошибок et, сгенерированных при помощи датчика псевдослучайных чисел. Однако такой подход к оценке параметров распределения доходности портфеля требует колоссального объема вычислений и затрат машинного времени.
Автор считает,что в случае, когда стоимость платежей, которые должны быть получены в течение периода вложений, существенно меньше общей стоимости портфеля, целесообразно воспользоваться альтернативным методом. Пусть все купонные платежи, выплачиваемые по облигациям выпуска j, реинвестируются путем приобретения дополнительных облигаций выпуска j. Тогда доходность вложений в каждую облигацию определяется по формуле
, (2.3.4)где
– доходность облигации выпуска j за период времени m при реализации сценария перемещения временной структуры процентных ставок s, – цена облигации выпуска j через промежуток времени m при реализации сценария перемещения временной структуры процентных ставок s, – коэффициент увеличения размера позиции по облигации выпуска j в результате реинвестирования полученных купонных платежей, – цена облигации выпуска j в начальный момент времени.Поправочный коэффициент
зависит не только от размеров и сроков выплаты купонных платежей по облигации, но и от состояния временной структуры процентных ставок в момент реинвестирования, которое определяет цену покупки новых облигаций. Используя предположение об устойчивости динамики временной структуры процентных ставок, можно определить , (2.3.5)где Ct – размер купона, выплачиваемого через промежуток времени t,
– цена облигации выпуска j через промежуток времени t при реализации сценария перемещения временной структуры процентных ставок, который определяется той же комбинацией квантилей распределения будущих значений главных компонент, что и сценарий состояния временной структуры процентных ставок в конечный момент времени.Полученная выборка доходностей облигаций позволяет оценить математическое ожидание E(h(m)) и среднеквадратическое отклонение s(h(m)) распределения доходности портфеля за период m по формулам
, (2.3.6) , (2.3.7)где xj – доля вложений в облигации выпуска j в рыночной стоимости портфеля в начальный момент времени, S – число сценариев перемещения временной структуры процентных ставок, J – число выпусков облигаций, включенных в состав портфеля.
Рис.2.3.1. Методика сценарного анализа процентного риска
портфеля государственных облигаций.
Методика сценарного анализа процентного риска, разработанная диссертантом, дает возможность ответить на ряд вопросов, имеющих как прикладное, так и теоретическое значение. Во-первых, она позволяет измерить ожидаемую доходность и риск портфелей государственных облигаций и сопоставить их с характеристиками альтернативных объектов вложений. Во-вторых, она позволяет оценить характер соотношения между доходностью и риском для различных портфелей облигаций и определить структуру эффективных портфелей, обеспечивающих наибольшую ожидаемую доходность при заданной степени риска. В-третьих, она позволяет выяснить, как изменяются значения показателей доходности и риска при увеличении срока вложений инвестора.