Санкт-Петербургский государственный университет
Кафедра математической теории игр и статистических
решений
Арасланова Валерия Альбертовна
Выпускная квалификационная работа бакалавра
Применение копула-функций в анализе данных
Направление 010400
Прикладная математика и информатика
Научный руководитель,
доктор технич. наук,
профессор
Буре В.М.
Санкт-Петербург
2016
Содержание
Введение
3
Постановка задачи
4
Обзор литературы
5
1
Математическая теория по копула-функциям
1.1. Копула-функция и некоторые ее свойства . . . . . . . . . . . .
1.2. Плотность копула-функции . . . . . . . . . . . . . . . . . . . .
1.3. Виды взаимодействия случайных величин . . . . . . . . . . .
7
7
10
10
2
Методология анализа данных с помощью копула-функции
12
3
Практическое применение копула-функции для анализа временных рядов
3.1. Анализ финансового рынка России . . . . . . . . . . . . . . .
3.2. Анализ секторов российского финансового рынка . . . . . . .
3.3. Взаимодействие рядов курсов акций при различных временных лагах . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
14
14
16
25
Выводы
28
Заключение
29
Список литературы
30
Приложение
35
2
Введение
В настоящее время все более актуальным становится вопрос о получении, подготовке и анализе данных. Копула-функция представляет собой
инструмент для выполнения последнего.
Метод копула-функции определяет характер взаимодействия между
случайными величинами, что позволяет сформировать более адекватную
модель описания совместного распределения нескольких переменных, чем
многомерный нормальный закон [1]. Важным преимуществом функции является то, что при моделировании многомерного распределения, помимо
частных распределений она учитывает характер их взаимодействия. Копулафункция может быть применена для описания совместного распределения
переменных имеющих разные частные распределения.
Широкое применение метод находит в исследовании финансовых рынков, где важно следить за динамикой курсов акций различных компаний, а
также за существующей между ними корреляцией. Копула-функция применима к оценке состояния какого-либо экономического сектора или экономики страны в целом. Помимо этого, с помощью нее возможно исследование
на поиск временных лагов.
В данной работе представлены возможные решения перечисленных
задач с использованием копула-функции. Областью применения был выбран
финансовый рынок России.
3
Постановка задачи
Основной задачей работы является поиск разнообразных применений
копула-функции, а также реализация метода копулы с получением практических результатов.
В связи с поставленной задачей были выделены подзадачи:
1. Проанализировать общее состояния финансового рынка в России.
2. Проанализировать характер взаимодействия компаний внутри отраслей. Сравнить поведение различных секторов на временных промежутках.
3. Применить метод копула-функции к задаче поиска временного лага.
Результатом проделанного исследования является:
1. Характеристика состояния финансового рынка России;
2. Анализ взаимодействия фирм внутри отрасли, поиск компаний и
секторов наиболее устойчивых к негативному влиянию со стороны рынка;
3. Выявление временного промежутка, по истечению которого изменение цены на нефть оказывает влияние на стоимость бензина.
4
Обзор литературы
Понятие копула-функции относительно новое в математике. Функция
подробно изучалась в работах [2,3], кроме того, нужно отметить статью [4].
В русскоязычной литературе копула получила распространение только в начале 21 века. Для ознакомления с теорией по рассматриваемой функции
лучше всего подходят работы [1, 5], также в [1] описаны типы копулафункций и методы оценивания их параметров. Обзор наиболее распространенных копул приведен в [6]. Многие алгоритмы подразумевают правильный выбор модели копула-функции, методы подбора функций разобраны
в [7].
Копула-функция имеет широкое применение в финансовом моделировании. На сегодняшний день является актуальной задача формирования
портфеля ценных бумаг, разработки по данной теме предложены в [8–12].
В частности в [13] построение портфеля реализуется с помощью парных
копул, соответствующая теория по последним приведена в статье [14]. В
этой области важным является правильная оценка границ потерь портфеля
(так же здесь содержится теория по рискам) [15]. Предложенная методика в [16] применяется для оценки риска кредитного портфеля московского
коммерческого банка.
Задача формирования портфеля ценных бумаг является более узкой
задачей экономики. В широком смысле копула-функция используется для
оценки взаимосвязи курсов акций [17], для обнаружения структурных сдвигов [18] или для анализа общего состояния финансового рынка [19].
Помимо экономической сферы с помощью копула-функции реализуется классификация [20, 21], осуществляется статистический подход в области физики (исследование массового состава) [22], копула используется в
5
страхования [23], сельском хозяйстве [24, 25], а также функция нашла свое
применение в задачи внедрения сетей сотовой подвижной связи для оценки
суммарного помехового сигнала [26].
Все перечисленное показывает насколько разнообразен спектр применения копулы, что доказывает целесообразность ее изучения. В данной
работе будут разобраны некоторые задачи с применением рассматриваемой
функции.
6
1
1.1.
Математическая теория по копула-функциям
Копула-функция и некоторые ее свойства
Определение 1: [1] Копула-функция C(u1 , u2 .., un ) - это функция n
переменных, определенная на единичном кубе U n = [0, 1]n (т.е. ui ∈ [0, 1],
i = 1, 2, ..., n), обладающая следующими свойствами:
1) область значений функции - единичный интервал [0,1];
2) если ui = 0 по крайней мере для одного i ∈ 1, 2, ..., n, то
C(u1 , u2 , ...un ) = 0;
3) C(1, ..., 1, ui , 1, ..., 1) = ui , для ∀ ui ∈ [0, 1];
4) C(u1 , u2 .., un ) является n-возрастающей функцией в том смысле,
что для ∀ (a1 , a2 , ..., an ), (b1 , b2 , ..., bn ) ∈ [0, 1]n с ai ≤ bi справедливо неравенство:
2
X
i1 =1
...
2
X
(−1)i1 +...+in C(u1i1 , ..., unin ) ≥ 0,
in =1
где uj1 = aj и uj2 = bj для ∀ j ∈ 1, ..., n.
Теорема 1 (Теорема Скляра): [1] Пусть H(·) - n-мерная функция распределения, с частными распределениями F1 (·), ..., Fn (·). Тогда существует
n-мерная копула-функция C(·) такая, что для всех действительных x1 , ..., xn :
H(x1 , ..., xn ) = C(F1 (x1 ), ..., Fn (xn ))
(1)
Если F1 , ..., Fn непрерывны, то копула-функция определена единственным образом;
Обратно, если C(·)-копула-функция, а F1 (·), ..., Fn (·)- функции распределения, то функция H(·), определяемая выражением (1), является совместной функцией распределения с частными распределениями F1 (·), ..., Fn (·).
7
Доказательство: См. [3–5].
Следствие 1: [1] Пусть F1−1 (·), ..., Fn−1 (·) - обратные функции частных
распределений. Тогда для любого (u1 , ..., un ) ∈ U n существует единственная
копула-функция C : [0, 1] × ... × [0, 1] → [0, 1] такая, что:
C(u1 , ..., un ) = H(F1−1 (u1 ), ..., Fn−1 (un ))
Доказательство: См. [4].
Теорема 2: [1] Рассмотрим n случайных величин X1 , ..., Xn , зависимость между которыми определяется копула-функцией C(·). Если преобразования Ti : R → R, i = 1, .., n определяются строго возрастающими
функциями (т.е. ∂x Ti > 0), то структура зависимости случайных величин
T1 (X1 ), ..., Tn (Xn ) определяется той же самой копула-функцией C(·).
Доказательство: См. [27, 28], для двумерного случая — [29].
Границы Фреше –Хёффдинга: [1]
W (u1 , ..., un ) ≤ C(u1 , ..., un ) ≤ M (u1 , ..., un ),
где
W (u1 , ..., un ) = max(u1 + ... + un − n + 1, 0);
M (u1 , ..., un ) = min(u1 , ..., un )
При n = 2 границы являются копула-функциями, однако при n > 2 нижняя
граница не будет n-возрастающей.
Доказательство: См. [4, 5]
8
Условия Липшица для копула-функции: [5]
Любая копула удовлетворяет условию Липшица:
| C(s1 , ..., sn ) − C(t1 , ..., tn ) |≤| s1 − t1 | +...+ | sn − tn |
Доказательство: См. [5]
Из свойства Липшица и теорем о монотонных функциях следует, что
любая копула-функция почти всюду дифференцируема по каждой из переменных [1, 4, 5].
Марковское свойство копула-функции: [5]
Рассмотрим копулы A(s1 , ..., sm ) и B(t1 , ..., tn ) размерностей m и n.
Пусть
Z
A0m (u1 , ..., um−1 , x)B10 (x, um+1 , ..., um+n−1 )dx,
C(u1 , ..., um+n−1 ) =
0≤x≤um
где
∂
Am (u1 , ..., um−1 , x)
∂x
∂
B1 (x, um+1 , ..., um+n−1 )
B10 (x, um+1 , ..., um+n−1 ) =
∂x
A0m (u1 , ..., um−1 , x) =
что кратко можно записать как результат произведения копул:
Cm+n−1 (·) = Am (·) ◦ Bn (·)
Пусть A(u | t1 , ..., tn ), uT = (u1 , ..., up ) - копула для частных распределений X(t1 ), ..., X(tn ). Тогда X(t) является марковским процессом тогда и
только тогда, когда для любого n ≥ 2 и любых t1 < ... < tn :
A(u | t1 , ..., tn ) = A(· | t1 , t2 ) ◦ A(· | t2 , t3 ) ◦ ... ◦ A(· | tn−1 , tn )
Этот результат можно найти в работе [30].
9
1.2.
Плотность копула-функции
Определение 2: Если для копула-функции C(·) существует неотрица-
тельная функция c(·), для которой выполняется:
Zu1 Zu2
C(u1 , ..., un ) =
Zun
...
0
0
c(s1 , ..., sn )ds1 ds2 ...dsn
0
тогда эта функция c(·) называется плотностью копула-функции.
Из следствия к теореме Склара выводится каноническое представление плотности: [1]
n
n
Y
∂ n C(F1 (x1 ), ..., Fn (xn )) Y
· fi (xi ) = c(F1 (x1 ), ..., Fn (xn ))· fi (xi ),
fH (x1 , ..., xn ) =
∂F1 (x1 )...∂Fn (xn )
i=1
i=1
где
c(F1 (x1 ), ..., Fn (xn )) - плотность копула-функции;
fH (x1 , ..., xn ) - плотность совместной функции распределения H(·).
тогда
fH (x1 , ..., xn )
.
c(F1 (x1 ), ..., Fn (xn )) = Qn
i=1 fi (xi )
1.3.
Виды взаимодействия случайных величин
Следующие определения представлены в работе [17].
1) Случайные величины независимы, если их копула-функция независима:
C ⊥ (u1 , u2 , ..., un ) = u1 u2 ..un ;
2) Если рост одной случайной величины сопровождается ростом другой случайной величины, то они комонотонны и им соответствует копулафункция:
Cmax (u1 , u2 , ..., un ) = min(u1 , u2 , ..., un );
10
3) Если рост одной случайной величины сопровождается убыванием
другой случайной величины, то они контрмонотонны и им соответствует
копула-функция:
Cmin (u1 , u2 , ..., un ) = max(u1 + u2 + ... + un − n + 1, 0).
11
2
Методология анализа данных с помощью копулафункции
Пусть имеются наблюдения двух временных рядов (xk , yk ), k = 1..n,
1
N . За оценку значения
n(i,j)
n , где n(i, j)-количество
где n-число наблюдений. Выберем шаг решетки
копула-функции C( Ni , Nj ) принимаем величину
наблюдений (xk , yk ), для которых выполняются неравенства:
xk ≤ x([ in
N ])
yk ≤ y([ jn ])
N
где x(k) , y(k) -порядковые статистики. Для того, чтобы выяснить характер взаимодействия временных рядов вычислялись дискретные приближения
N
N X
X
| C(
i1 =1 j1 =1
i j
i j
, ) − C ∗( , ) |
N N
N N
(2)
поочередно используя в качестве эталонной функции C ∗ копула-функции
C ⊥ , Cmax , Cmin . Из трех отклонений выбиралось наименьшее. Если наименьшему отклонению соответствовало C ⊥ , то временные ряды считались
независимыми, если Cmax , то зависимость комонотонная, если Cmin , то
контрмонотонная. [17]
Для более подробного анализа характера взаимодействия случайных
величин введем обозначения:
dmin - отклонение, соответствующее копуле Cmin
d⊥ - отклонение, соответствующее копуле C ⊥
dmax - отклонение, соответствующее копуле Cmax .
Затем определим классы, основанные на ранжировании полученных
отклонений:
1 класс: dmax ≤ d⊥ ≤ dmin ;
12
2 класс: dmax ≤ dmin ≤ d⊥ ;
3 класс: d⊥ ≤ dmax ≤ dmin ;
4 класс: d⊥ ≤ dmin ≤ dmax ;
5 класс: dmin ≤ dmax ≤ d⊥ ;
6 класс: dmin ≤ d⊥ ≤ dmax .
Построенные классы, по мнению [19], выражают относительную интенсивность проявления свойств того или иного случая взаимодействия. Так
же авторы утверждают, что использование классов актуально при формировании портфеля ценных бумаг. В этом случае необходимо включать в портфель акции компаний, у которых копулы пар доходностей относятся к 6
классу, затем к 5 классу, и т.д..
13
3
Практическое применение копула-функции для
анализа временных рядов
В представленных ниже исследованиях используется алгоритм копула-
функций, основанный на подсчете отклонений полученной копула-функции
от эталонных по формуле (2) из предыдущего параграфа.
3.1.
Анализ финансового рынка России
В данном исследовании для анализа финансового рынка производится
подсчет количества контрмонотонных, независимых и комонотонных связей между курсами акций компаний произвольных отраслей, в частности:
Аптека 36и6, Магнит, Газпром, Лукойл, Роснефть, Сургутнефть, Татнефть,
Аэрофлот, Ютейр, Возрождение, ВСМПО-АВИСМА, ГМКНорНикель, СеверСталь, ПолюсЗолото, НЛМК, Иркут-3, ММК, Распадская, ИркутскЭнерго, МосЭнерго, ЭнелРос, Новатек, Ростелеком, МТС, Разгуляй.
Таблица 1: Количество независимых, контр-, комонотонных связей
2007
2008
2009
2010
2011
2012
2013
2014
2015
1
82
288
271
136
208
149
99
77
42
2
0
0
0
0
0
0
0
0
0
3
130
12
5
127
35
68
121
132
153
4
72
0
5
32
32
68
55
79
95
5
0
0
0
0
0
0
0
0
0
6
16
0
19
5
25
15
24
12
10
В качестве данных выбраны ежедневные цены закрытия акций компаний за период 01.01.2007-31.12.2015. Материалы взяты с сайта [31]. Расчеты
14
выполнены в среде Матлаб.
Вычисления производились с параметром N=5 исходя из ежегодного
объема выборки (в среднем 250 наблюдений). Результаты представлены в
Таблице 1.
Из Таблицы 1 видно, что в 2007 году комонотонный и независимый характер взаимодействия курсов акций имеет почти равное количество.
Действительно, в этом году темп роста экономики России имел лучший показатель с конца 90-х [32]. Однако этот период докризисный, о чем свидетельствует комонотонный характер. В 2008-2009 годах в связи с мировым
кризисом имеет место явное преобладание 1 класса зависимости. К 2010
году ситуация на российском рынке улучшилась, что подтверждает рост
независимого характера взаимодействия российских компаний.
В обзоре [33] авторы утверждают, что в 2011 году произошел большой
отток капитала из страны, что повлияло на фондовый рынок, а также приводят статистику о том, что фондовые индексы акций российского рынка
потеряли 20%. В полученных результатах эти события выражаются резким
возрастанием числа комонотонных связей. В 2012-2013 годах в основном
наблюдается независимый характер с преобладанием комонотонности, в это
время падение демонстрируют отрасли электроэнергетики, металлургии и
горнодобычи. Наиболее доходным становится сектор потребительских товаров, также нефтегазовая отрасль, которая в 2012 году занимает 50% капитализации [34, 35]. В последующий год в стране происходят негативные
события: конфликт в Украине, санкции, падение цены на нефть, однако показатели фондового рынка остаются практически на том же уровне, что и
в предыдущие года. По мнению авторов статьи [36] это обусловлено повышенными оборотами в марте и декабре. Однако как и в 2015 году наблюдается преобладание 3 класса, что указывает на склонность к комонотонному
15
характеру. Важно отметить, что 2 и 5 классы взаимодействия не свойственны российскому рынку.
Исходя из полученных результатов можно сделать вывод о том, что в
основном российский рынок характеризуется комонотонной зависимостью.
Это обосновывается неготовностью российской экономики игнорировать
негативные события, произошедшие за рассматриваемый период.
3.2.
Анализ секторов российского финансового рынка
В данном параграфе приведены результаты сравнения акций компа-
ний, принадлежащих к одной рыночной отрасли и попавших в сотню лидирующих по величине капитала за 2015 год.
а) 1 квартал 2014 год
б) 2 квартал 2014 год
в) 3 квартал 2014 год
г) 4 квартал 2014 год
Рис. 1: Относительная динамика цен акций компаний продовольственного
сектора
16
Рассматриваются 5 секторов: сектор потребительских товаров, сектор нефтегазодобычи и нефтепереработки, сектор телекоммуникации, сектор строительства и недвижимости и сектор электроэнергетики. Также в
конце параграфа сделано сравнение экономики этих отраслей.
а) 1 квартал 2015 год
б) 2 квартал 2015 год
в) 3 квартал 2015 год
г) 4 квартал 2015 год
Рис. 2: Относительная динамика цен акций компаний продовольственного
сектора
В качестве данных выбраны акции за период 2014-2015 года следующих компаний:
• сектор потребительских товаров: Дикси(ПТ1), Лента(ПТ2), М.Видео(ПТ3),
X5 Retail Group(ПТ4), Магнит(ПТ5);
• сектор нефтегазодобычи и нефтепереработки: Газпром(Н1), Лукойл(Н2),
Роснефть(Н3), Сургутнефть(Н4), Татнефть(Н5);
17
• сектор телекоммуникации: СТС Медиа(Т1), Ростелеком(Т2), Мегафон(Т3),
МТС(Т4), VimpelCom(Т5);
• сектор строительства и недвижимости: Группа ЛСР(С1), ГАЛС Девелопмент(С2), МОСТОТРЕСТ(С3), Группа компаний ПИК(С4), AFI
Девелопмент(С5);
• сектор электроэнергетики: ФСК ЕЭС(Э1), РусГидро(Э2), Э.ОН Россия(Э3), ИркутскЭнерго(Э4), МосЭнерго(Э5), ОГК-1(Э6), Россети(Э7),
ТГК-1(Э8).
а) 1 квартал 2014 год
б) 2 квартал 2014 год
в) 3 квартал 2014 год
г) 4 квартал 2014 год
Рис. 3: Относительная динамика цен акций компаний строительного сектора
Сравниваются ежедневные цены акций за период 2014-2015 годов деленных на кварталы. Средний объем выборки 60, исходя из этого выбираем
18
а) 1 квартал 2015 год
б) 2 квартал 2015 год
в) 3 квартал 2015 год
г) 4 квартал 2015 год
Рис. 4: Относительная динамика цен акций компаний строительного сектора
параметр N=3. Источники данных [31,37]. Расчеты выполнены в среде Матлаб. Результаты приведены в Приложении (Таблицы 3-41).
Как мы видим из Таблиц 3-4 в первом полугодии 2014 года для потребительского рынка характерна комонотонная зависимость, действительно,
из Рис.1 мы видим схожую динамику цен: к концу марта подъем, затем постепенное снижение. Во втором полугодии (Таблицы 5-6) возникает значительное количество независимых связей, однако по прежнему сохраняется
комонотонная зависимость.
За 2015 год мы наблюдаем существенное количество независимости.
Возникает контрмонотонная связь с компанией М.Видео(ПТ3), которая, по
нашему мнению, связана с различием товарного ассортимента. Примеча19
а) 1 квартал 2014 год
б) 2 квартал 2014 год
в) 3 квартал 2014 год
г) 4 квартал 2014 год
Рис. 5: Относительная динамика цен акций компаний нефтяного сектора
тельна преимущественная независимость акций магазина Магнит(ПТ5) в
течении двух лет, что скорее всего, стало следствием отличающегося менеджмента, проводимого компанией. Наглядно данные изображены на Рис.2,
результаты полученных вычислений приведены в Таблицах 7-10.
Важно отметить, что в 2014 году в строительном секторе (Таблицы
11-14) нельзя выделить какой-либо характерный вид связи. Здесь в равном
количестве присутствует независимый, ко- и контрмонотонный тип взаимодействия. Аналогичная ситуация наблюдается в 2015 году (Таблицы 15-18).
Однако примечательна контрмонотонная связь в 4 квартале 2014 года (Таблица 14), которую доказывает Рис.3. Группа компаний ПИК(С4) терпит падение стоимости акций, в то время как у остальных фирм происходит рост.
20
а) 1 квартал 2015 год
б) 2 квартал 2015 год
в) 3 квартал 2015 год
г) 4 квартал 2015 год
Рис. 6: Относительная динамика цен акций компаний нефтяного сектора
Так же эта компания показывает независимый характер во 2 квартале 2014
и в трех кварталах 2015 года.
Нефтяной сектор в 2014-2015 годах принимает комонотонной тип взаимодействия (Таблицы 19-26), что так же видно на графиках (Рис.5, Рис. 6).
Эти года были не легкими для этой отрасли. Сказалось резкое падение цены
на нефть, как следствие отток капитала. Изредка независимость показывали
лидеры российского рынка Газпром(Н1) и Лукойл(Н2).
Для сектора телекоммуникации 2014 год характеризуется комонотонной связью (Таблицы 27-29). Здесь во втором квартале отличается динамика VimpelCom(Т5), а в третьем квартале Ростелеком(Т2). Так же Ростелеком(Т2) проявляет контрмонотонность во втором квартале 2015 года, в то
21
а) 1 квартал 2014 год
б) 2 квартал 2014 год
в) 3 квартал 2014 год
Рис. 7: Относительная динамика цен акций компаний телекоммуникационного сектора
время как остальные фирмы имеют комонотонную зависимость. В остальные кварталы 2015 года в рассматриваемом секторе наблюдается независимый характер. Динамика курса акций проиллюстрирована на Рис.6, Рис.7.
Сектор электроэнергетики в первом полугодии 2014 года насыщен комонотонной связью (Таблицы 34-41), далее появляется независимость протяженностью в квартал и вновь возникает комонотонность вплоть до 2 квартала 2015 года. Важно отметить, что несмотря на периоды выраженной комонотонности компания ТГК-1(Э8) демонстрирует независимость на протяжении двух лет (за исключением 4 квартала 2014 года). Во втором квартале
2014 года у фирмы возникает контрмонотонные связи, которые выражают22
а) 1 квартал 2015 год
б) 2 квартал 2015 год
в) 3 квартал 2015 год
г) 4 квартал 2015 год
Рис. 8: Динамика цен акций компаний телекоммуникационного сектора
ся быстрым ростом стоимости акций в конце квартала при незначительном
спаде курса акций других компаний этого сектора.
Из приведенного анализа секторов можно сделать вывод о том, что на
российском рынке 2014 года преобладала комонотонная зависимость, однако в 2015 году большинство связей имеют независимый характер. Сравнивая сектора можно утверждать, что наиболее независимой отраслью является сектор строительства и недвижимости, в то время как на нефтяном рынке
явно выражено комонотонное взаимодействие курсов акций. Несмотря на
это из некоторых секторов удалось выделить компании, политика которых
позволяет им держать индивидуальную динамику стоимостей акций. Таковыми являются Магнит(ПТ5) из потребительского рынка, Группа компаний
23
а) 1 квартал 2014 год
б) 2 квартал 2014 год
в) 3 квартал 2014 год
г) 4 квартал 2014 год
Рис. 9: Относительная динамика цен акций компаний электроэнергетического сектора
ПИК(С4) строительной отрасли, Ростелеком(Т2) и VimpelCom(Т5) - телекоммуникация и ТГК-1(Э8) из электроэнергетического сектора. При этом
компании Магнит(ПТ5) и VimpelCom(Т5) являются лидерами по капитализации в своих отраслях за 2014 год, а Магнит(ПТ5) и Группа компаний
ПИК(С4) за 2015 год по данным [38], что демонстрирует положительный
результат политики этих фирм.
24
а) 1 квартал 2015 год
б) 2 квартал 2015 год
в) 3 квартал 2015 год
г) 4 квартал 2015 год
Рис. 10: Относительная динамика цен акций компаний электроэнергетического сектора
3.3.
Взаимодействие рядов курсов акций при различных
временных лагах
В данном параграфе решается задача поиска временных лагов между
курсами акций. Рассматривается период с 2008 по 2015 года.
Для анализа были выбраны ежедневные цены таких товаров как нефть
сорта Brent и бензин. Выяснялась зависимость стоимости бензина от стоимости нефти. Брался месячный период цен нефти и месячный период цен
бензина, который в последствии сдвигался на различный промежуток времени.
Алгоритм реализован в среде Матлаб. Средний объем выборки 25,
25
исходя из этого выбираем параметр N=2. В приложении представлены таблицы отклонений от комонотонной копула-функции (Таблица 42-48).
Рис. 11: Динамика цен нефти сорта Brent
В таблице 2 подсчитано количество минимальных значений отклонений от комонотонной копула-функции при различных лагах. Из нее следует,
что наибольшую комонотонность цены нефти и бензина показывают при
отсутствии лага. Таким образом, можем сделать вывод, что цена бензина
реагирует на изменение цены на нефть практически мгновенно. Интервью
с главой ФАС - Артемьевым И.Ю. [39] подтверждает этот результат для
возрастающей динамики цен нефти: "Когда цена на нефть росла, лага никакого не было - в тот же день бензин дорожал...". В нашем исследовании
мы получили аналогичный результат и для убывающей динамики. Однако
в 2012 и 2014 годах количество минимальных значений, соответствующих
отсутствию лага, уменьшается. Как мы видим из Рис.11 в эти годы наблюдается сильное падение цены на нефть. Отсюда можно сделать вывод, что
26
отсутствие лага более характерно для возрастающей динамики цен, нежели
для убывающей.
Таблица 2: Количество минимальных значений отклонений от комонотонной копулы при различных лагах
Лаг
0
1
7
14
30
60
120
240
299
2008
10
3
0
1
1
2
0
1
2
2009
12
3
1
0
1
0
0
2
0
2010
10
2
1
0
2
2
0
1
1
2011
11
2
1
0
1
1
1
3
1
27
2012
6
5
1
2
0
1
0
3
3
2013
11
2
0
0
0
2
2
0
0
2014
9
3
0
4
0
2
1
2
2
Сумма:
69
20
4
7
5
10
4
12
9
Выводы
В ходе проделанной работы были сделаны следующие выводы:
1. Российский рынок характеризуется комонотонной зависимостью, что происходит из-за большой концентрации событий, произошедших за рассмотренный период времени и оказывающих негативное влияние на состояние
фондового рынка.
2. Российский рынок 2014 года характеризуется комонотонной связью, но в
2015 году становится выражен независимый характер.
3. Наиболее независимой отраслью показал себя сектор строительства и
недвижимости, в то время как на нефтяном рынке явно выражено комонотонное взаимодействие курсов акций.
4. Компании, политика которых позволяет держать индивидуальную динамику стоимостей акций: Магнит(ПТ5), Группа компаний ПИК(С4), Ростелеком(Т2) и VimpelCom(Т5), ТГК-1(Э8)
5. Между ценами на нефть и бензин отсутствует временной лаг.
6. Возрастающей и убывающей динамике цен на нефть соответствуют разные изменения цены бензина.
28
Заключение
В данной работе реализован метод копула-функции для анализа российского финансового рынка и его секторов. Использованный подход позволил указать некоторые особенности динамики курсов акций, в частности выявлено различие этой динамики в докризисный и кризисный период,
также выделены компании наиболее устойчивые к изменениям в рыночной
среде.
Помимо этого алгоритм применен к задаче поиска временного лага
между ценами нефти и бензина. Было выявлено отсутствие временного лага, а также различие реакции динамики стоимости бензина в зависимости
от положительного или отрицательного изменения цены на нефть.
29
Список литературы
[1] Фантаццини Д. Моделирование многомерных распределений с использованием копула-функций. l. Прикладная эконометрика, (№ 2(22)),
2011.
[2] W. Hoeffding.
Scale-invariant correlation theory.
Schriften des
Mathematischen Seminars und des Instituts fur Angewandte Mathematik
der Universitat, 1940. 5, p.181-183.
[3] A. Sklar. Fonctions de repartition a n dimensions et leurs marges. Publ.
Inst. Statist. Univ. Paris, 1959. 8, p.229–231.
[4] Nelsen B. R. An intriduction ti copulas. Lection notes in statistics, 2-end
edition. New York: Springer. U., 1999.
[5] Благовещенский Ю.Н. Основные элементы теории копул. Прикладная
эконометрика, (№ 2(26)), 2012.
[6] Пеникас Г.И. Модели "копула"в приложении к задачам финансов. Журнал новой экономической ассоциации, (№ 7), 2010.
[7] Фантаццини Д. Моделирование многомерных распределений с использованием копула-функций. lll. Прикладная эконометрика, (№ 4(24)),
2011.
[8] Коршунова Г.Я. Анализ валютного портфеля организации при помощи
копула-функций. Technical report, ФГБОУ ВПО "Рязанский государственный радиотехнический университет".
[9] Соложенцев Е.Д. Алексеев В.В., Шоколов В.В. Логико-вероятностное
30
моделирование портфеля ценных бумаг с использованием копул. Эконометрика, (№ 3(7)), 2006.
[10] Пеникас Г.И. Иерархические копулы в моделировании рисков инвестиционного портфеля. Прикладная эконометрика, (№ 35(3)), 2014.
[11] Калимуллина Л.Д. Бронштейн Е.М. Использование копула-функций
при формировании портфелей ценных бумаг. Риски финансовых рынков, (№ 3(27)), 2011.
[12] Симакова В.Б. Пеникас Г.И. Управление процентным риском на основе
копулы-garch моделей. Прикладная эконометрика, page №1(13), 2009.
[13] Травкин А.И. Конструкции из парных копул в задаче формирования
портфеля акций. Прикладная эконометрика, (№ 32(4)), 2013.
[14] Фантаццини Д. Моделирование многомерных распределений с использованием копула-функций. ll. Прикладная эконометрика, (№ 3(23)),
2011.
[15] Фантаццини Д. Эконометрический анализ финансовых данных в задачах управления риском. Прикладная эконометрика, (№ 2(10)), 2008.
[16] Бологов Я.В. Оценка риска кредитного портфеля с использованием
копула-функций. Прикладная эконометрика, (№ 29(1)), 2013.
[17] Бронштейн Е.М., Прокудина Е.И., Герасимова А.С., and Дубинская К.Г.
Оценка взаимосвязей временных рядов курсов акций. Прикладная эконометрика, (№ 2(22)), 2011.
[18] Сафарян И.А. Бродский Б.Е., Пеникас Г.И. Обнаружение структурных
сдвигов в моделях копул. Прикладная эконометрика, (№ 4(16)), 2009.
31
[19] Зинурова А.Н. Бронштейн Е.М. Копулы специального вида и их применение при анализе состояния финансового рынка. Прикладная эконометрика, (№ 3(27)), 2012.
[20] Кузнецов М. П. Построение интегрального индикатора в ранговых
шкалах с использованием копул для анализа совместного распределения критериев. Машинное обучение и анализ данных, (№ 1(4)), 2012.
[21] Крылов В.А. Моделирование и классификация многоканальных дистанционных изображений с использованием копул. Инофрматика и
ее применения, том 4(№ 4):с. 33–37, 2010.
[22] Кириллов И. А. Кириллов А. А. О массовом составе при ультравысоких энергиях по данным оже-коллаборации (статистический подход с
использованием копулы). Известия РАН. Серия физическая, том 79(№
3):с. 371–373, 2015.
[23] Журов А.Н. Модели оптимального распределения капитала страховой компании. PhD thesis, ФГОБУВПО "Финансовый университет при
правительстве Российской Федерации 2012.
[24] Barry K. Goodwin Jieyuan Zhao.
Dynamic cross-hedge ratios: An
application of copula models. Association’s 2012 AAEA Annual Meeting,
Seattle, Washingto, August 12-14 2012.
[25] Xiaodong Ming, Wei Xu, Ying Li, Juan Du, Baoyin Liu, and Peijun Shi.
Quantitative multi-hazard risk assessment with vulnerability surface and
hazard joint return period. Stoch Environ Res Risk Assess, (№ 29):p. 35–44,
2015.
[26] Поскакухин В.Н. Проблемы внедрения сетей сотовой подвижной связи
32
третьего поколения в московском регионе. Технологии информационного общества, 2009.
[27] Wolff E. F. Schweizer B. Sur une mesure de dependence pour les variables
aleatoires., 1976. Comptes Rendus de l’Academie des Sciences de Paris.
Ser. A, 283, p. 659 – 661.
[28] Wolff E. Schweizer B. On non-parametric measures of dependence for
random variables., 1981. Annals of Statistics, 9, p. 879 – 885.
[29] Sornette D. Malevergne Y. Extreme financial risks (from dependence to
risk management). Springer: Heidelberg., 2006.
[30] Olsen E. T. Darsow W. F., Nguyen B. Copulas and markov process. Illinois
Journal of Mathematics, (№36 (4), 600–642.):p. 600–642., 1992.
[31] www.finam.ru.
[32] Посохов Ю. Е. Трегуб А. Я. Российский фондовый рынок: 2007. События и факты. НАУФОР.
[33] Посохов Ю. Е. Трегуб А. Я. Российский фондовый рынок: 2011. События и факты. НАУФОР.
[34] Посохов Ю. Е. Трегуб А. Я. Российский фондовый рынок: 2012. События и факты. НАУФОР.
[35] Посохов Ю. Е. Трегуб А. Я. Российский фондовый рынок: 2013. События и факты. НАУФОР.
[36] Посохов Ю. Е. Трегуб А. Я. Российский фондовый рынок: 2014. События и факты. НАУФОР.
33
[37] www.finanz.ru.
[38] www.riarating.ru.
[39] www.ria.ru.
34
Приложение
Реализация алгоритма в среде Matlab.
% N-шаг алгоритма; q-количество анализируемых компаний; n-количество
% выборки; T-матрица, нечетные столбцы которой содержат данные о коти% ровках, а четные далее содержат ранг соответствующий порядковой ста% тистике; Y-матрица для ранжирования;
1.
2.
3.
4.
5.
N=...;
q=...;
n=...;
T=zeros(n,2*q);
Y=zeros(n,1);
% Заполняем матрицу Y числами от 1 до n;
6.
7.
8.
for j=1:n
Y(j)=j;
end
% Ранжирование: упорядочиваем данные нечетных столбцов матрицы T в
% матрицу S; в ячейку I(i) записываем индексы j элементов S(i), соответ% ствующих элементам T(j); записываем в четные столбцы T(i+1) ранг эле% ментов T(i);
9. for l=0:q-1
10.
[S, I] = sort(T(:,2*l+1));
11.
T(I(:),2*l+2)=Y(:);
12. end
% C1-матрица копула-функций, соответствующих независимости;
% C_min-матрица копула-функций, соответствующих контрмонотонности;
% C_max-матрица копула-функций, соответствующих комонотонности;
35
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
C1=zeros(N,N);
C_min=zeros(N,N);
C_max=zeros(N,N);
for i=1:N
for j=1:N
C1(i,j)=i*j/N^2;
C_min(i,j)=max(i/N+j/N-1,0);
C_max(i,j)=min(i/N,j/N);
end
end
% Оценка копула-функции временных рядов: h-матрица оценок копул;
% h1-матрица числителей оценок копула-функции;
% строки 28-32, 33-37: поиск соответствующих порядковых статистик;
% строки 38-42: вычисление числителя оценки копула-функции рядов;
% d1-отклонения копул от эталонных независимого характера;
% d_min-отклонения копул от эталонных контрмонотонного характера;
% d_max-отклонения копул от эталонных комонотонного характера;
% строки 50-56: вычисление d1, d_min, d_max;
23. for s=1:2:2*(q-1)-1
24.
for p=s+2:2:2*q-1
25.
h=zeros(N,N);
26.
h1=zeros(N,N);
27.
for i=1:N
28.
for j=1:N
29.
k1=1;
30.
while (T(k1,s+1)~=fix(i*n/N))
31.
k1=k1+1;
32.
end
33.
ps_s=T(k1,s);
34.
k2=1;
35.
while (T(k2,p+1)~=fix(j*n/N))
36
36.
37.
38.
39.
40.
41.
42.
43.
44.
45.
46.
47.
48.
49.
50.
51.
52.
53.
54.
55.
56.
57.
end
58. end
k2=k2+1;
end
ps_p=T(k2,p);
for k=1:n
if ((T(k,s)<=ps_s)&&(T(k,p)<=ps_p))
h1(i,j)=h(i,j)1+1;
end
end
h(i,j)=h1(i,j)/n;
end
end
d1=0;
d_min=0;
d_max=0;
for i=1:N
for j=1:N
d1=d1+abs(h(i,j)-C1(i,j));
d_min=d_min+abs(h(i,j)-C_min(i,j));
d_max=d_max+abs(h(i,j)-C_max(i,j));
end
end
37
В Таблицах 3-41 представлены результаты анализа секторов российского финансового рынка (параграф 3.2). В них наглядно показан характер
взаимодействия между компаниями одного экономического сектора в различные периоды времени.
Таблица 3: Сектор потребительских товаров, 1 квартал 2014 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
+
ПТ3
+
+
ПТ4
+
0
+
ПТ5
+
0
+
+
Таблица 4: Сектор потребительских товаров, 2 квартал 2014 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
+
ПТ3
+
+
ПТ4
+
+
+
ПТ5
+
+
+
+
Таблица 5: Сектор потребительских товаров, 3 квартал 2014 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
0
ПТ3
0
+
38
ПТ4
0
+
0
ПТ5
0
0
0
0
Таблица 6: Сектор потребительских товаров, 4 квартал 2014 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
0
ПТ3
+
0
ПТ4
+
0
+
ПТ5
+
0
+
0
Таблица 7: Сектор потребительских товаров, 1 квартал 2015 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
+
ПТ3
+
+
ПТ4
+
+
+
ПТ5
0
0
0
0
Таблица 8: Сектор потребительских товаров, 2 квартал 2015 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
0
ПТ3
0
0
ПТ4
0
+
0
ПТ5
+
0
0
0
Таблица 9: Сектор потребительских товаров, 3 квартал 2015 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
+
ПТ3
0
0
39
ПТ4
+
+
-
ПТ5
0
0
0
0
Таблица 10: Сектор потребительских товаров, 4 квартал 2015 года
ПТ1
ПТ1
ПТ2
ПТ3
ПТ4
ПТ2
0
ПТ3
0
ПТ4
0
0
0
ПТ5
0
0
0
0
Таблица 11: Сектор строительства и недвижимости, 1 квартал 2014 года
С1
С1
С2
С3
С4
С2
0
С3
+
0
С4
0
+
-
С5
+
0
+
0
Таблица 12: Сектор строительства и недвижимости, 2 квартал 2014 года
С1
C1
C2
C3
C4
С2
0
С3
+
0
С4
0
0
0
С5
+
0
Таблица 13: Сектор строительства и недвижимости, 3 квартал 2014 года
С1
C1
C2
C3
C4
С2
0
С3
0
-
40
С4
+
+
-
С5
+
0
0
0
Таблица 14: Сектор строительства и недвижимости, 4 квартал 2014 года
С1
C1
C2
C3
C4
С2
+
С3
+
+
С4
-
С5
+
+
+
-
Таблица 15: Сектор строительства и недвижимости, 1 квартал 2015 года
С1
C1
C2
C3
C4
С2
+
С3
+
+
С4
0
0
0
С5
0
0
Таблица 16: Сектор строительства и недвижимости, 2 квартал 2015 года
С1
C1
C2
C3
C4
С2
+
С3
0
-
С4
0
0
0
С5
+
+
0
+
Таблица 17: Сектор строительства и недвижимости, 3 квартал 2015 года
С1
C1
C2
C3
C4
С2
-
С3
+
-
41
С4
+
-
С5
+
+
0
Таблица 18: Сектор строительства и недвижимости, 4 квартал 2015 года
С1
C1
C2
C3
C4
С2
+
С3
+
+
С4
0
0
0
С5
0
0
0
Таблица 19: Сектор нефтедобычи и нефтепереработки, 1 квартал 2014 года
Н1
Н1
Н2
Н3
Н4
Н2
+
Н3
+
+
Н4
+
+
+
Н5
0
+
0
0
Таблица 20: Сектор нефтедобычи и нефтепереработки, 2 квартал 2014 года
Н1
Н1
Н2
Н3
Н4
Н2
+
Н3
+
+
Н4
+
+
+
Н5
+
+
+
+
Таблица 21: Сектор нефтедобычи и нефтепереработки, 3 квартал 2014 года
Н1
Н1
Н2
Н3
Н4
Н2
+
Н3
+
+
42
Н4
0
+
+
Н5
0
+
+
0
Таблица 22: Сектор нефтедобычи и нефтепереработки, 4 квартал 2014 года
Н1
Н1
Н2
Н3
Н4
Н2
0
Н3
+
0
Н4
+
0
+
Н5
+
0
+
+
Таблица 23: Сектор нефтедобычи и нефтепереработки, 1 квартал 2015 года
Н1
Н1
Н2
Н3
Н4
Н2
+
Н3
+
+
Н4
0
+
+
Н5
+
+
+
+
Таблица 24: Сектор нефтедобычи и нефтепереработки, 2 квартал 2015 года
Н1
Н1
Н2
Н3
Н4
Н2
0
Н3
+
+
Н4
+
+
+
Н5
0
+
+
+
Таблица 25: Сектор нефтедобычи и нефтепереработки, 3 квартал 2015 года
Н1
Н1
Н2
Н3
Н4
Н2
0
Н3
0
0
43
Н4
0
0
+
Н5
0
0
+
+
Таблица 26: Сектор нефтедобычи и нефтепереработки, 4 квартал 2015 года
Н1
Н1
Н2
Н3
Н4
Н2
0
Н3
+
+
Н4
+
0
0
Н5
+
0
+
+
Таблица 27: Сектор телекоммуникации, 1 квартал 2014 года
Т1
Т1
Т2
Т3
Т4
Т2
+
Т3
0
+
Т4
+
+
+
Т5
+
+
+
+
Таблица 28: Сектор телекоммуникации, 2 квартал 2014 года
Т1
Т1
Т2
Т3
Т4
Т2
+
Т3
+
+
Т4
+
+
+
Т5
0
0
0
0
Таблица 29: Сектор телекоммуникации, 3 квартал 2014 года
Т1
Т2
Т1
Т2
Т3
Т4
44
Т3
Т4
Т5
0
0
+
0
+
+
Таблица 30: Сектор телекоммуникации, 1 квартал 2015 года
Т1
Т1
Т2
Т3
Т4
Т2
0
Т3
0
+
Т4
0
0
0
Т5
0
0
0
+
Таблица 31: Сектор телекоммуникации, 2 квартал 2015 года
Т1
Т1
Т2
Т3
Т4
Т2
-
Т3
+
-
Т4
+
+
Т5
0
0
+
Таблица 32: Сектор телекоммуникации, 3 квартал 2015 года
Т1
Т1
Т2
Т3
Т4
Т2
+
Т3
0
0
Т4
0
0
0
Т5
0
0
0
0
Таблица 33: Сектор телекоммуникации, 4 квартал 2015 года
Т1
Т1
Т2
Т3
Т4
Т2
0
Т3
Т4
0
0
Т5
0
0
0
45
Таблица 34: Сектор электроэнергетики, 1 квартал 2014 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
0
Э3
+
+
Э4
+
+
0
Э5
+
+
+
+
Э6
+
+
+
+
+
Э7
+
0
+
+
+
+
Э8
0
0
0
0
0
+
0
Таблица 35: Сектор электроэнергетики, 2 квартал 2014 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
0
Э3
0
+
Э4
0
+
+
Э5
0
+
+
+
Э6
0
+
+
+
+
Э7
0
+
+
+
+
+
Э8
+
0
0
0
-
Таблица 36: Сектор электроэнергетики, 3 квартал 2014 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
+
Э3
+
0
Э4
0
0
0
46
Э5
0
0
0
0
Э6
0
0
0
0
0
Э7
+
+
+
0
0
0
Э8
0
0
0
+
+
0
0
Таблица 37: Сектор электроэнергетики, 4 квартал 2014 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
+
Э3
0
+
Э4
0
0
+
Э5
+
+
+
0
Э6
+
+
+
+
+
Э7
+
+
+
+
+
+
Э8
+
+
+
+
+
+
+
Таблица 38: Сектор электроэнергетики, 1 квартал 2015 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
+
Э3
+
0
Э4
+
+
0
Э5
+
0
+
+
Э6
+
0
+
0
+
Э7
+
+
+
+
+
+
Э8
0
0
0
0
0
0
0
Таблица 39: Сектор электроэнергетики, 2 квартал 2015 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
0
Э3
0
0
Э4
0
0
0
47
Э5
Э6
0
0
0
0
Э7
0
0
+
0
Э8
0
0
0
0
0
0
-
Таблица 40: Сектор электроэнергетики, 3 квартал 2015 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
0
Э3
0
0
Э4
0
+
0
Э5
0
0
0
Э6
0
0
0
+
Э7
0
0
0
0
+
+
Э8
0
0
0
0
+
+
+
Таблица 41: Сектор электроэнергетики, 4 квартал 2015 года
Э1
Э1
Э2
Э3
Э4
Э5
Э6
Э7
Э2
0
Э3
0
0
Э4
0
0
+
Э5
0
0
0
0
Э6
0
+
0
0
0
Э7
0
0
0
0
+
0
Э8
0
0
0
0
0
0
0
Ниже представлены таблицы показывающие величину отклонения,
соответствующую комонотонному характеру взаимодействия между ценами нефти и бензина при различных лагах сопоставленные разным месяцам
рассматриваемого года. Цветом выделены минимальные значения за месяц.
Выводы по Таблицам 42-48 представлены в параграфе взаимодействие рядов курсов акций при различных временных лагах (параграф 3.3).
48
49
0.35
0.35
0.20
0.25
0.30
60
120
240
299
0.25
7
30
0.15
1
0.30
0.05
0
14
1
Лаг
0.20
0.10
0.45
0.10
0.05
0.15
0.10
0.05
0.05
2
0.30
0.30
0.25
0.30
0.15
0.30
0.25
0.15
0.05
3
0.45
0.23
0.50
0.41
0.18
0.00
0.18
0.05
0.00
4
0.12
0.12
0.50
0.45
0.21
0.21
0.21
0.12
0.12
5
0.40
0.26
0.26
0.36
0.31
0.31
0.40
0.31
0.12
6
0.09
0.09
0.36
0.00
0.09
0.05
0.09
0.05
0.00
7
0.40
0.40
0.31
0.26
0.26
0.21
0.36
0.26
0.17
8
0.26
0.40
0.31
0.21
0.31
0.36
0.40
0.21
0.26
9
0.54
0.37
0.28
0.46
0.20
0.11
0.11
0.11
0.07
10
Таблица 42: Отклонение от копулы комонотонной зависимости, 2008 год
0.13
0.39
0.45
0.29
0.39
0.24
0.24
0.18
0.18
11
0.27
0.23
0.18
0.32
0.32
0.27
0.36
0.14
0.05
12
50
0.30
0.25
0.35
0.35
0.35
60
120
240
299
0.35
7
30
0.25
1
0.25
0.20
0
14
1
Лаг
0.45
0.29
0.29
0.45
0.29
0.34
0.45
0.24
0.13
2
0.36
0.00
0.36
0.23
0.05
0.27
0.14
0.05
0.00
3
0.21
0.45
0.50
0.40
0.40
0.45
0.45
0.21
0.17
4
0.15
0.25
0.40
0.30
0.45
0.10
0.05
0.05
0.05
5
0.32
0.32
0.23
0.36
0.23
0.27
0.27
0.18
0.14
6
0.18
0.41
0.23
0.09
0.41
0.23
0.14
0.05
0.05
7
0.26
0.36
0.36
0.26
0.26
0.26
0.26
0.26
0.21
8
0.26
0.40
0.36
0.36
0.21
0.36
0.36
0.31
0.21
9
0.23
0.05
0.27
0.32
0.27
0.32
0.18
0.05
0.05
10
Таблица 43: Отклонение от копулы комонотонной зависимости, 2009 год
0.20
0.20
0.20
0.25
0.25
0.20
0.30
0.20
0.10
11
0.18
0.23
0.27
0.41
0.23
0.45
0.27
0.09
0.05
12
51
0.50
0.29
0.39
0.29
0.50
60
120
240
299
0.18
7
30
0.13
1
0.24
0.08
0
14
1
Лаг
0.10
0.10
0.35
0.35
0.05
0.15
0.10
0.05
0.00
2
0.28
0.33
0.37
0.28
0.33
0.28
0.28
0.24
0.15
3
0.40
0.26
0.26
0.31
0.21
0.40
0.26
0.26
0.21
4
0.25
0.40
0.40
0.00
0.00
0.35
0.20
0.00
0.00
5
0.14
0.05
0.18
0.18
0.09
0.36
0.32
0.09
0.00
6
0.45
0.40
0.21
0.26
0.40
0.36
0.17
0.26
0.17
7
0.27
0.14
0.27
0.36
0.36
0.32
0.09
0.09
0.05
8
0.36
0.21
0.26
0.21
0.36
0.26
0.31
0.26
0.12
9
0.34
0.18
0.26
0.34
0.38
0.42
0.34
0.22
0.18
10
Таблица 44: Отклонение от копулы комонотонной зависимости, 2010 год
0.08
0.12
0.35
0.35
0.38
0.38
0.35
0.27
0.23
11
0.26
0.18
0.46
0.14
0.34
0.22
0.18
0.14
0.18
12
52
0.17
0.17
0.29
0.38
0.17
60
120
240
299
0.29
7
30
0.13
1
0.13
0.08
0
14
1
Лаг
0.09
0.05
0.09
0.32
0.09
0.18
0.05
0.09
0.05
2
0.27
0.23
0.15
0.31
0.23
0.23
0.23
0.19
0.15
3
0.41
0.24
0.41
0.24
0.33
0.37
0.28
0.20
0.11
4
0.27
0.27
0.27
0.23
0.15
0.27
0.27
0.15
0.15
5
0.12
0.08
0.12
0.31
0.23
0.46
0.27
0.08
0.08
6
0.22
0.38
0.22
0.42
0.22
0.38
0.30
0.22
0.18
7
0.24
0.35
0.20
0.24
0.35
0.43
0.35
0.17
0.13
8
0.26
0.46
0.46
0.42
0.18
0.42
0.22
0.14
0.10
9
0.35
0.12
0.04
0.12
0.38
0.38
0.38
0.19
0.15
10
Таблица 45: Отклонение от копулы комонотонной зависимости, 2011 год
0.38
0.15
0.12
0.38
0.42
0.31
0.31
0.12
0.08
11
0.20
0.33
0.24
0.41
0.41
0.41
0.41
0.28
0.20
12
53
0.38
0.21
0.25
0.17
0.21
60
120
240
299
0.38
7
30
0.38
1
0.42
0.38
0
14
1
Лаг
0.08
0.33
0.04
0.42
0.17
0.08
0.04
0.00
0.04
2
0.23
0.23
0.23
0.19
0.23
0.15
0.15
0.15
0.19
3
0.14
0.50
0.30
0.26
0.22
0.14
0.26
0.14
0.14
4
0.35
0.38
0.15
0.35
0.12
0.08
0.12
0.08
0.04
5
0.34
0.14
0.30
0.34
0.34
0.46
0.26
0.14
0.18
6
0.04
0.12
0.42
0.23
0.15
0.15
0.19
0.08
0.08
7
0.39
0.31
0.28
0.39
0.43
0.46
0.24
0.13
0.09
8
0.08
0.42
0.38
0.17
0.13
0.25
0.21
0.17
0.13
9
0.39
0.13
0.24
0.39
0.35
0.35
0.24
0.17
0.13
10
Таблица 46: Отклонение от копулы комонотонной зависимости, 2012 год
0.18
0.54
0.54
0.10
0.14
0.54
0.34
0.10
0.10
11
0.30
0.30
0.34
0.22
0.18
0.38
0.18
0.14
0.10
12
54
0.26
0.50
0.26
0.42
0.38
60
120
240
299
0.10
7
30
0.10
1
0.18
0.06
0
14
1
Лаг
0.41
0.33
0.41
0.41
0.28
0.37
0.33
0.33
0.24
2
0.30
0.30
0.30
0.18
0.30
0.26
0.30
0.30
0.22
3
0.42
0.27
0.08
0.19
0.15
0.42
0.27
0.12
0.08
4
0.31
0.23
0.27
0.19
0.35
0.31
0.27
0.15
0.12
5
0.38
0.42
0.42
0.18
0.22
0.46
0.34
0.18
0.10
6
0.35
0.38
0.19
0.31
0.15
0.35
0.31
0.12
0.12
7
0.31
0.23
0.38
0.38
0.46
0.50
0.27
0.19
0.19
8
0.38
0.10
0.50
0.42
0.14
0.26
0.14
0.10
0.06
9
0.31
0.31
0.20
0.13
0.28
0.35
0.31
0.20
0.13
10
Таблица 47: Отклонение от копулы комонотонной зависимости, 2013 год
0.46
0.42
0.04
0.17
0.13
0.29
0.13
0.08
0.04
11
0.27
0.31
0.15
0.19
0.23
0.27
0.38
0.12
0.08
12
55
0.21
0.25
0.17
0.33
0.25
60
120
240
299
0.25
7
30
0.21
1
0.17
0.17
0
14
1
Лаг
0.24
0.46
0.46
0.41
0.37
0.20
0.28
0.20
0.15
2
0.19
0.15
0.19
0.23
0.38
0.31
0.19
0.19
0.12
3
0.10
0.10
0.46
0.46
0.22
0.46
0.34
0.14
0.10
4
0.21
0.21
0.50
0.46
0.21
0.29
0.25
0.08
0.04
5
0.15
0.15
0.38
0.35
0.42
0.46
0.27
0.04
0.04
6
0.22
0.38
0.22
0.30
0.22
0.26
0.22
0.26
0.18
7
0.22
0.34
0.22
0.22
0.38
0.14
0.30
0.30
0.26
8
0.23
0.15
0.42
0.12
0.15
0.12
0.19
0.31
0.31
9
0.08
0.08
0.12
0.08
0.12
0.08
0.15
0.08
0.08
10
Таблица 48: Отклонение от копулы комонотонной зависимости, 2014 год
0.15
0.20
0.41
0.50
0.20
0.15
0.20
0.15
0.11
11
0.15
0.23
0.31
0.46
0.38
0.08
0.08
0.04
0.08
12
Отзывы:
Авторизуйтесь, чтобы оставить отзыв