УДК 336.763
Адаптивное управление риском позиции на основе фьючерсных контрактов
Аннотация. Страхование или хеджирование позиции является одной из основных задач портфельного менеджмента. Традиционный метод оценки оптимального числа фьючерсных контрактов предполагает оценку регрессии изменения цены спот на изменения фьючерсной цены [1, 4, 6].
В статье рассматривается динамическое страхование открытой позиции фьючерсными контрактами с учетом неоднородности изменений цен как по уровню, так и по дисперсии. При этом нестационарность корреляционной связи точно так же принимается во внимание. Для прогноза волатильности цен используется простая модель динамики изменений, основанная на экспоненциальном сглаживании, которая позволяет прогнозировать ожидаемую доходность и волатильность временных рядов.
Ключевые слова: экспоненциальное сглаживание; фьючерсные контракты; коэффициент хеджирования; динамическое страхование позиции.
Abstract. Insurance or hedging a position is one of the main tasks of portfolio management. The traditional method for estimating the optimal number of futures contracts involves an assessment of the regression of price changes on the spot futures price [1, 4, 6]. In this paper we consider the dynamic insurance of open positions in futures contracts, taking into account heterogeneity of price changes in both level and dispersion. Nonstationarity of correlation is also taken into account. A simple model of the dynamics of changes based on exponential smoothing used for forecasting allows predicting the expected return and volatility of time series.
Keywords: exponential smoothing; futures contracts; hedge ratio; dynamic insurance of a position.
Калинин А. Т.,
студент магистратуры Финансового университета Н andrey.kaLinin.maiL@gmaiL.com
E(t) = ax(t) + pE(t-1), р = 1 -а
(1)
1. Модель динамики волатильности
Экспоненциальное среднее временного ряда х можно определить равенством [3, 5]:
Е()= х()+ х( -1)+р2х(-2) + р3 х(-3)+ •••.
Здесь Ъ^) — значение экспоненциальной средней в момент времени ^ а — параметр сглаживания, 0 < а < 1, р = 1 - а. В этой формуле предполагается, что исходный временной ряд определен для всех моментов времени, предшествующих текущему моменту t. При этом веса ю. = арг экспоненциально спадают до нуля, и их сумма равна единице. Экспоненциальное среднее можно также определить рекуррентно согласно формуле:
Чем ниже значение а, тем меньший вес придается текущей цене и тем сильнее эффект сглаживания. В техническом анализе вместо вещественного параметра а используется временное окно w — целочисленный параметр, определяемый из равенства а = 2/(w + 1). В дальнейшем экспоненциальное среднее временного ряда x в момент времени t с окном w обозначается через E(t; w,x). В этом обозначении часть аргументов может отсутствовать, если это не вызывает недоразумений. Экспоненциальное сглаживание предполагает задание начальной оценки S (0), то есть оценки экспоненциальной средней на нулевой момент времени. Обычно в качестве такой оценки используют среднее арифметическое первых w значений временного ряда (w - размер временного окна).
Рассмотрим следующую модель временного ряда, учитывающую неоднородность по уровню и дисперсии:
х(t) = E(t;x,wj) + e(t), e2(t) = E(t;e2,w2) + e(t), (2)
здесь: x(t) — абсолютные или относительные изменения цен; случайная величина e(t) имеет нулевое среднее и ди-
Научный руководитель: Керимов А.К., кандидат физико-математических наук, доцент кафедры «Прикладная математика» Финансового университета.
сперсию, вообще говоря, зависящую от времени £ остатки е(г) имеют нулевое среднее и постоянную дисперсию; w1, w2 е параметры модели, подлежащие оценке.
Уравнения (2) можно представить в более наглядной форме, позволяющей определить прогнозы на один шаг вперед уровня исходного ряда х и его дисперсии. Используя равенство (1), получим:
х (,) = а1 х (,) + РЕ(, -1;х,1) + е (,).
Откуда легко получаем равенство: х(,) = Е(, - 1; х, w1) + е(,)/Рг Заменяя t на t+1, последнее равенство перепишем в виде:
х (( +1) = E (t; х ,w 1) +
+^=£ajP; +e ((+j)
i=0
Pi
Pi
(3)
=Sa 2p2 +
<(+j)
P2 '
(4)
xf ((+1) = M {x(t +1)| х (s),s < t известны } = = E ((;w 1, x),
(5)
(индекс t означает, что рассматривается условное математическое ожидание). Поскольку среднее случайной величины e(t + 1) равно нулю, то отсюда следует, что прогнозы дисперсии и стандартного отклонения на один шаг вперед определяются равенствами:
Df (t +1) = EMA (T;w2,e2), of ((+1) = ^EMA (T;w 2,e2)
(6)
Последнее равенство показывает, что динамика исходного ряда представляется авторегрессией бесконечного порядка с коэффициентами, зависящими от окна Точно так же показывается, что второе уравнение системы (2) можно представить в виде:
е2 ((+1)= Е ((;е2) + е (+^ =
2. Оценка параметров
Опишем метод оценки параметров w1, модели (3), (4) по эмпирическим ценовым данным 1 < t < Т. Для заданных w1, w2 проводится следующая последовательность вычислений. Ряд х определяется как абсолютные или относительные изменения ценовых данных, то есть:
(,) = Ур (,) =р (,)-р (, -1) или х (,)=Ур (,) / р (, -1).
Для заданных значений окон w1, w2 прогнозы уровня и стандартного отклонения изменений определяются согласно алгоритму:
а) вычислить квадраты отклонения доходностей от их среднего уровня согласно формуле:
е2(,) = (х(,)-ЕМА (^ 1,х))2,
б) в качестве оценки дисперсии доходностей принимается экспоненциальное среднее отклонений е2 с параметром сглаживания w2, то есть полагается:
Таким образом, уравнения (2) эквивалентны уравнениям (3) и (4).
Прогноз случайной величины х^ + 1), сделанный в текущий момент времени t, определяется как условное математическое ожидание этой случайной величины при условии, что вся предыстория процесса известна вплоть до текущего момента, то есть известны все значения х(б) для 5 < I Из представления (3) следует, что прогноз хД + 1) значения ряда х на один шаг вперед определяется равенством:
o
(t)= EMA (t;w 2,e2)
поскольку по условию случайная величина е(г) имеет нулевое математическое ожидание.
Из равенства (4), в силу предположения на остатки е(г), следует, что условное математическое ожидание е2^ + 1) совпадает с Е ^; w2 е2), то есть:
М, [е2 (+1) = (Е(2,е2),
В результате получаем ряд ст(г), t = 2,3, ...,Т, ст(г) -это оценка стандартного отклонения на момент времени I Прогнозы ожидаемого значения ряда х и его стандартного отклонения на один шаг вперед, сделанные в момент определяются равенствами (6).
Окна сглаживания w1, w2 можно выбирать из условия, что все значения х(г) лежат в коридоре:
[х, ((+1)-2а/ ((+1),х, ((+1) + 2а/ ((+1)]^ , = 2,3,...Г.
Периодически по мере поступления новых данных значения окон переоцениваются.
Проверка адекватности модели исходным данным заключается в проверке гипотез относительно остатков е2(г) и е(г), которые определяются равенствами (2). В частности, остатки должны удовлетворять условиям белого шума. Проверка этой гипотезы осуществляется стандартными методами.
i=0
Прогнозирование на основе экспоненциального сглаживания — простая и достаточно эффективная процедура, особенно в случае, когда неоднородности по среднему или дисперсии не слишком резкие. При возникновении скачка в ценах, а следовательно, и в изменениях цены оценки прогнозов довольно быстро адаптируются к новому уровню.
Пример 1. Динамика цен закрытия вместе с динамикой прогнозов волатильности на основе модели (2) приведена на рис. 1. Исходными данными являются ежедневные цены закрытия акций в период 29.12.1101.06.12 (всего 101 отсчет). Выбранная модель характеризуется окнами w1 = 18, = 22. Штриховые линии определяют прогнозную верхнюю и нижнюю границу коридора. Верхняя граница определяется равенством х^ + 1) + 2ст^ + 1), нижняя — х(г + 1) - 2ст^ + 1). Сплошная тонкая линия представляет относительные изменения цены, сплошная жирная линия — прогнозы относительных изменений цен на один шаг вперед. Из рисунка видно, что практически все изменения цен (96%) лежат в пределах установленных границ.
3. Оценка корреляции на основе экспоненциального сглаживания
В общем случае ковариация или корреляция между относительными изменениями доходностей двух различных активов заметно меняется во времени. Динамика корреляционной связи между относительными изменениями цен и фьючерсных контрактов на них, как правило, довольно высокая, порядка 0,60,7 в начале эмиссии серии фьючерсных контрактов, и растет по мере приближения к дате исполнения. За месяц до исполнения уровень корреляции имеет порядок 0,85-0,9 (см. рис. 2).
Ниже приводится краткое описание модификации известной схемы адаптивной оценки корреляции на основе экспоненциального сглаживания, приведенное, например, в работе [3]. Пусть х(^, у{И) — значения абсолютных или относительных изменений цен для t = 2, ..., Т. Задача заключается в получении прогноза ковариации и корреляции на момент времени t+1, по ценовым данным, известным до текущего момента t. Вначале методом, изложенным в предыдущем разделе, оцениваются средние уровни и стандартные отклонения этих временных рядов, тем самым окна №х1, wx2) и №у1, wу2), по которым определяются эти величины, можно считать известными. В случае относительных изменений цен спот и фьючерсных цен эти пары окон, как правило, совпадают в силу их сильной корреляции. Далее определяется временной ряд 52:
52 ()=(х ()-Шх ())(у ()-my ())
где
Шх () = ЕМА (г^хрх), Шу ()= ЕМА 1,у )-
средние временных рядов х и у на момент времени t соответственно. В качестве оценки ковариации на текущий момент времени принимается скользящее среднее ряда 52:
Соу {),у (?)} = ЕМА ^).
Параметр сглаживания w в этой формуле полагается равным максимуму из окон wх2, wу2. Уровень корреляции на текущий момент времени оценивается по формуле:
Согг{х у (?)} =
Соу{х (),у (?)}
х С)у (') '
где стх(^, сту(^ - оценки стандартных отклонений рядов х(^ и у{И) соответственно на момент времени ^ полученные методом предыдущего раздела. Прогнозы ковариации и корреляции на следующий момент времени определяется равенствами:
Соуг ( +1) = ЕМА (( ,s2), Соггг ( +1) = ^ЕМА ( ,^) .
Пример 2. Динамика корреляционной связи между относительными изменениями цен акций ОАО «Газпром» и фьючерсных контрактов на них в интервале 22.01.12 — 30.04.12 показана на рис 2. Окна сглаживания, по которым оценивалась корреляция, равны w1 = 18, w2 = 22. Как видно из рисунка, уровень корреляции в целом довольно высокий и в то же время заметно растет по мере приближения к дате исполнения фьючерсного контракта.
4. Оптимальное число фьючерсных контрактов
Рассмотрим портфель, содержащий 0, единиц определенного актива (длинная позиция) и к фьючерсных контрактов на этот актив, при этом каждый контракт содержит ц} единиц базового актива. Обозначим через и Е^) - цену спот и фьючерсную цену (в денежном выражении) на момент времени I Изменение стоимости портфеля за единицу времени представляется в виде:
УЖ ( +1) = ГО (+ 1)е5 + к VР (+1), (7)
где через VW(t), - обозначают абсолют-
ные изменения соответствующих величин за единицу времени, то есть:
УЖ( +1) = Р( +1)-Р(), ГО(+1) = = Б( +1)-Б(), УР( +1) = Р( +1)-Р().
Рис. 1. Прогнозы и доверительный коридор изменений цен в пределах 96%
Обычно, если не оговорено противное, в качестве единицы измерения времени выбираются сутки. Отметим, что вероятностные характеристики изменения рассматриваемых величин меняются во времени. В дальнейшем прогнозы вероятностных характеристик рассматриваемых случайных величин:
УЖ ( +1), УБ ( +1), У Г ( +1),
сделанные в момент времени t, определяются при условии, что история цен спот и фьючерсных цен известна вплоть до этого момента времени. Соответствующие прогнозы снабжаются индексом t. Например, через М{(УБ) обозначается условное математическое ожидание случайной величины УБ^ + 1), которое интерпретируется как прогноз изменения цены, сделанный в момент I Из равенства (7) следует, что:
М, (УЖ) = М, (У Б) + ку{ М (УБ) (8)
0,7 -
Рис. 2. Динамика корреляций между относительными изменениями цен акций ОАО «Газпром» и фьючерсных контрактов серии SPFB.GAZR-12.12
в, (УЖ) = (УБ) + кЦБ, (УГ) + +2Соу, ( УБ, куг У Г )_ = 02Б, (УБ) + кЦБ, (УГ) +
+2кОгд,Соу, (УБ,УГ). (9)
Задача страхования заключается в определении числа фьючерсных контрактов к в портфеле таким образом, чтобы изменчивость стоимости портфеля была минимальной, то есть к следует выбирать из условия минимума дисперсии случайной величины У^. Выражение для дисперсии (9) представляет квадратный трехчлен относительно к, минимум которого достигается при:
_ ОРоу, (У Б,У Г) _ ре5
к0
У{Б> (УГ) дг
СоМУБУГ). (10)
Б, (УГ) v ;
Поскольку число контрактов - целое число, полученное значение к0 необходимо округлить до целого
Рис 3. Динамика изменений стоимости портфеля при динамическом страховании (ряд 1) и незастрахованной позиции (ряд 2)
числа. Минимальное значение дисперсии получается подстановкой в равенство (9) найденного значения к0. Знак минус перед к0 показывает, что страхователь должен занять противоположные позиции по хеджируемому активу и фьючерсным контрактам.
Представление для коэффициента хеджирования (10) полностью совпадает с традиционным, однако при выводе не требуется никаких предположений об однородности временных рядов изменений цен спот и фьючерс. При этом ожидаемые значения изменений интерпретируются как прогнозы изменений, а дисперсии - как характеристики ошибки прогноза на будущий период. Эти величины оцениваются в текущий момент времени по известным на этот момент истории цен методом, изложенным выше. По мере поступления данных прогнозы и их дисперсии корректируются.
5. Динамическое управление риском портфеля
Пусть t - текущий момент времени. Состояние портфеля на текущий момент определяется парой (О, к), где 0, > 0 - количество базового актива (предполагается постоянным), к = к^ -1) определено в предыдущем t-1 периоде. Задача заключается в определении числа фьючерсных контрактов к = к(г) на текущий момент времени. При этом предполагается, что на этот момент времени известна вся информация о спотовых и фьючерсных ценах на рассматриваемый актив. Управление риском позиции отождествляется с выбором количества фьючерсных контрактов к = к(г) на основе доступной текущей информации о ценах с целью уменьшения риска и получения (по возможности) дополнительного дохода. Ниже приводится схема адаптивного управления риском портфеля на основе прогнозов дисперсии изменений цен и ковариации между изменениями цен спотовых и фьючерсных цен.
1. На текущий момент времени получить прогнозы дисперсии изменений фьючерсной цены и уровень кова-риации Соу(ЧБ, у Я) на следующий момент времени t + 1.
2. Определить коэффициент хеджирования и оптимальное число фьючерсных контрактов согласно формулам (10).
Пример 3. Для оценки эффективности страхования использовались ежедневные данные по спотовым и фьючерсным ценам по акциям ОАО «Газпром» на интервале 11.01.2012-19.02.2012 г [7]. Исходная позиция содержит 1000 акций (Qs = 1000). На рис. 3 сплошной линией (ряд 1) - результат изменения стоимости портфеля (Qs = 1000, k) относительно начальной даты при адаптивном управлении риском. Штриховая линия (ряд 2) представляет незастрахованную позицию, то есть изменение стоимости 1000 акций. Из графика видно, что незастрахованная позиция к моменту t = 19.02.12 приведет к потере около 13 000 руб. Динамическое страхование приводит к потерям за весь период страхования не более 2000 руб.; при этом к концу рассматриваемого интервала доход составит около 800 руб. Из сравнения этих графиков видно, что динамическое страхование эффективно снижает риск больших потерь. Однако оно не приводит к существенному доходу.
Литература
1. Буренин А. М. Хеджирование фьючерсным контрактами фондовой биржи РТС. М.: НТО им. С. И. Вавилова, 2008.
2. Бокс Дж., Дженкинс Г. Анализ временных рядов: прогноз и управление. М.: Мир, 1974.
3. Керимов А. К. Методы анализа и прогнозирования ценовых данных. М.: РУДН, 2003.
4. Керимов А. К. Финансовые фьючерсные контракты. М.: Финансовый Университет, 2013.
5. Лукашин Ю. П. Адаптивные методы краткосрочного прогнозирования. М.: Статистика, 1981.
6. Халл Дж. К. Опционы, фьючерсы и другие финансовые инструменты. М: Вильямс, 2010.
7. Сайт фондовой биржи РТС. URL: http://www.rts.ru/ (дата обращения: 05.06.2015).