CC BY
299№2(22) 2011
Д. Фантаццини
Моделирование многомерных распределений с использованием копула-функций. I
Проблематика копула-функций, их свойств, способов подбора под конкретные исходные данные, оценивания, прикладных возможностей крайне скупо представлена в мировой специальной литературе и почти никак — в отечественной. При этом уже существуют впечатляющие примеры их прикладного использования в ситуациях, когда построение, статистическое оценивание и анализ многомерных распределений вероятностей оказывается необходимым инструментом исследования, а использование для таких задач модели многомерного нормального (гауссовского) закона не отражает специфики имеющихся исходных статистических данных. Есть основания утверждать, что модели, основанные на копула-функциях, окажутся особенно востребованными в прикладных эконометрических исследованиях, посвященных задачам оценки, анализа и управления финансовыми и страховыми рисками, доходностями разных финансовых инструментов. Предлагаемый в данном номере журнала материал является, по существу, фрагментом готовящегося к изданию учебника С. А. Айвазяна и Д. Фантаццини «Методы эконометрики. Продвинутый уровень».
Перевод на русский язык осуществлен А. В. Кудровым под научной редакцией С. А. Айвазяна.
Ключевые слова: многомерное распределение, эллиптические копула-функции, архимедовы копула-функции, иерархические копула-функции.
JEL classification: С69, С49.
Содержание
Введение
1. Копула-функции
1.1. Определения и основные свойства
1.2. Плотности копула-функций
1.3. Эмпирические приложения в статистическом пакете R: двумерные границы Фреше-Хёффдинга и копула-функция, отвечающая случаю независимости
2. Эллиптические копула-функции
2.1. Обзор
2.2. Нормальная копула-функция
2.3. Г копула-функция (копулаСтьюдента)
2.4. Численные примеры: оценивание эллиптических копула-функций
2.5. Эмпирические приложения в статистическом пакете R: эллиптические копула-функции
3. Архимедовы копула-функции
3.1. Основные понятия и определения
3.2. Некоторые однопараметрические двумерные архимедовы копула-функции
3.3. Многомерные перестановочные архимедовы копула-функции
3.4. Вложенные архимедовы копула-функции
3.5. Моделирование наблюдений из архимедовых копула-функций
3.6. Эмпирические приложения в статистическом пакетеR: архимедовы копула-функции Список литературы
Введение
Доказательства того, что не все совместные распределения экономических характеристик могут быть описаны при помощи многомерного нормального закона, широко представлены в различных работах, начиная с (Mills, 1927). Наиболее распространенными статистиками, позволяющими охарактеризовать степень отклонения от нормального распределения, являются коэффициенты эксцесса и асимметрии. Последние исследования показывают, что для финансовых данных также наблюдается отклонение их совместного распределения от многомерного нормального закона, а именно, отклонение в форме зависимостей компонент анализируемого многомерного признака. Один из таких примеров — асимметрия зависимостей доходностей рыночных активов, которая больше в периоды спада на рынке, чем в периоды рыночного роста, см. (Erb et al., 1994; Longin, Solnik, 2001; Ang, Chen, 2002; Patton, 2004, 2006a, б) и др.
Основной вывод этих статей заключается в том, что многомерное нормальное распределение не является хорошей моделью для описания совместного распределения многих экономических и финансовых переменных. Это приводит к проблеме поиска более адекватных многомерных моделей. Теория копула-функций — один из возможных способов ее решения.
Начало теории копула-функций было положено работами (Hoeffding, 1940) и (Sklar, 1959), но статистическое моделирование с использованием копула-функций в приложениях появилось относительно недавно и относится к концу 1990-х гг. Количество статистических исследований на тему копула-функций увеличивается достаточно быстро, см. (Wang, 1998; Frees, Valdez, 1998; Embrechts et al., 1999). Сегодня копула-функции успешно применяются в финансах и страховании, биостатистике (например, (Lambert, Vandenhende, 2002)), гидрологии (например, (Zhang, Singh, 2006)) и климатологии (например, (Salvadori, De Michele, 2007)). Среди учебников по теории копула-функций следует отметить (Joe, 1997) и (Nelsen, 1999), где представлено введение в теорию, тогда как в (Cherubini et al., 2004) и (Malevergne, Sornette, 2006) излагается использование техники копула-функций в финансовых приложениях.
Копула-функция является функцией, агрегирующей всю информацию относительно структуры зависимости между компонентами случайного вектора. Когда в качестве компонент копула-функции берутся частные функции распределения, которые необязательно принадлежат одному и тому же семейству распределений, получаем многомерную функцию распределения. Как следствие, эта теория позволяет достаточно гибко моделировать структуру зависимости между различными переменными, которые могут иметь разные частные распределения.
Несмотря на то что теория копула-моделей исследована относительно полно, проблема оценивания и статистические выводы для копула-моделей, в определенном контексте, все еще требуют дальнейших исследований, см. (Genest, Favre, 2007). За последние годы предложены различные методы оценивания параметров копула-функций, начиная с параметрических (Jondeau, Rockinger, 2003; Patton, 2004, 2006а, б; Fantazzini, 20096), полупараметрических (Genest et al., 1995; Breymann et al., 2003; Fantazzini, 2010), и заканчивая непараметрическими методами (Fermanian, Scaillet, 2003). Более того, недавние исследования направлены на исследование «смеси» указанных выше методов, которые позволили бы сэкономить время на вычислениях (Bouyé et al., 2000; Marshal, Zeevi, 2002; Cherubini et al., 2004).
№2(22) 2011
В данной статье будет представлен некоторый обзор теории копула-функций, основные методы оценивания, а также сделана попытка ответить на, вероятно, самый важный вопрос теории копула-функций: как следует подбирать корректную копула-функцию?
1. Копула-функции
1.1. Определения и основные свойства
С помощью копула-функций описываются законы многомерного распределения вероятностей. Они определяются частными одномерными распределениями анализируемого многомерного закона и характером зависимостей, существующих между компонентами рассматриваемой многомерной случайной величины.
Определение 1. Функция С(м1,и2,...,ип) от п переменных, определенная на единичном гиперкубе 1п = [0,1]п (т. е. м; Е [0,1], / = 1,2,...,п), называется копула-функцией, если она обладает следующими свойствами:
1) область значений функции — единичный интервал [0,1];
2) если и = 0 покрайнеймередляодного / Е {1,2,...,п], то С(и1,и2,...,ип) = 0;
3) С(1,...,1,и;,1,...,1) = и длялюбыхи; Е[0,1];
4) С(и1,и2,...,ип) является п-возрастающей функцией в том смысле, что для всех (а1,...,ап), | (Ь1 ,...,Ьп) Е [0,1]п с а1 < Ъ1 справедливонеравенство:
I §
| ^ =1
!
о где ип = а. и и.2 = Ъ. длявсех у Е {1,___,п].
м
(-1)^С (иципп) > О,
§ Например, если п = 2, четвертое свойство принимает следующий вид: длявсех (а1, а2), (61, Ь2) Е [0,1]2 с а1 < Ъ1 имеем,
0 >8
1
ш
§ 2 2 5 Wí_1Vl+i:
<з
та а.
л
! 4=1
с(ич,и2Ч)>о
¡1=1 ¡2 =1 2
О 2 [(-1)4+1 С(их.,и21) + (-1)>+2С(иц,и22) > о
§ О (-1)2С(ип,и21) + (-1)3С(ип,и22) + (-1)3С(и12,и21) + (-1)4С(«12,и22) > о
§ ^ С(и11, и21) — С(и11, и22) — С(и12, и21) + С(и12, и22) > 0,
II
та о о
^ где и}1 = а у и м. 2 = 6. при у Е {1,2] . Таким образом, свойство 4 равносильно тому, что
С(а1, а2) — С(а1, Ъ2) — С(Ь1, а2) + С(61, Ъ2) > 0 . (1)
100 =
Консультации •
№2(22) 2011
Следовательно, (1) можно интерпретировать как: |
§
Й
Р(0 < и, < а,, 0 < м2 < а2) - Р(0 < м, < а,, 0 < м2 < b2) - ig
в
-P(0 < u, < ¿j, О < м2 < а2) + P(0 < u, < b,, О < и2 < b2) > 0. c¿
Теперь сформулируем теорему Склара, которая показывает роль копула-функций в описании многомерных распределений вероятностей.
Теорема 1 (теорема Склара). Пусть H(•) — n-мерная функция распределения с частными распределениями F,,...,Fn. Тогда существует n-мерная копула-функция C(•) такая, что для всех действительных x,,...,xn:
H(xl,..., хи) = C(FX (xj),..., Fn (xn)). (2)
£слм все частные функции распределения непрерывны, то копула-функция определена единственным образом; в противном случае C(•) определена единственным образом лишь на области определения Ran F1 X Ran F2 X.. .X Ran Fn, где Ran — область значений частных функций распределения. Обратно, если C (•) — копула-функция, a F1,..., Fn — функции распределения, то функция H(•), определяемая выражением (2), является совместной функциейраспределения с частнымираспределениями FJ,..., Fn.
Доказательство. См. (Sklar, 1959, 1996; Joe, 1997;Nelsen, 1999).
Последнее утверждение теоремы представляет большой интерес в задаче моделирования многомерных функций распределения, поскольку из него следует, что можно связывать вместе любые n (n > 2) одномерных функций распределения разного типа (не обязательно из одного семейства), используя любую копула-функцию, для того чтобы получить двумерные или многомерные функции распределения.
Следствие. Пусть F1 (_1)(•),..1 (•) — обратные (в обобщенном смысле) функции частных распределений. Тогда для каждого (u1,...,un) из единичного n-мерного куба существует единственная копула-функция C : [0,1] X.. .X [0,1] ^ [0,1] такая, что
C u,..., un) = H (F )(u1 ),..., FT \un)). (3)
Доказательство. См. теорему 2.10.9 в (Nelsen, 1999).
Таким образом, копула-функция — это такая функция, которая, с использованием знания об одномерных частных распределениях, позволяет получить многомерную функцию распределения, поскольку функция распределения случайного вектора полностью описывает его вероятностную структуру, куда, в частности, входит структура зависимости его компонент. Копула-функции дают возможность разделить описание распределения случайного вектора на две части: частные распределения компонент и структура их зависимостей.
№2(22) 2011
Перечислим кратко варианты использования результатов теоремы Склара и рассмотренного следствия:
• извлечение копула-функций из хорошо известных многомерных распределений, например, из многомерного нормального распределения, из многомерного распределения Стьюдента ит.д.;
• конструирование новых многомерных распределений, связывая произвольные одномерные распределения при помощи некоторой копула-функции;
• в то время как соотношение (2) обычно служит отправной точкой в задаче моделирования, соотношение (3) является теоретическим инструментом для получения копула-функ-ции из многомерной функции распределения.
Следует отметить, что для дискретных распределений копула-функции имеют не столь простую форму, как в случае непрерывных распределений. Покажем это на следующем примере.
Пример 1. Рассмотрим две одинаковые кости, на которых выпадают два числа Х1 иХ2 из множества {1,2,...,6}. Предположим, что мы знаем значение Х1 и можем делать ставки на значения Х2, которые выпадают на второй кости. Важным вопросом является следующий: существует ли взаимосвязь или зависимость между этими двумя случайными величинами, т. е. возможно ли получить некоторую информацию относительно значений Х2, зная Х1. Мы знаем, что каждая случайная величина полностью опи-'! сывается своей функцией распределения ^ (х) = Р{Х1 < х). В нашем случае мы имеем | ^ (•) = ^ (•) = ^(•).
■&■ Однако в общем случае одномерная функция распределения не дает нам никакой инфор-Ц мации относительно совместного распределения двумерной случайной величины (Х1,Х2). | Правда, в нашем случае предполагается независимость ( Х1,Х2), поэтому совместное рас-§ пределение (Х1,Х2) имеет вид:
Р(Х1 < х1, Х2 < х2) = ^(х1) • ^(х2).
§ Для того чтобы получить полное описание (Х1,Х2), мы использовали два элемента: чао стные распределения и тип взаимосвязи (в нашем случае — независимость). Возможность | разделения между частными распределениями и зависимостью — это именно то, что ут-| верждает теорема Склара. Поскольку мы рассматриваем дискретные частные распределе-5 ния, то копула-функция будет определена единственным образом только на Яап ^ ХЯап .
§ В случае независимости имеем:
та
а
! С(т' т\ = 7 '7"' = 1>---,6 •
(I
и
0
1 §
6 6! 6 6
Любая копула-функция, удовлетворяющая этому ограничению, является подходящей: ¡5 например, С (и, у) = и -V будет естественным выбором.
^ Важное свойство, которым обладают все копула-функции — это инвариантность отно-| сительно строго возрастающих преобразований случайных величин (см. приведенную ни-| же теорему 2).
102 у~ =
Консультации •
№2(22) 2011
Теорема 2. Рассмотрим п случайных величин X1,...,Xn, зависимость между кото- | рыми определяется копула-функцией C(•). Если преобразования Tt : R ^ R, i = 1,...,n g определяются строго возрастающими функциями (т. е. dxTt > 0), то структура зависимости случайных величин T1(X1 ),...,Tn(Xn) определяется той же самой копула-функ- в цией C(•). 4
Доказательство. См. (Schweizer, Wolff, 1976, 1981), для двумерного случая — (Male-vergne, Sornette, 2006).
Последнее утверждение означает, что строго возрастающие преобразования не меняют структуру зависимости. Это значит, что при заданных частных распределениях сохраняется структура зависимости, которая определяется копула-функцией.
Более того, можно показать, что любая копула-функция ограничена так называемыми
границами Фреше-Хёффдинга:
w(üi,...,Un) < C(Mi,...,Un) < M(иг,...,Un),
где
W(u1,...,un) = max(u1 + ...+ un -n + 1,0), M(u1,...,un) = min(u1,...,un).
В качестве примера рассмотрим две равномерно распределенные на [0,1] случайные величины U1 и U2. Если U1 =U2, то эти две случайные величины имеют копула-функции вида
C (u1, u2) = P (U1 < u1 ,U1 < м2) = min(u1, u2)
и называются абсолютно зависимыми.
Ту же самую копула-функцию можно получить, если взять X2 = T(X1), где T(•) — монотонно возрастающее преобразование. Такие случайные величины X1 и X2 называются комонотонными. Прямо противоположным понятию комонотонности является противо-монотонностъ случайных величин. В качестве примера возьмем равномерные случайные величины (U1 ,U2), для которых U2 =1 — U1. Копула-функция случайного вектора (U1,U2) равна
C(u1, u2) = P(U1 < u1, 1 — U1 < u2) = P(U1 < u1, 1 — u2 <U1) = u1 + u2 — 1 при u1 + u2 >1.
В противном случае она равна нулю (рис. 1).
В двумерном случае границы Фреше-Хёффдинга являются копула-функциями, но для случаев большей размерности нижняя граница Фреше-Хёффдинга Wy:xe не является n-возрастающей (см. свойство 4 в Определении 1). Однако, неравенство в левой части нельзя улучшить, поскольку для любого фиксированного u из единичного n-куба существует C(•) такая, что W(u) = C(u), см. доказательство теоремы 2.10.12 в (Nelsen, 1999).
Копула-функции равномерно непрерывны и (почти всюду) дифференцируемы по каждому из аргументов — см. теоремы 2.2.4 и 2.2.7 в (Nelsen, 1999). Более детальная информация относительно этих и других свойств приведена в (Joe, 1997; Nelsen, 1999, 2006).
№2(22) 2011
Lower Frechet-Hoeffiding bound W(u,v) = max (u + v - 1,0)
Upper Frechet-Hoeffiding bound M(u,v) = min (u,v)
Product copula
P = u*v
й
>S
s
I
та
t
Si §
!
та о о
та
л §
С о X
0 >s
5
1 §
6
0 та а.
Й л
1
о. ф
§
о
0
1
5
!
та о
о
6
§ 5
Рис. 1. Графики коиула-фуикций Ж и М, а также коиула-фуикция Р, соответствующая
независимым случайным величинам
1.2. Плотности копула-функций
Для копула-функций, аналогично функциям распределения, можно определить понятие плотности. В частности, плотность е(и1,и2,...,ип), ассоциированная с копула-функцией С(и1,и2,...,ип), определяется соотношением
с(и1, u2,..., un) =
dnC (ui ,u2 ,...,un) du1 ...dun
Плотность может быть использована для того, чтобы определить так называемые абсолютно непрерывную АС и сингулярную БС компоненты копула-функции С(•):
u u2 un -¡n
AC (ui> u2 v, un ) = ff-f
dnC C?1, «2, Sn)
0 0 0 • • • ^Sn
ds1 ... dsn,
ИС (и1, и2,..., ип) — С (и1, и2,..., ип) АС (и1, и2,..., ип).
Копула-функция, для которой С(•) = АС (•) на 1п, называется абсолютно непрерывной. В противном случае копула-функция называется сингулярной.
Существуют копула-функции, которые имеют как абсолютно непрерывную, так и сингулярную компоненты. Например, копула-функция Р = и1и2и3 является абсолютно непрерывной, поскольку
104
u
u
u
№2(22) 2011
дС (^, ^, ^з ) |
— 1 s
ds1ds2ds3 Э
123 ¡5
идлялюбого (u1,u2,u3) Е13 имеем: ®
et
AC (u1,u2,u3) = fff-1' 2' 3 ds1ds2ds3 = u1u2u3> = P.
000 dSidS2dS3
Более того, применяя следствие теоремы Склара и рассматривая непрерывные случайные величины, можно видеть, что плотность копула-функции с(^(х1),...,Кп(хп)) ассоциирована с плотностью совместной функции распределения Н(•), обозначенной как /н (•), следующим образом (каноническое представление):
ЭпГс(К(хД...,Кп(х )) 1 -Нт гт
/н (х.-, хп)-п / (х)=С( к (х, к (хп ))-п у; (х),
эК1 (х1 )...зКп (хп) ^ ^
где
сК х),...,К (хпУ) = . (4)
П / (х;)
г=1
1.3. Эмпирические приложения в статистическом пакете Я: двумерные границы Фреше-Хёффдинга и копула-функция, отвечающая случаю независимости
Если необходимо изобразить линии уровня двумерных границ Фреше-Хёффдинга и копула-функции, соответствующей случаю независимости, то для этого можно использовать следующий R код (результаты представлены на рис. 2):
library (fCopulae)
# Generate Grid:
N = 50; uv = grid2d (x = (0:N)/N); u = uvx; v = uvy
# Compute Frechet and Product Copulae:
W = matrix (apply (cbind (u+v-1, 0), 1, max), ncol = N+l) Pi = matrix (u*v, ncol = N+l)
M = matrix (apply (cbind (u, v), 1, min), ncol = N+l)
# Create Perspective Plots:
persp (z=W, theta = -40, phi = 30, main = "Lower Frechet", cex = 0.5, ticktype = "detailed", col = "steelblue") persp (z = Pi, theta = -40, phi = 30, main = "Pi Copula", cex = 0.5, ticktype = "detailed", col = "steelblue") persp (z=M, theta = -40, phi = 30, main = "Upper Frechet", cex = 0.5, ticktype = "detailed", col = "steelblue")
№2(22) 2011
# Create Contour Plots: contour (W, xlab = "u", ylab = contour (Pi, xlab = "u", ylab = contour (M, xlab = "u", ylab =
"v", main = "Lower Frechet")
"v", main = "Pi Copula") "v", main = "Upper Frechet")
Рис. 2. Графики линий уровня двумерных границ Фреше-Хёффдинга и копула-функции,
отвечающей случаю независимости
>8
2. Эллиптические копула-функции
та
I
§ 2.1. Обзор
§
ф
I Класс эллиптических распределений включает, главным образом, класс симметричных
О распределений, который весьма популярен в актуарной математике, страховании и финан-
е; сах. Этот класс включает интересные примеры многомерных распределений, и многие
g из них имеют некоторые общие свойства с многомерным нормальным распределением. Они
о позволяют моделировать многомерные экстремальные события, формируя зависимость,
| не совпадающую с зависимостью многомерного нормального распределения, и использо-
| вать распределение с большим эксцессом, чем эксцесс нормального распределения. Кроме
5 того, с их помощью можно моделировать феномен «тяжелых хвостов», который часто на-
5 блюдается для финансовых данных, см. (Embrechts et al., 1999; Schmidt, 2002). £ Эллиптические копула-функции — это попросту копула-функции многомерных распре-
jg делений эллиптического типа. Если следовать работе (Fang et al., 1987), то эллиптические
g распределения определяются следующим образом.
ф
Определение 2 (эллиптические распределения). Пусть X есть «-мерный случайный век-
§
о и
0
1 §
| тор и 2Е МпХп — симметричная и неотрицательно определенная матрица. Если существует
щ /¿б Мп и функция у: М ^ М+ такие, что характеристическая функция X — /и имеет вид
! §
5
0X-, с )=exP0Y т м) и- 21т ^ I'
106
№2(22) 2011
• = \ln f g ( z )dz
В качестве примеров распределений эллиптического типа следует привести нормальное распределение, распределение Коши и распределение Стьюдента, логистическое распределение, распределение Лапласа, распределение Котца и экспоненциальное распределение, см. табл. 1.
Таблица 1. Список генерирующих функций g(t) и нормирующих констант у для некоторых копула-функций эллиптического типа
Распределение Генератор Константа
Нормальное е~"2 (2jt)4
Коши (l+2iT3/2 (2 л)"1
-v-2
2t\ 2 „
Стьюдента |1Н--1 (2 л)'
Логистическое е' (1 + е ') 2 л"1
Лапласа е(2л)'1
Котца е^лп г(2я)"'
Экспоненциальное е г(''2 ¡г11" (2яГ(1 / л))"
Отметим, что в некоторых работах используются иные определения генерирующих функций g(•), где ' / 2 заменено на ', что влияет на значение нормирующих констант у.
107
et
для любого t Е М", то X называют случайным вектором, имеющим распределение эллипти- |
ческого типа с параметрами 2. g Копула-функция C(•) называется эллиптической, если она отвечает распределению эл- £ липтического типа. Функция плотности распределения эллиптического типа, если она су- в ществует, имеет вид
/(х)=ViiT^(x2"1(хxем"
+ с00
для некоторой функции g : М ^ М , которая удовлетворяет условию: J о g(x)dx <ю. Функцию g(•) называют оператором плотности распределения эллиптического типа или генерирующей функцией, см. (Fang et al., 1987). Нормирующая константа у может быть определена явным образом с использованием перехода к полярным координатам, см. (Landsman, Valdez, 2003). В результате получим:
№2(22) 2011
Пример 2. Рассмотрим двумерный случай. Предположим, что X = (Х,У)т и
2 =
'1 р^ \Р 1 /
|2| = 1-р2, -1<р<1.
Тогда X имеет функцию плотности распределения
/Р (х, У) =
г-
р
^ х2 — 2 рху + у2 ^
(л, У) <
Следовательно, в соответствии с теоремой Склара, двумерная копула-функция эллиптического типа имеет вид:
Ре1 (и) Ре1 (и2) / 2 о ,2 \
Г( ■ л г V С I х -2 рху + У С (и1, «2 ;р) = г,—г J J
р
— го — го
1-р2
ёх dy.
(5)
Далее более детально рассмотрим две наиболее распространенные копула-функции эллиптического типа, а именно, нормального распределения и распределения Стьюдента.
'! 2.2. Нормальная копула-функция
Г
Копула-функция двумерного нормального распределения с коэффициентом корреляции Ц р может быть вычислена из уравнения (5) при g = е~"2 и у = 1 / (2ж):
§ 1 Ф1 и) Ф1 и)
ц С" С«1 , и2 ;р) =-/ / ехр
2 Жу/1 - р
^ х2 — 2 рху + у2 ^
— го — го
2(1-Р2)
ёх dy,
§
§ что может быть представлено в виде:
| С (и, и2 ;р) = Ф2 (Ф"1 и ),Ф-1 (и2 );р),
5 »
| где Ф2 (•) обозначает функцию двумерного нормального распределения с единичными дис-5 персиями компонент, корреляцией р и нулевым средним значением, а Ф"1 (•) —обратная
с
о
® персиями компонент, корреляцией р и нулевым средним значением, я ^ 1' функция для одномерного стандартного нормального распределения.
В общем случае многомерная нормальная копула-функция определяется следующим
а
ф §
о и
0
1 §
| где Фп (•,..., •; 2) — это п-мерная нормальная функция распределения с нулевыми средними
щ значениями и ковариационной матрицей 2 .
I
® Плотность нормальной копула-функции может быть получена из уравнения (4) с исполь-о
образом:
С Ц, и2 ,...,ип ;2) = Фп (Ф"1 (М1 ),Ф"1 (и 2 (ип );2)
зованием канонического представления: 108
№2(22) 2011
fG— (х,..„x ) (2л)-"121 2Г1'2 exp(-1x^"1 x) с(Ф(х1 ),...,Ф(X")) = J—-= V-^-1 =
о о I
П ¿""(xj П (2^)"1/2exp(-ix2) |
i=i ¡=1 . $
= exp(-2 (2"1 -I)?|, *
где £ = (Ф-1 (u),...,Ф"1 (un))т — вектор, компонентами которого являются значения обратной функции для стандартного одномерного гауссовского распределения в точках ut =Ф(х;), i = 1,...,п ; I — единичная матрица.
Представим основные шаги алгоритма моделирования нормальной копула-функции (более детальное описание моделирования копула-функции см. в (Cherubini et al., 2004; McNeil et al., 2005)):
• найти разложение Холецкого A для корреляционной матрицы 2 ;
• смоделировать из стандартного нормального распределения п независимых случайных величин z = (z1,..zn f;
• взять x = Az ;
• вычислить компоненты ut =ф(xt), i = 1,...,n, где ф(-) — одномерное стандартное нормальное распределение.
В результате получим вектор (u1,...,un)т, являющийся реализацией случайного вектора из n-мерной гауссовской копула-функции.
2.3. Г копула-функция (копула Стьюдента)
Копула-функция для двумерного распределения Стьюдента с корреляцией р может быть получена из формулы (5) при _ _
/ 2t Г^ 1
*Т+ VI • у=
Для Г копула-функции справедливо соотношение:
С^ (щ, ы2 ;р,у) = Т2 (С (щ), С1 (и2 );р,у),
где Т2—функцияраспределения двумерногораспределенияСтьюдентас корреляцией р иг степенями свободы; С^ (•) —обратная функция одномерного распределения Стьюдента. При г^<х> копула-функция Стьюдента сходится к нормальной копула-функции, а при у^ 1 сходится к копула-функции распределения Коши.
В общем случае многомерная Т копула-функция определяется следующим образом:
С (щ, Ы2,..., ип = Т" (Г1 (щ (и2),..., Г1 (ып ),
а плотность Т копула-функции может быть получена из уравнения (4) и канонического представления:
109
№2(22) 2011
р/ v+n \
c(tv(х),...,tv(хп)) = |2|"1/2 V 2 '
г(*)
П (i+|2)
y+n 2
v+1
2
где £ = (t„1 (u1),..., t„1 (un ))т — вектор, компонентами которого являются значения обратной функции распределения Стьюдента в точках ut = tv (xt).
Приведем алгоритм моделирования наблюдений, подчиняющихся распределению с Т копула-функцией:
• найти разложение Холецкого A матрицы 2;
• смоделировать n независимых случайных величин из стандартного нормального распределения, z = (z1zn )т;
• смоделировать случайную величину s из ^^-распределения, не зависящую от z ;
• определить вектор y = Az;
4v
• положить x = —j= y ;
<S
• определить компоненты ui = tv (xi), i = 1,..., n, где tv(-) —одномерное стандартное pac- пределение Стьюдента с v степенями свободы.
'! Вектор (u^...,un)т смоделирован из n-мерной Tкопула-функции.
^ Нормальная и T копула-функции весьма популярны при анализе финансовых данных, * поскольку они позволяют моделировать портфели высокой размерности, являясь при этом Ц относительно легко оцениваемыми. Кроме того, они позволяют генерировать наблюдения | при помощи относительно простых алгоритмов. На рисунке 3 представлены графики плот-§ ностей некоторых эллиптических копула-функций, построенных с использованием функ-I ции fcopulae из пакета R.
Щ
о
п
г? Обобщение: сгруппированные Г копула-функции. Сгруппированные T копула-функ-
g ции впервые были предложены в работах (Daul et al., 2003) и (Demarta, McNeil, 2005). Калиб-
o ровка таких копула-функций не сложнее, чем калибровка T копула-функций, но для первых
| имеется возможность определять различные структуры зависимости для подгрупп компо-
| нент вектора из T копула-функции.
<f Пусть Z Е Nn (0, R), где R — ковариационная матрица размерности n X n с единичны-
§ ми диагональными элементами. Далее рассмотрим случайную величину U, равномерно рас-
С5
£ пределенную на отрезке [0,1] и не зависящую от Z. Пусть через Gv(-) обозначена функция
Л распределения случайной величины^ v / %2V , где %2V — случайная величина, имеющая хи-о квадрат распределение с v степенями свободы. Кроме того, задано разбиение множества
о
s §
g
О §
5
1,___,п на т подмножеств с размерами ^,...,5т (т. е. ^,...,5т >1, ^ + ... + 5т = п) и для каждого к = 1,2,...,т задано число степеней свободы ук. Возьмем Жк = G~1 (Ц), к = 1,___,т,
| и У = (+1, ). Случайная величина У имеет так назы-
ваемое сгруппированное ¿-распределение. Наконец, определим
U = (t,1 (11 ),...,t„1 Y ),tVi (YSi+1),-,tVi (Y^ ),...,t% Yn)).
110 J
№2(22) 2011
Рис. 3. Плотности эллиптических копула-функций: нормальное распределение, распределение Стьюдента, логистическое и экспоненциальное распределения. Параметры указаны над соответствующими им графиками
Случайный вектор и имеет распределение на [0,1]", а его компоненты имеют равномерное распределение на [0,1]. Будем называть распределение вектора и сгруппированной Т копула-функцией.
Отметим, что (7р..., 75) имеет многомерное Т-распределение с г1 степенями свободы, идля к = 1,..., т-1 вектор (7^ + + +1 + ) имеет многомерное Т-распределение с ук+1 степенями свободы. Соответствующее распределение подвектора вектора и представляет собой Ткопула-функцию с Ук+1 степенями свободы, для к = 0,...,т — Для сгруппированной Т копула-функции нет простого выражения для плотности. Однако существует весьма полезная аппроксимация для автокорреляции (в предположении ее постоянства), полученная в работе (Оаи1 е1 а1., 2003):
Рц (^, ^ ) ~ 81П(ЖСд (щ, и] ) / 2),
где г и у принадлежат различным группам, а г^ — парный коэффициент ранговой корреляции Кендалла. Эта аппроксимация позволяет производить оценку параметров методом максимального правдоподобия независимо для каждой подгруппы. В работе (РаЩаггЫ,
№2(22) 2011
2009а) представлено обобщение модели из (Оаи1 е1 а1., 2003), в котором рассматривается динамическая структура корреляционной матрицы Ы, .
2.4. Численные примеры: оценивание эллиптических копула-функций
Нормальная копула-функция. Выше было показано, как, применяя теорему Склара и используя соотношение между функцией распределения и ее плотностью, можно получить плотность нормальной (гауссовской) копула-функции:
/Оаиаыап / \ л
^ ,..., хп ) 1 си, п2,..., ип = —п-= з;т ехР
| | ^ Саш$1ап ^ х ^
- 2 £ (2"1 -Щ
где £ = ( Ф 1 (и1),..., Ф 1 (ип ))т ,аи( =Ф( х;); логарифмическая функция правдоподобия имеет вид:
Iе
\в) = 4ш|2|-22С:(2"1 .
>8 I
та
I §
§
ф £
та о о
та
л §
С о X
0 >8
1
ф §
§■
0 та а.
Й л
1
а.
ф §
о и
0
1 §
ф
!
та о
о &
§ 5
Если логарифмическая функция правдоподобия дифференцируема по параметрам в, а решение уравнений дв /(в) = 0 дает точку глобального максимума функции /(в), можно
получить оценку максимального правдоподобия для ковариационной матрицы в ж =2 гауссовской копула-функции:
Т 1
2 2-
а, следовательно,
Э2
л Т 1 Т
¡Оа^ап (б) = _ _ ^ ?Т= ^
*=1
-1
1 Т
^=1, £'
Т 11
т. е. получаем классическую оценку для ковариационной матрицы. Следовательно, в качестве приближенной оценки корреляционной матрицы можно использовать матрицу:
Ы = (Шаё 2)_1/2 2 (Шаё 2)
4-1/2
Г копула-функция. Плотно стьГ копула-функции равна:
(6)
сЦ, и2,..., ип ;2):
П Л^ (х)
1 ) гй
И * г0 )
1+
" /
П(1+
у+1 ?
V
а логарифмическая функция правдоподобия имеет вид 112
!=1
,=1
г=1
г=1
№2(22) 2011
,"(в)=_П11 ГШ—nT in Щ-T ид-— yuL^)^ у Уш^Ш |
r(i) Г© 2 11 2 £ [ v j 2 ££ \ v J i[
<0 в
В случае T копула-функции невозможно получить аналитическое выражение для оценки 4 максимального правдоподобия, поэтому необходимо использовать численные алгоритмы максимизации функции правдоподобия. Однако в случае большой размерности T копула-функции такие алгоритмы являются слишком затратными по времени. В связи с этим предлагаются альтернативные алгоритмы, многошаговые параметрические или полупараметрические подходы. Рассмотрим три наиболее часто применяемых метода: первый из них предложен в работе (Bouye et al., 2000) и основан на процедуре рекурсивной оптимизации для корреляционной матрицы; второй представлен в работе (Marshal, Zeevi, 2002) и базируется на использовании оценки ранговой корреляции Кендалла; третий является смешанным параметрическим методом и основан на методах моментов и максимального правдоподобия (Chen et al., 2004).
Метод 1 (Bouye et al., 2000).
1) Используя эмпирические функции одномерных частных распределений, преобразуем {Xt>x2t v,xnt), t = 1,...,T, к приблизительно равномерно распределенным (и 1t,и2t,...,иnt).
2) Для каждого значения степени свободы v из рассматриваемого интервала значений оцениваем корреляционную матрицу Rv:
а) длякаждогофиксированного t (t -I, ...,T )положим:
tt = (CU ),..., C(unt ))T;
б) пусть 2 — корреляционная матрица для гауссовской копула-функции, оцениваемая по формуле (6), положим 2v,i = 2;
в) пусть 2v,*+1 удовлетворяет следующей рекурсивной схеме:
. T _
2 v^ + IV v)"1;
T t=1
г) перемасштабируем элементы матрицы следующим образом:
(2 „¿+1).;.
(2 v,k+\)i
в результате получим, в частности, единичные диагональные элементы;
д) повторяем шаги в)-г) до тех пор, пока 2^+1 = 2^, и положим 2„ = 2•
3) Оценим V путем максимизации логарифмической функции правдоподобия плотности Т копула-функции:
№2(22) 2011
T
V = arg max V ln cstudent (U u,..., U „ ,t; 2v, v).
v
t=l
Метод 2 (Marshall, Zeevi, 2002).
Предыдущая процедура является вычислительно трудоемкой в случае большой обучающей выборки. Более того, она может давать неустойчивые результаты в ситуациях, когда существует к (к > 1), для которого матрица 2v,k близка к вырожденной. Для преодоления этой проблемы в работе (Marshall, Zeevi, 2002) предлагается следующий алгоритм:
1) Используя эмпирические функции одномерных частных распределений, преобразуем (x1t,x2t,...,xnt), t = 1,_,Гк приблизительно равномерно распределенным (и 1t,и2t,...,иnt).
2) Оценим 2 с использованием непараметрической оценки рангового коэффициента Кендалла (Lindskog et al., 2002):
2j = sin(f ? j
i, j = 1,..., n.
Как указано в работах (Lindskog, 2000; Genest et al., 1995), если имеется два временных ряда Xt и Yt, t =!,...,T, то состоятельнойоценкойкоэффициентаКендаллабудет:
>s s
I 2 ^
I = T(T_1)2S1gn[(^i -)Y-Yj-)]' i j = 1'-'n,
t 'j T(T-1) tj j j IL r
с [ 1, änee x > 0
§ äää sign(x) = ]
§ 1-1, änee x < 0.
Ф L '
!
|5 3) Оценим v путем максимизации логарифмической функции правдоподобия плотно-| сти T копула-функции:
T
£ = argmax J lncStudent (Uu,...,Unt;2;v).
л
is
с u s
u
1 t=1 t Ф
is
¡F Метод 3 (Chen et al., 2004).
§ 1) Используя эмпирические или параметрически оцененные функции одномерных ча-
С5
£ стных распределений, преобразуем (x1t,x2t,...,xnt), t = 1,...,T к приблизительно равномерно
jg распределенным (и 1t,и2t,...,иnt).
® 2) Пусть R — корреляционная матрица^для гауссовской копула-функции, оцененной о
0
1 S
|| 3) Оценим v путем максимизации логарифмической функции правдоподобия плотно-щ сти T копула-функции:
I г
ч T
§ v = argmax ^ ln cSt"dent (и 1t,..., и nt ^oamsian ;v).
с использованием уравнения (6). Положим 2 Gaussian = R
t=1
114
I №
2(22) 2011
4) Пусть £ = (tv (u\t),...,i (unt)) .Наконец,используяуравнение(6),вычислим: g
!
1 J ¡5
^ = 1 у г ^т -с
Student 1 ¿J vt^> W ф
1 t=l cj.
Здесь также могла быть применена итеративная процедура, но авторы работы (Chen et al., 2004) не делали этого. Однако после первого шага разница между фактической ковариационной матрицей и ее оценкой весьма мала.
2.5. Эмпирические приложения в статистическом пакете R\ эллиптические копула-функции
В модуль «АСори1ае» входит множество интересных функций, позволяющих исследовать эллиптические функции. Приведем некоторые команды:
fcopulae функция Описание
ellipticalList список доступных эллиптических копула-функций
ellipticalParam установка параметров по умолчанию
ellipticalRange диапазон допустимых значений коэффициента корреляции
ellipticalCheck проверка принадлежности коэффициента корреляции допустимому диапазону
Для функции распределения и плотностей имеются следующие команды:
fcopulae функция Описание
rellipticalCopula моделирование наблюдений из эллиптической копула-функции
pellipticalCopula вычисление значения эллиптической копула-функции
dellipticalCopula вычисление значения плотности эллиптической копула-функции
rellipticalSlider график случайных величин, смоделированных из эллиптической копула-функции
pellipticalSlider график функции распределения
dellipticalSlider график функции плотности
Если необходимо изобразить график функции плотности и линий уровня нормальной копула-функции, нужно использовать следующие команды (результат см. на рис. 4):
# Generate the 2D Grid:
N = 50; х = (0:N)/N; uv = grid2d (x); u = uv$x; v = uv$y
# Compute the Normal Copula Density:
c. uv = dellipticalCopula (u, v, rho = 3/4, type = "norm", output = "list")
# Create a Perspective Plot:
persp (c. uv, theta = -40, phi = 30, ticktype = "detailed", col = "steelblue", main = "Normal Copula Density", cex = 0.5)
# Create a Contour Plot:
contour (c. uv, nlevels = 20, main = "Normal Copula Density", cex = 0.5)
№2(22) 2011
Elliptical Copula Density No: 1 - Normal rho = 0.5
Normal Copula Density
Рис. 4. График плотности нормальной копула-функции и ее линий уровня
Для представления о графике эллиптической копула-функции и ее плотности при различном выборе параметров рекомендуем воспользоваться «слайдером» функций:
# Start Slider for Perspective Distribution Plots: — pellipticalSlider () 5 llcm!
3
t # Start Slider for Perspective Density Plots: dellipticalSlider ()
та
I
§ В случае, если необходимо оценить эллиптические копула-функции и сгенерировать coll ответствующие им наблюдения, можно использовать следующий код:
!
та о о
та
¡о # Estimation of 100 Samples:
С est = NULL о s u >s
5
® ans = ellipticalCopulaFit (R, type = "t" ® *
4 ans=ans$par if
6 est = rbind (est, ans)
0
ü # Print the Result:
л
1
a.
<u
§ ...
о densl=density (est [,1]) plot (densl,main="kernel density of estimated Rho" о ...
* dens2=density (est [,2]) plot (dens2,main="kernel density of estimated NU")
!
та о
о &
s
5
for (i in 1:100)
R = ellipticalCopulaSim (n = 100, rho = 0.6, param = 4, type = "t")
for (i in 1:2) print (c ( mean = mean (est [, i]), sd=sd (est [, i] # Make a kernel density of the 100 estimates
116
№2(22) 2011
3. Архимедовы копула-функции
3.1. Основные понятия и определения
та
е
Архимедовы копула-функции обеспечивают аналитическую гибкость и широкий спектр различных мер зависимости. По следующим причинам эти копула-функции могут быть использованы в широком диапазоне приложений:
• Архимедовы копула-функции могут быть представлены в явном аналитическом виде, в отличие от семейства эллиптических копула-функций, которые определяются в неявной форме.
• Архимедовы копула-функции допускают относительно простое построение, включая вычислителънуюреализацию (см. табл. 2).
• Многие параметрические семейства копула-функций принадлежат этому классу.
• Архимедовы копула-функции не ограничены обязательным наличием радиальной симметрии, что свойственно случаю нормальной копула-функции или общему случаю эллиптических копула-функций. Это является преимуществом, т. к. во многих финансовых и страховых приложениях наблюдается более сильная зависимость между большими убытками, чем между большими доходами.
Двумерные архимедовы копула-функции могут быть определены следующим образом.
Определение 3. Рассмотрим непрерывную, строго убывающую и выпуклую функцию ср(и) с неотрицательными значениями, определенную при и Е [0,1] и удовлетворяющую условию <р(1) = О.Определимпсевдообратнуюфункцию <р[_1](г) соотношением:
называют архимедовой копула-функцией с генератором <р(-). Более того, если <р(0) = « , то псевдообращение дает обычную обратную функцию (т. е. <р[_1] (•) = ср~1 (•)), а ср(-) и С(•) в этом случае называют строгим генератором и строгой архимедовой копула-функцией, соответственно (в противном случае они называются нестрогими), см. рис. 5.
В работе (Кекеп, 1999) представлен список из 22 однопараметрических двумерных архимедовых копула-функций. Этот список вместе с генерирующими функциями ср(-) и соответствующими диапазонами параметра а воспроизведен в табл. 2.
Копула-функции с номерами 3-6, 9, 10, 12-14, 17, 19, 20, 22 являются строгими копу-ла-функциями, тогда как копула-функции с номерами 1 и 16 — строгие только для а> 0 . Остальные копула-функции из табл. 2 не являются строгими. Некоторые из архимедовых копула-функций имеют специальные названия. Например, копула-функцию под номером 1 называют копула-функцией Клейтона, под номером 4 — копула-функцией Гумбеля,
) =
<р 1 (г) для О < г < <р(0); 0 для г > ср(0),
где (р 1 (г) — обычная обратная функция к функции ср(и). Тогда функцию С: [0,1]2 ^ [0,1], определенную как:
С (и, и2) = ср1 1][^(и1) + <ри)],
(7)
№2(22) 2011
Strict generator
Non strict generator
1- \ <p
Inverse
Inverse
>S S
I
та
t
Si §
ф
!
та о о
та
л §
С о X
0 >s
5
1
ф
§
6
0 та а.
Й л
1
а
<u §
о
0
1 §
ф
!
та о
о &
§ 5
Рис. 5. Строгие и нестрогие генераторы и их обращения
а под номером 5 — копула-функцией Франка. Учитывая, что последние три копула-функ-ции довольно популярны, они будут рассмотрены ниже более подробно.
Таблица 2. Список архимедовых копула-функций
Номер
Ca(u,v)
Генератор <ра(() Диапазон параметра а
1 2
3
4
5
6
7
8
118
max([u~a + v^-lf1'", 0) max(1 - [(1 - u)a + (1 - v)a ]1M, 0)
1-«(1- u )(1- v) exp(-[(-ln u)a + (-ln v)a ]1M)
_1hL ^-1)
« { e~a-1 ,
1 - [(1 - u f + (1 - v)a - (1 - u У (1 - v)a fa
max(auv + (1 — a)(u + v — 1), 0) a2uv — (1 — u)(1— v)
ln
(ra-1)/ a
(1-t)" 1-a(1-1)
t
(-In t)a
, e-at -1 -ln-
a2 - (a-1)2 (1- u)(1- v)
,0
e --1 -ln[1- (1-1)" ] -ln[ at + (1-a)] 1-1
1 + (a-1)t
[-!.») 1 {0} [1,»)
[-1,1)
[1,»)
\ {0}
[1,») (0,1] [1,»)
uv
№2(22) 2011
Окончание табл. 2
I
та
е Ct
Номер
Ca(u,v)
Генератор <pjt) Диапазон параметра а
9
10 11 12
13
14
15
16
17
18
19
20
21 22
uv ехр(—а ln u ln v) uv[1 + (1- u a )(1- va )Г1/а max([uavа - 2(1 - u)a (1 - v)a]1M,0)
(1 + [(u^1 -1) a + (v1 -1)а ]1/а )_1 exp(1 - [(1 - lnu)a + (1 - ln v)a - 1]1M) (1 + [(u~1/a! -1)а + (v"1/a -1)а]1Mya max ({1 - [(1 - ulla )a + (1 - vl,a )a fa}a, о)
2(S + VS2 + 4«), S = u + v-1-«(u"1 + v1 -1)
1 +
[(1 + ur-1][(1 + v)"a-1]
max (l + а/ ln[e
1
*'("-1) + ga '(v-1)
1
]. 0)
a /In(ea'u + ea'" - ea) [ln(exp(u~a) + exp(v~a) - e)]"1'а
1- (1-{ max([1- (1- u)a fa +
+[1- (1- v)a ]1M-1,0)}a fa
max (|~1- (1- uа1- (1- uа)2 -
/-]1'« \
-(1- v * W1- (1- v" )2 J , oj
ln(1 — aln t) ln(2ra-1) ln(2 —ta) (Г1 -1)" (1- In t) a-1
(r1M-rf (1- tlla у
(«Г1 +1)(1-1)
-ln t1^-1 2~a —1
e«/(u-1)
ea" - ea ехр(Га) - e
1-[1- (1-1)a ]1M
arcsin(1 — ta)
(0,1] (0,1] (0, 1 ] [1,»)
[1,») [1,»)
[o,«o \ {0}
[2,») [o,«o (0,~)
[1,»)
(0,1]
Плотность архимедовых копула-функций может быть вычислена по формуле:
c(u1, u2) =
d2C (u1, u2) du1du2
или, если генератор <p(-) дважды непрерывно дифференцируем,
С (u1, u2 ) = ■
(^УЦ) + <P(u2 )])"
Это означает, что, зная генератор ср(-), его обратную функцию <р[_1] (•), ее первую и вторую производные, можно получить архимедову копула-функцию и ее плотность. Однако, за исключением некоторых простых случаев, плотности копула-функций оцениваются с использованием таких программ, как Maple или Mathematica.
119
№2(22) 2011
3.2. Некоторые однопараметрические двумерные архимедовы копула-функции
Копула-функция Клейтона. Семейство копула-функций Клейтона было впервые предложено в работе (Clayton, 1978). Рассмотрим генератор <p(t ) = (t~a -1) la с а£[-1,») \ {0}. Тогда <р~1 (t) = (1 +1)_1/а . Используя (7), получим:
\ — 1/а
Cclay,on (щ,u2;а) = max (u"a + u~а-1)" 0
Если а> 0, то <р(0) = °° и вышеприведенное выражение для копула-функции примет вид: _уа
С (м1, и2) = (и^а + и2а — 1) ,
а плотность копула-функции Клейтона будет равна:
е(и1,и2) = (1 + «)(и1и2)"а_1 (и~а + у~а -1)~"~2.
Отметим, что в отличие от случая а> 0, при —!<а< 0 копула-функция Клейтона не яв-'! ляется строгой, поскольку ср(0) = 11а. При а^О получим копула-фунцию, отвечающую ^ случаю независимости, а при а^ю — копула-функцию, соответствующую комонотонно-сти. При а = — 1 получим нижнюю границу Фреше-Хёффдинга. Таким образом, копула-Ц функция позволяет интерполировать промежуточные структуры зависимости между про-
| тивомонотонностью, комонотонностью и независимостью. §
ф
| Копула-функция Гумбеля. Семейство копула-функций Гумбеля было впервые представ-
|5 лено в работе (ОитЬе1,1960), но поскольку подробный анализ этого семейства был приведен
| в работе (Но^аагё, 1986), то его иногда называют семейством Гумбеля-Хоугарда. Генератор
о (р 1 (t) = ехр(—t1/a ). Таким образом, копула-функции из этого семейства имеют вид:
g этого семейства имеет вид: <p(t) = (—In t)а, где а>1. Обратная функция генератора равна
IM образом, копула-функции из этого ceiv
C (u1, u2 ;а) = expl— [ (— ln u1)a + (— ln u2)a]
-i1/a|
с функцией плотности:
>8
I
Ф §
§■
о та а.
Й 2«/
5 е(и1, и2 ;а) = С (и1, и2) •и~1и21| (— 1п и1 )а + (— 1п и2) а\ х
§
о и
0
1 §
| При а = 1 получаем копула-функцию, соответствующую случаю независимости; а при
щ а^со копула-функция Гумбеля С(и1,и2;а) стремится к копула-функции, отвечающей ко-
^ монотонности. Таким образом, копула-функция Гумбеля интерполирует структуру зависи-
® мости, промежуточную между независимостью и абсолютной положительной зависимо-
§ стью.
120
х [ln u1ln u2 f1 - {[(-In u1 )a + (-In u2 у ]"1/a + a-1}.
J
№2(22) 2011
Копула-функция Франка. Впервые копула-функция Франка была представлена в работе (Frank, 1979), а более подробный ее анализ проведен в работе (Genest, 1987). Эта копула-функция имеет генератор:
<p(t) = ln
1
-1
I
та
e et
Обратная функция этого генератора равна:
(р~1 (') = - -1п(1 + е' (е~
а
а сама копула-функция и ее плотность имеют вид:
■1)),
C (м1, u2 ;а) =--ln
а
c(u1, u2 ;а) = -
1 ! ^—au1 л\ / —au~ л \ \
1 +
{e~aUl -1)-(e~au2 -1)
e~a — 1
-a[e~a-1]e "a(U1+U2}
{[e"aU1 -1][e"aU2 -1] + e~a — 1}2
для аФ 0.
При а = 0 копула-функция Франка имеет вид копула-функции, отвечающей случаю независимости; а при а^ — <» и а^ + ю копула-функция Франка принимает вид, соответственно, нижней и верхней границ Фреше-Хёффдинга.
Для того чтобы продемонстрировать различия между копула-функциями Клейтона, Гум-беля и Франка, на рис. 6 они представлены графически. Видно, что изображенные копула-функции имеют различный характер поведения в окрестностях точек (0,0) и (1,1). Например, копула-функция Гумбеля в окрестности (1,1) демонстрирует резкий рост, тогда как рост в окрестности точки (0,0) менее резок. Как будет показано далее, хорошим измерителем характера хвостов этих копула-функций является коэффициент «хвостовой зависимости». Копула-функция Гумбеля имеет «верхнюю хвостовую зависимость», копула-функция Клейтона — «нижнюю хвостовую зависимость», а копула-функция Франка не имеет хвостовой зависимости.
Archimedian Copula No: 1 - Clayton Archimedian Copula No: 4 - Gumbel-Hougard Archimedian Copula No: 5 - Frank
[-1|Inf] alpha = 2 tau = 0.5 rho = 0.582 [1|Inf] alpha = 2 tau = 0.5 rho = 0.682 [-Inf | Inf] alpha = 5 tau = 0.451 rho = 0.634
Рис. 6. Плотности копула-функций Клейтона, Гумбеля и Франка
121
e
№2(22) 2011
3.3. Многомерные перестановочные архимедовы копула-функции
В работе (Nelsen, 2006) показано, что «-мерную копула-функцию для случая независимости можно записать в следующем виде:
P(Wj,u2,...,un) = uj • u2 •...•un = exp(-[(-lnuj) + (-lnu2) + ...+ (-lnun)]).
Это наблюдение подсказывает возможность обобщения (7) на n-мерный случай:
C (uj, u2,..un) = J] [^(uj) + (p{u2) + ...+ (p(un)]. (7a)
Функции C(•) в (7a) называют серийно-итеративными функциями (Schweizer, Sklar, 1983), основанными на двумерной архимедовой копула-функции с генератором ср(-). Если положить C(uj,u2) = <p[_j][<p(ul) + <p(u2)],тодля n >3 получимследующееобобщение:
C (uj, u2,..., un ) = C (C (uj, u2un-j ), un )•
Двумерные архимедовы копула-функции обладают свойствами симметричности и ассоциативности (см. теорему 4.1.5в (Nelsen, 2006)), а именно:
'SS
S
Ig C (u, v) = C(v, u) Vu, v e [0,1] (neüäöSe^fTnöü),
C(C(u,v),w) = C(u,C(v,w)) Vu,v, w E [0,1] (ännTöeäöeäfTnöü).
§ Копула-функции вида (7а) также известны как перестановочные архимедовы копу-
§ ла-функции: случайные величины X и У называются перестановочными, если векто-
| ры (X,У) и (У,X) являются одинаково распределенными. В нашем случае перестано-
(5 вочность копула-функции эквивалентна ее симметричности (см. теорему 2.7.4 в (№1зеп,
§ 2006)).
§ Для того чтобы С(•) была копула-функцией при п > 3, требуются некоторые допол-
о нительные свойства функций ср(-) и <р[_1] (•). Для этого понадобится следующее опреде-
| ление. §
5 Определение 4. Неотрицательная функция ё(?) называется абсолютно монотонной
§ на интервале если она непрерывна на этом интервале и имеет производные любого по-
£ рядка, причем
о и
0
1 §
и
I
щ Например, функция g(?) = е *, t Е Я является абсолютно монотонной, т. к. ) = —е*, ё"(() = ,...,г>(*) = (-1)*е"г,... .
для всех внутренних точек t Е J и к = 1,2,....
о
® Следует отметить, что если ё^) абсолютно монотонна на [0,») и ё(с) = 0 для некоторо-§ го с > 0 , то ё(•) = 0 на [0,»). Более того, если псевдообратная функция <р[_1] (•) генератора
122 у =
Консультации •
№2(22) 2011
архимедовой коиула-фуикции <р(-) абсолютно монотонна, то на интервале [0,») функция |
<р[_1] (•) положительна,функция ср(-) абсолютномонотоннаи <р[_1] (•) = ср~1 (•). g
Ниже приведем необходимое и достаточное условие, налагаемое на строгий генератор £
ср(-), для того чтобы функция (7а) была «-мерной копула-функцией при n > 2 (Kimberling, в
1974). ^
Теорема 3. Пусть <р\ [0,1] ^ [0,ю] — непрерывная и строго убывающая функция, такая, что <р(0) = ю и <р( 1) = 0, а ср~1 (•) — обратная функция для (р(-), C : [0,l]n ^ [0,1] определяется соотношением (7а). Тогда необходимым и достаточным условием для того, чтобы C (•) была n-мерной копула-функцией для всех n > 2, является требование абсолютной мо-нотонностифункции ср~1 (•) на [0,»).
Доказательство. См. (Kimberling, 1974; Schweizer, Sklar, 1983; Aisina et al., 2005).
Однако требование абсолютной монотонности в вышеприведенном условии является весьма жестким, что приводит к ограничениям структуры зависимости. В работах (Genest, Rivest, 1993; Nelsen, 2005; Müller, Scarsini, 2005) было предложено заменить это условие более слабым, в котором для функции ср(-) требуется существование производных лишь до некоторого конечного порядка. Недавно в (McNeil, Neslehova, 2009) представлено необходимое и достаточное условие для того, чтобы функция ср(-) обладала свойством «-монотонности (^-монотонности в терминах этой работы). Это условие показывает, что существуют n-мерные архимедовы копула-функции, не имеющие плотностей. Кроме того, из результатов работы следует существование точных нижних границ для множества всех n-мерных архимедовых копула-функций относительно так называемого конкордационного упорядочивания (для более детальной информации см. (McNeil, Neslehova, 2009)).
Представленное выше многомерное обобщение архимедовых копула-функций является весьма ограничительным, поскольку в своей спецификации оно использует лишь один генератор, не зависящий от рассматриваемой размерности. Как следствие, все ^-мерные частные распределения ( k < п) имеют идентичный вид. Для того чтобы обеспечить большую гибкость, было предложено множество подходов, некоторые из которых мы рассмотрим ниже.
3.4. Вложенные архимедовы копула-функции
Построение полностью вложенных архимедовых копула-функций. Описание и исследование одного из обобщений для функций (7а) имеется в работах (Joe, 1997; Embrechts et al., 2003; Whelan, 2004; Savu, Trede, 2009; McNeil, 2008). Структура такого обобщения для четырехмерного случая показана на рис. 7.
Узлы u1 и u2 связываются при помощи копула-функции C1 (•), узел u3 связывается с C1 (u1,u2) при помощи копула-функции C2 (•) и, наконец, узел u4 связывается с C2 u 2 C1 (u j u2)) при помощи копула-функции C3 (•). Таким образом, четырехмерный случай требует трех двумерных копула-функций C1 (•), C2 (•) и C3 (•), генераторы которых равны ср1 (•), ср2(•) и <р3(•) соответственно. Несмотря на то, что структура «связывания» довольно проста, ее выражение в виде аналитического уравнения является весьма непростой задачей:
№2(22) 2011
Рис. 7. Построение полностью вложенных архимедовых копула-функций
С(ихип) = <p\i (<pn_i о<pj_2[...(<р2 о^1 Ц) + ^(м2)] + (р2(u3))
+... + <?2 (Un_i )] + ^n-i (Un)) ■ (8)
С(^1 1 ) ? ^ — С3 4' С2 (^3 ? С^ (^1 ? ))) — = ^з1 {^3 (и4 ) + <рз (<р21 {^2 (иэ ) + <Р2 (^Г1 (и1) + 01 (и2 )})})}•
- Другими словами, пары (и1, и3) и (и2, и3) имеют одну и ту же копула-функцию С2 (•) спа-'! раметром зависимости а2, тогда как пары (и1,и4), (и2,и4) и (и3,и4) имеют копула-функ-^ цию С3 (•) спараметромзависимости а3 .Вобщем п-мерномслучае:
та
I §
§ £
та о о
та
¡5 Эта копула-функция называется полностью вложенной архимедовой копула-функцией,
§ поскольку получается путем пошагового увеличения размерности на единицу. Такая струк-
о тура построения копула-функций является более общей, чем (7а). Кроме того, она обладает
| лишь частичной перестановочностью. Полностью вложенная архимедова копула-функция
| имеет п(п — 1) / 2 различных двумерных частных распределений, но только п — 1 несовпа-5 дающих двумерных копула-функций со свободной спецификацией, остальные двумерные
§ функции неявным образом определяются из структуры построения.
С5
£ Выражение (8) будет «-мерной копула-функцией только тогда, когда, в дополнение к свой! ству абсолютной монотонности для обратных функций генераторов, «композитная» функ-® ция (рт ° ср~1 будет иметь абсолютно монотонные первые производные для вложений всех
о и
0
1 §
уровней г.
Как оказывается, если все генераторы однотипны, степень зависимости, выражае-| мая параметром копула-функции, должна убывать с увеличением уровня вложения, т. е. щ а1 >а2 >...> ап_1. Если же генераторы принадлежат разным семействам, то ограничения ^ на параметры будут еще более жесткими, а список генераторов, которые могут быть свя-® заны, не столь велик (более детальную информацию см., например, в (Whelan, 2004; Savu, § Trede, 2009; McNeil, 2008)).
124 у =
Консультации •
№2(22) 2011
Построение частично вложенных архимедовых копула-функций. Построение час- | тично вложенных архимедовых копула-функций — альтернативный способ обобщения g многомерных копула-функций вида (7а). Эта структура связывания изначально была пред- £ ложена в работе (Joe, 1997), а дальнейший подробный анализ представлен в (Whelan, 2004; в McNeil et al., 2005; McNeil, 2008). Структура этого метода связывания представляет собой 4 смесь структур, используемых при построении перестановочной и полностью вложенной копула-функций, и носит название частично вложенной структуры. Наименьшая размерность, для которой частично вложенная структура связывания дает копула-функции, равна четырем. Для нее:
С(uj, u2 J U3 J U4 ) = C3 (C1 (U1 J U2 )' C2 (U3' U4 )) = ^^
= (p~1 {<Рз ((p~X {<Pi (Ui ) + <Pi (u2)}) + <Рз (<p~1 {<P2 (из ) + (P2 (UA)})}.
Рис. 8. Частично вложенные архимедовы копула-функции
Как и в предыдущих случаях, имеем громоздкие аналитические выражения, но весьма ясное описание: на первом шаге связываем пары (щ,и2) и (и3,и4) копула-функциями Cx (•) и C2 Q, генераторы которых равны , соответственно. Эти две копула-функции связы-
ваются при помощи третьей копула-функции C3 (■). Случайные величины щ и и2 являются перестановочными, так же, как и3 и и4, но все другие пары перестановочными не являются. Однако пары (щ,и3), (щ,и4), (и2,и3) и (и2,и4) имеют копула-функцию C3(•) (см. рис. 8). Для того чтобы получить «-мерную копула-функцию с использованием такого способа связывания, требуются ограничения, похожие на рассмотренные выше при построении полностью вложенных архимедовых копула-функций.
Общий случай: иерархически вложенные архимедовы копула-функции. Этот класс копула-функций впервые был предложен в работе (Joe, 1997), а в дальнейшем подробно изучался в (Whelan, 2004). В работе (Savu, Trede, 2009) представлена первая попытка разработки иерархически вложенных архимедовых копула-функций в общем виде. Такая структура связывания является расширением структуры частичного вложения в том смысле, что участвующие в процедуре связывания копула-функции не обязательно двумерные.
Основная идея структуры связывания, предложенная в (Savu, Trede, 2009), заключается в использовании иерархического вложения. При этом предполагается, что имеется L уровней, по п1 различных объектов на каждом уровне / (объектом может быть либо ко-
№2(22) 2011
s
пула-функция, либо переменная). На уровне / = 1 переменные и1,...,ип разбиваются на п1 групп, в каждой из которых переменные связываются при помощи перестановочной архимедовой копула-функции. В свою очередь, на уровне / = 2 эти копула-функции связываются при помощи п2 копула-функций и т. д. Более формально, на уровне / = 1 величины и1ип связываются при помощи п1 многомерных архимедовых копула-функций С1 .(•), у = 1,___,п1, имеющих вид:
-i
Q, j ("и) = <p¡,J
2n j к J )
\ "и
где <р1 .(•) —генератор копула-функции С1;.0); и1;. (у = 1,...,^) — подмножество элементов и1ип, связываемых при помощи С1 ](•). Копула-функции С11 (•),.„,С1 п (') могут принадлежать разным семействам архимедовых копула-функций, таким как семейство Франка или Гумбеля. Копула-функции уровня / = 1, в свою очередь, агрегируются при помощи копула-функций уровня / = 2 .В результате получаем п2 обобщенных архимедовых копула-функций С2 . (•), у = 1,___, п2, структура зависимости которых является частично перестановочной. Их компоненты — копула-функции предыдущего уровня. Более формально:
2 <Ри (С.,)
C2,J (C2J ) =<Plj
>S S в-
■g где j = i,..., n2; (p2j (•) — генератор копула-функции C2 . (•); C2 . — множество всех копула-
* функций уровня / = i, агрегируемых копула-функцией C2 . (•). Далее продолжаем шаг за ша-
<i т '
^ гом до тех пор, пока не достигнем уровня L, на котором получим единственную иерархи-
| ческуюархимедовукопула-функцию CLi (•). §
ф
| В работе (Savu, Trede, 2009) используются несколько отличные от используемых в данной
о статье, но эквивалентные, обозначения: посколькуу'-ая копула-функция /-го (/ = i,...,L) уров-
еЦ ня C¡j(•) имеет в качестве аргументов u7 ., т.е. подмножество из uiun, элементы кото-
g poro в качестве аргументов входят явным или неявным образом в C¡ ■ (•), или (что эквива-
о лентно) аргументы С, . (•) (т. е. множество всех копула-функций уровня I — 1, участвующих
| в С1;]О), то С1;](С/;]) и С1;](и7;]) суть одно и то же.
| Для того чтобы иметь хорошо определенную иерархию, число используемых на каж-
5 дом уровне копула-функций должно убывать с увеличением уровня, т. е. п1 < п1_1 для всех
5 1 = 2,...,Ь, а на самом верхнем уровне должен использоваться лишь один агрегирующий объ-
С5
£ ект Сь 1 (•), т. е. пь =1. Кроме того, размерность используемых копула-функций на каждом
! следующем уровне должна увеличиваться. На самом верхнем уровне размерность СЬ1 (•)
® равна размерности исходных данных п,т. е. пЬ1 = п .
g ^^^^^^ ^^^^^^^^^^^^^^^ диш хш/v /t^ i, w. Li
o
U
0
1 §
Пример 3. Четырехмерные частично вложенные архимедовы копула-функции, описан-¡5 ные в (9) (см. рис. 8), имеют вид:
s
та о
о &
§ 5
C2,i(") _ C2,i(Ui' U2 ' U3> U4 4 — C2,i (Ci,i(Ui) U2 Ci,2 CU3' U4)) _
= <РЦ (<P2,i ° <P\l [<Pi,i (Ui ) + ^i,i(U2)]+ ^2,i [<Pi,2 (U3 ) + <Pi,2 (U4 )])) ■•
№2(22) 2011
Пример 4. Рассмотрим 9-мерную структуру, определяемую на рис. 9. Копула-функция |
для нее имеет вид: §
¡5
I
С (Wj,..и9) — С4д (С3д (С2ц (С1;1 (Wj, u2), u3, u4), u5, u6), C22(u7, C1;2 (u8, u9))).
Рис. 9. Иерархически вложенная архимедова конструкция
В работе (Savu, Trede, 2009) указаны некоторые условия для того, чтобы иерархическое вложение давало многомерную функцию распределения. Во-первых, требуется, чтобы все обратные функции генераторов (р^ О были абсолютно монотонными. Во-вторых, для / = 1,...,L и j = 1,___,nl, i = 1,___,п1+1 сложные функции <pl+1J ^^ должны иметь абсолютно монотонные производные при всех / = 1,2,...,L и j = 1,___,п1. В работе (Ет-
brechts et al., 2003) показано, что в случае полного вложения копула-функций либо из семейства Гумбеля, либо семейства Клейтона степень зависимости, выраженная параметром а, должна убывать с увеличением уровня иерархических вложений. В примере 4 это означает, что а1Д > а2Д > а3д > а4Д и а12 > а2 2 > а4Д.
Если же агрегируемые копула-функции принадлежат разным семействам архимедовых копула-функций, то, как показано в (McNeil, 2008), две архимедовы копула-функции из двух различных семейств а и b могут быть иерархически вложены только в том случае, если производная «композитной» функции сра о (р~1 является абсолютно монотонной. Вопрос о том, какие семейства копула-функций могут быть агрегированы, в некоторой степени рассмотрен в работе (Joe, 1997).
В силу достаточно сложной структуры иерархических копула-функций, их плотности задаются весьма непростыми выражениями. В (Savu, Trede, 2009) используется рекурсивный подход для дифференцирования «-мерной копула-функции CL1 (•) верхнего уровня по аргументам u L1 .А именно, в силу того, что
№2(22) 2011
C1,1 (U1,1) C1,1 (C1,1 ) C1,1 (C 1-1,1' • " ' C1-1,n1_1 ) C1,1 (C1-1,1 (CL-1,1 )' • " ' C1-1,n1_1 L-l,nL-i ^ ' = CL,1 (C 1-1,1 (uL-1,1 )»-» CL-1,n1_1 (UL-1,nи )),
плотностьдля CL1 (uL1) можетбытьвычисленаследующимобразом:
3nCL, (u ,1 ) v C,, (u l, )
C1,1 (U 1,1) =■
=2-
x
"1-1 _
П 2
r=1 u={o, ,...,or }
дщ ...dun ^dC^^dC^
dWlK-i,r (u l,i ) 3|t7rlCi_i;r (и,,)
X
3a,
3a r
(10)
где внешняя сумма берется по всем к1,...,к ЕNи{0} таким, что шах^. к^ <п1_ и ^ ^ = я _ г' Для всех г'= 0,^, я - пь_,.
>s
S
I
та
t §
§
ф
!
та о
0
та
л §
С U S
U
>s
S
1
ф
5
5?
6
U
та а.
л I
а
ф §
о
0
1
5 ф
!
та о
о
6
S 5
Вторая часть формулы включает в себя производные копула-функций уровня 1 — 1 с аргументами uL_Xj, j = 1,...,nL_,. Суммирование во второй части производится по r Е (0,1,...,nL_1} различным подмножествам [о,,...,оr} множества uL_, r, а произведение берется по всем порядковым номерам копула-функций уровня L — 1.
Алгоритм нахождения плотности n-мерной копула-функции CL1 (•) является рекурсивным: плотность для CL, (•) определяется через частные производные копула-функций CL_,, (•),... ,CL_, i (•), которые в свою очередь могут быть вычислены с использованием (10) для копула-функций, используемых на более низких уровнях. Рекурсивная процедура завершается на самом низком уровне, когда требуется найти частные производные (различных порядков) лишь стандартных архимедовых копула-функций.
В силу рекурсивной природы алгоритма нахождения плотности для CL, (•) количество вычислительных операций растет вместе с усложнением копула-функции. На практике для вычисления плотности иерархической n-мерной копула-функции приходится прибегать к компьютерным алгебраическим системам, таким как Mathematica или R (функция D). Нет необходимости упоминать про то, что время, необходимое для вычисления оценок максимального правдоподобия, растет вместе с ростом размерности n (более детально см. (Aas et al., 2009)).
Пример 5. Используя формулу (10), получим, что плотность четырехмерной иерархической архимедовой копула-функции, представленной на рис. 9, определяется соотношением:
34C,
du, ...du4
+
Ci
dC,2,dC122 du, du2
33C21 d2Cv dC1,2
dC, 1dC122 du,du2 du3
dC,~ dC,
du, du.
d3C
dC„ dCn d2C,
dC, 1dC1 2
du.
du2 du3du4
- +
3Q. du.
d2Cn
32C,, 32C,
dC, 1dC12
du1du2
du3du4
№2(22) 2011
3.5. Моделирование наблюдений из архимедовых копула-функций g
!
Прямое применение традиционных методов генерирования многомерных наблюдений, £ подчиняющихся заданному закону распределения, в принципе, возможно и к распреде- в лениям, описанным с помощью копула-функций, см. (Frees, Valdez, 1998; McNeil, 2008). 4 К сожалению, эти методы оказываются чрезмерно трудоемкими в случае вложенных архимедовых копула-функций и многомерных перестановочных архимедовых копула-функ-ций. Для последних моделирование наблюдений целесообразно проводить с использованием метода условной инверсии, который также можно использовать для любой копула-функции.
Если имеется копула-функция C = C(u1,...,un) (не обязательно архимедова), и необходимо сгенерировать наблюдения (ulv..,un) n-мерного распределения случайного вектора (U1,...,Un), копула-функция которого равна C(•), а частные распределения являются равномерными на отрезке [0,1], то можно воспользоваться методом условного распределения. Пусть Ck(u1 ,u2,...,uk) = C(u1,u2,...,uk,1,...,1), к = 1,...,n ; при этом C1 (u1) = u1 и Cn(u1,u2,...,un) = C(u1 ,...,un). Условное распределение Uk, к = 2,...,n при заданных U1,^,Uk_1 равно:
Ck(ut lu1 u) = p(Uk £Щ = «1 ,...,Uk_1 = uk_1) = = Ck (u1 'u2 v>uk V Ck-1 (u1 'u2 V- Iuk-1 )•
Приведем алгоритм моделирования:
• смоделировать n независимых равномерно распределенных случайных величин v1,..., vn;
• положить u1 = v1;
• аналитически или с использованием численных методов найти обратные функции ус-ловныхраспределений Ck(u1,u2,...,uk) при к = 2,...,п.
Результатом этого алгоритма является вектор наблюдений (u1,...,un), смоделированный в соответствии с распределением C(•). Несмотря на то что моделирование с использованием условных распределений весьма элегантно, сама процедура может быть весьма трудоемкой с вычислительной точки зрения.
В случае архимедовых копула-функций, условное распределение Uk при заданных значениях U^...,ик-1 определяется теоремой 4.
Теорема 4. Пусть C(u1,u2,...,un) = cp~l(<p(u1) + <p(u2) + ... + (p(un)) есть n-мерная архимедова копула-функция с генератором <р(;). Тогда для к = 2,..., п
<p~l( k-1>(<p(u 1) + ф2) + ... + <p(uk))
Ck \uk lu1 'u2 V5u*-1 > - -Kk-1)( , 4, , 4, , , \\ •
<P (<P(u 1 ) + <P(u 2 ) + --- + <P(uk~1 ))
Доказательство. См. (Cherubini et al., 2004, c. 183).
В работе (Marshall, Olkin, 1988) предложен метод построения архимедовых копула-функ-ций с использованием преобразования Лапласа и обратной к нему функции. Этот метод
129
№2(22) 2011
особенно полезен, учитывая то, что каждое абсолютно монотонное отображение [О,»] в [0,1] можно представить в терминах преобразования Лапласа и обратной к нему функции. Пусть О(•) — функция распределения на М+ такая, что С(0) = 0, а ее преобразование Лапласа-Стилтьеса определяется соотношением
G (t) = /<
e~txdG(x), t > 0.
Можно показать, что G: [0,^] ^ [0,1] является непрерывной, абсолютно монотонной и строго убывающей функцией. Поэтому G(•) — хороший выбор для обратной функции генератора архимедовой копулы.
Как показано в работах (Frees, Valdez, 1998; Cherubini et al., 2004; McNeil et al., 2005), с использованием преобразования Лапласа-Стилтьеса можно предложить следующий алгоритм моделирования наблюдений из распределения, описанного многомерной архимедовой копула-функцией:
• смоделировать случайную величину V из распределения G(•) такого, что G(•) (преобразование Лапласа функции G(•)) является обратной функцией генератора ср(-) копула-функции, из которой моделируются наблюдения. Например:
— для копула-функции Клейтона Vимеет гамма-распределение Ga( 1 /а, 1), а>0,
^ a G(t) = (1 +1)~l1a ; отметим, что G = (t~a — 1) отличается от генератора под номером 1 | из табл. 2 лишь на константу;
Л — для копула-функции Гумбеля V имеет устойчивое распределение St (1 / а,1,у, 0)с
I у = cosа(па/2),д>1и G(t) = exp(-t1/a); §
s — длякопула-функцииФранка Гимеетдискретноераспределение
i .
S P(V = к) = (1-)* /(ка) при к = 1,2,... и а>0.
ч • сгенерировать независимые одинаково распределенные случайные величины X1,...,Xn.
u = Gпт^, i = 1,...,
о с о 5
0 • ВЗЯТЬ Ы1 »
| Вектор (и1;...,ип)т смоделирован из «-мерной архимедовой копула-функции. Отметим,
Э1 что для рассмотренных выше случаев алгоритм предполагает необходимость моделирова-
§ ния случайных величин из гамма-распределения, устойчивого распределения или дискретен а.
л
1
а.
ф §
о и
0
1 §
ного распределения. Более подробную информацию см. в (Cherubini et al., 2004; Schoutens, 2003; McNeil et al., 2005; McNeil, 2008).
3.6. Эмпирические приложения в статистическом пакете R: s архимедовы копула-функции
^ Для того чтобы вычислить значения функции распределения двумерной архимедовой ® копула-функции и ее плотности, необходимо воспользоваться некоторыми процедурами § из модуля геори1ае в пакете Я:
130
I №
2(22) 2011
fcopulae функция
Описание
!
i
m
е et
rarchmCopula
parchmCopula darchmCopula
rarchmSlider
parchmSlider darchmSlider archmCopulaFit
моделирование случайных величин из архимедовой копула-функции
вычисление значения архимедовой копула-функции
вычисление значения плотности архимедовой копула-функции
график наблюдений, смоделированных из архимедовой копула-функции
график функции распределения
график функции плотности
оценка параметров архимедовой копула-функции
Например, для того чтобы смоделировать наблюдения из двумерной копула-функции Гумбеля, а затем по ним оценить параметры той же копула-функции, необходимо использовать код:
# Random Variâtes:
R = archmCopulaSim (n = 1000, alpha = 1, type = "4")
# Fit:
fit = archmCopulaFit (u=R [, 1], v=R [, 2], type = "4") fit
Если же рассматривать многомерную перестановочную архимедову копула-функцию, то необходимо прибегнуть к использованию модуля copula. Ниже представлен пример кода для трехмерной перестановочной архимедовой копула-функции (для более детальной информации рекомендуем обратиться к руководству статистическим пакетом R), а результаты моделирования изображены на рис. 10.
rcopula(clayton.cop, 1000)[,1]
Рис. 10. График смоделированных наблюдений из трехмерной копула-функции Клейтона
№2(22) 2011
library(copula)
#Some graphical examples
v <- rcopula(claytonCopula(2, dim = 3), 1000)
scatterplot3d(v)
frank.cop <- frankCopula(3)
persp(frank.cop, dcopula)
gumbel.cop <- archmCopula("gumbel", 5)
contour(gumbel.cop, dcopula)
# A 5-dim Frank copula
frank.cop <- frankCopula(3, dim = 5)
Для моделирования и оценки четырехмерной полностью вложенной архимедовой копу-лы можно было бы обратиться к модулю «copula GOF» (недавно разработанному Дэниелем Бергом). Однако к моменту написания данной работы этот модуль находился в состоянии «уЗ-версии» и не был общедоступным, т. е. он еще не был загружен в среду пакета R.
Список литературы
Aas К., Czado С., Frigessi A., Bakken Н. (2009). Pair-copula constructions of multiple dependence. Insurance: Mathematics andEconomics, 44 (2), 182-198.
'§ Alsina M., Frank J., Schweizer B. (2005). Problems on associative functions. Aequationes Mathemati-
f cae, 66,128-140.
■fr Ang A., Chen J. (2002). Asymmetric correlations of equity portfolios. Journal of Financial Economics,
! 63 (3), 443-494. с
g Bouyé E., Durrleman V., Nikeghbali A., Riboulet G., Roncalli T. (2000). Copulas for finance a reading
! guide and some applications. Groupe de Recherche Operationnelle, Credit Lyonnais, WorkingPaper. 5
и Breymann W., Dias A., Embrechts P. (2003). Dependence structures for multivariate high-frequency data
° infinance. QuantitativeFinance, 3, 1-14.
•e
§ Chen X., Fan Y., Patton A. (2004). Simple tests for models of dependence between multiple financial
5 time series, with applications to U. S. equity returns and exchange rates. London Economics Financial Maris ketsGroup, WorkingPapern. 483.
s
IE Cherubini U., Vecchiato W., Luciano E. (2004). Copula methods infinance. Wiley. Ц Clayton D. G. (1978). A model for association in bivariate life tables and its application in epidemiologist cal studies offamilial tendency in chronic disease incidence. Biometrika, 65, 141-151. о
a Daul S., De Giorgi E., Lindskog F., McNeil A. (2003). The grouped t-copula with an application to credit
3 risk. Risk, 1, 73-76.
® Demarta S., McNeil A. (2005). The t copula and related copulas. International Statistical Review, 73,
2 111-129.
о
§ Embrechts P., McNeil A. J., Straumann D. (1999). Correlation and dependency in risk management: Proper-
Is ties and pitfalls. Working paper, Department of Mathematik, ETHZ, Zurich. (Now in M. A. H. Dempster (ed.)
о (2002), RiskManagement: ValueatRiskandBeyond, 176-223. Cambridge: CambridgeUniversityPress).
6 Embrechts P., Lindskog F., McNeil A. (2003). Modelling dependence with copulas and applications to
4 risk management. In: Rachev S. T. (ed.), Handbook ofHeavy Tailed Distributions in Finance, Amsterdam: § Elsevier/North-Holland.
№2(22) 2011
Erb C., Harvey C., Viskanta T. (1994). Forecasting international equity correlations. Financial Analysts s
Journal, 50, 32-45. I
5
Fang K., Kotz S., Ng K. (1987). Symmetric multivariate and related distributions. London: Chapman Hall. E
Fantazzini D. (2009a). A dynamic grouped T copula approach for market risk management. In: G. Gre- § goriou (ed.), A VaR Implementation Handbook, 253-282, McGraw-Hill: New York. ct
Fantazzini D. (20096). The effects of misspecified marginals and copulas on computing the value at risk: AMonte Carlo study. ComputationalStatistics andDataAnalysis, 53 (6), 2168-2188.
Fantazzini D. (2010). Three-stage semi-parametric estimation of T-copulas: Asymptotics, finite-sample properties and computational aspects. Computational Statistics and Data Analysis, forthcoming.
Fermanian J., Scaillet O. (2003). Nonparametric estimation of copulas for time series. Journal of Risk, 5,25-54.
Frank M. J. (1979). On the simultaneous associativity of F(x, y) and x + y — F(x, y). Aequationes Mathematicae, 19, 194-226.
Frees E. W., Valdez E. (1998). Understanding relationship using copulas. North American Actuarial Journal, 2, 1-25.
Genest C. (1987). Frank's family ofbivariate distributions. Biometrika, 74 (3), 549-555.
Genest C., Favre A.-C. (2007). Everything you always wanted to know about copula modeling but were afraidto ask. Journal ofHydrologicEngineering, 12 (4), 347-368.
Genest C., Rivest L. (1993). Statistical inference procedures for bivariate archimedean copulas. Journal ofthe AmericanStatisticalAssociation, 88, 1034-1043.
Genest C., Ghoudi K., Rivest L.-P. (1995). A semiparametric estimation procedure of dependence parameters inmultivariate families ofdistribution. Biometrika, 82 (3), 543-552.
Gumbel E. J. (I960). Bivariate exponential distributions. Journal of the American Statistal Association, 55, 698-707.
Hoeffding D. (1940). Masstabinvariante Korrelationstheorie. Schriften des Mathematischen Seminars und des Institutsfur Angewandte Mathematik der Universität, 5, 181-233.
Hougaard P. (1986). A class ofmultivariate failure time distributions. Biometrika, 73, 671-678.
Joe H. (1997). Multivariate models and dependence concepts. London: Chapman Hall.
Jondeau E., Rockinger M. (2003). Conditional volatility, skewness, and kurtosis: existence, persistence, and comovements. Journal ofEconomic Dynamics and Control, 27, 1699-1737.
Kimberling C. H. (1974). A probabilistic interpretation of complete monotonicity. Aequationes Mathematicae, 10, 152-164.
Lambert P., Vandenhende F. (2002). A copula-based model for multivariate non-normal longitudinal data: analysis ofadose titration safety study on a new antidepressant. Statistics in Medicine, 21,3197-3217.
Lindskog F., McNeil A., Schmoch (2002). A note on Kendall's tau for elliptical distributions. ETHpreprint.
Lindskog F. (2000). Modelling dependence with copulas and applications to risk management. Master Thesis, ETH Zurich.
Longin F., Solnik B. (2001). Extreme correlation of international equity markets. Journal of Finance, 56 (2), 649-676.
Malevergne Y., Sornette D. (2006). Extreme financial risks (From dependence to risk management). Springer: Heidelberg.
MarshallA., Olkin I. (1988). Families ofmultivariate distributions. Journal of the American Statistal Association, 83,834-841.
№2(22) 2011
Marshall R., Zeevi A. (2002). Beyond correlation: Extreme co-movements between financial assets. Columbia University, Working Paper.
McNeil A. (2008). Sampling nested Archimedean copulas. Journal of Statistical Computation and Simulation, 78 (6), 567-581.
McNeil A., Neslehová J. (2009). Multivariate Archimedean copulas, d-monotone functions and Ll-norm symmetric distributions. Annals ofStatistics, 37 (5), 3059-3097.
McNeil A., Frey R., Embrechts P. (2005). Quantitative risk management: Concepts, techniques and tools. New Jersey: Princeton Series in Finance.
Mills F. C. (1927). The behaviour of prices. New York: National Bureau ofEconomic Research.
Müller A., Scarsini M. (2005). Archimedean copulae andpositive dependence. Journal ofMultivariate Analysis, 93, 434-445.
Nelsen R. (1999). An introduction to copulas. Lecture Notes in Statistics. New York: Springer-Verlag.
Nelsen R. B. (2005). Some properties of Schur-constant survival models and their copulas. Brazilian JournalofProbabilityandStatistics, 19, 179-190.
Nelsen R. B. (2006). An introduction to copulas. Lecture Notes in Statistics, 2nd Edition. New York: Springer-Verlag.
Patton A. (2004). On the out-of-sample importance of skewness and asymmetric dependence for asset allocation. JournalofFinancialEconometrics, 2 (1), 130-168.
Patton A. (2006a). Estimation of copula models for time series ofpossibly different lengths. Journal of AppliedEconometrics,2\, 147-173.
Patton A. (20066). Modelling asymmetric exchange rate dependence. International Economic Review,
>s s
r
* 47 (2), 527-556
i Sí
Salvadori G., De Michele C. (2007). On the use of copulas in hydrology: Theory and practice. Journal ofHydrologicEngineering, 12 (4), 369-380.
Savu C., Trede M. (2009). Hierarchies of Archimedean copulas. Quantitative Finance, forthcoming. Schmidt R. (2002). Tail dependence for elliptically contoured distributions. Mathematical Methods of
n Operations Research, 55 (2), 301-327.
SchweizerB., SklarA. (1983). Probabilistic metricspaces. New York: Elsevier.
Schweizer B., Wolff E. F. (1976). Sur une mesure de dependence pour les variables aléatoires. Comptes Rendus de l'Academie des Sciences de Paris. Ser. A, 283, 659-661.
Schweizer B., Wolff E. (1981). On non-parametric measures of dependence for random variables.^«-nals ofStatistics, 9, 879-885.
SklarA. (1959). Fonctions de répartition a n dimensions et leurs marges. Publ. Inst. Statis. Univ. Paris, 8,229-231.
SklarA. (1996). Random variables, distribution functions, and copulas: Personal look backward and forward. Lecture notes. Monograph series, 28, 1 -14.
Zhang L., Singh V. (2006). Bivariate flood frequency analysis using the copula method. Journal ofHy-drologicEngineering, 11 (2), 150-164.
Wang S. (1998). Aggregation of correlated risk portfolios: Model and algorithms. Proceedings of the Casualty ActuarialSociety, LXV, 848-893.
WhelanN. (2004). Sampling from Archimedeancopulas. QuantitativeFinance, 4 (3), 339-352.
134
