САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ
ФАКУЛЬТЕТ ПРИКЛАДНОЙ МАТЕМАТИКИ – ПРОЦЕССОВ УПРАВЛЕНИЯ
УДК 519.816
На правах рукописи
Ледовская Вероника Александровна
Выпускная квалификационная работа
Исследование мер риска в теории принятия решений
МК.3002.2014 «Прикладная математика и процессы управления»
Направление подготовки: 01.06.01 «Математика и механика»
Научный руководитель,
доктор физ.-мат. наук,
профессор
Колбин В. В.
Рецензент,
кандидат физ.-мат. наук,
директор ЦПИ ФГАОУ ВО «СПбПУ»
Козырев С. В.
Санкт-Петербург
2018
СОДЕРЖАНИЕ
Введение .................................................................................................................. 3
Постановка задачи .................................................................................................. 5
Обзор литературы .................................................................................................. 6
Глава 1. Анализ, оценка и управление риском ................................................... 8
1.1. Понятие риска ............................................................................................ 8
1.2. Методы управления риском ..................................................................... 9
1.3. Анализ мер риска на когерентность ...................................................... 10
1.4. Практическое применение ..................................................................... 24
Глава 2. Принцип ожидаемой полезности ......................................................... 30
2.1. Понятие ожидаемой полезности............................................................. 30
2.2. Функция полезности фон Неймана – Моргенштерна .......................... 30
2.3. Модель ожидаемой полезности ............................................................. 33
2.4. Метод ожидаемой полезности и меры риска ....................................... 34
2.5. Связь степеней рисковости и производных функции полезности...... 38
Глава 3. Концепция стохастического доминирования ..................................... 44
3.1. Принцип стохастического доминирования .......................................... 44
3.2. Стохастическое доминирование и меры риска .................................... 47
3.3. Стохастическое доминирование и степени рисковости ..................... 51
Глава 4. Сравнение методов ................................................................................ 53
Глава 5. Оценка интенсивности предпочтения ................................................. 66
5.1. Некоторые свойства интенсивности предпочтения ............................ 66
5.2. Степени и меры интенсивности предпочтения .................................... 71
5.3. Порог принятия решения ....................................................................... 83
5.4. Порог принятия решения на практике .................................................. 84
Результаты ............................................................................................................. 89
Заключение ........................................................................................................... 91
Список литературы .............................................................................................. 92
Приложение .......................................................................................................... 96
2
Введение
Людям каждый день приходится принимать решения, опираясь на свой
жизненный опыт и знания. При принятии решений, касающихся социальной,
экономической, политической жизни общества важным аспектом является
количественная и качественная оценка предпочтений лиц принимающих решения
(ЛПР), а также оценка риска различных альтернатив, из которых ЛПР стремятся
выбрать наилучшую. Жизнь человека с каждым годом становится все более
насыщенной событиями и информацией, темп жизни, особенно в больших
городах, невероятно высок. Поэтому ЛПР также стремится к тому, чтобы
принятие решения было быстрым и своевременным, а определение наилучшей
альтернативы было возможно для любых ситуаций.
Конечно, нет универсального средства, позволяющего решить любую
жизненную задачу, но ученые издавна стремятся к его созданию и рассматривают
ситуацию с различных точек зрения. Поэтому и существуют математические
направления
посвященные
исследованию
риска,
полезности,
отношений
предпочтения и т.д. Сложность задач, решаемых в данных исследованиях, состоит
в том, что они работают с эмпирическими данными, а, значит, не всегда есть
возможность применить одну теорию ко всем задачам.
Целью научной работы является исследование нескольких концепций
принятия
решений:
аппарата
мер
риска,
принципа
стохастического
доминирования, метода ожидаемой полезности и оценки интенсивности
предпочтений, а так же выявление и описание случаев, в которых применение
перечисленных методов будет наиболее эффективным.
Для достижения цели необходимо решить следующие задачи:
•
исследование мер риска на когерентность, выбор наиболее удобных мер
риска для применения к эмпирическим данным;
•
исследование таких подходов к принятию решения, как принцип
стохастического доминирования и метод ожидаемой полезности;
3
•
описание взаимосвязи понятия степени рисковости с методами
стохастического доминирования и ожидаемой полезности;
•
исследование аппарата для оценки интенсивности предпочтений;
•
разработка новых способов оценки интенсивности предпочтений.
4
Постановка задачи
Целью научной работы является изучение нескольких концепций принятия
решений, а также исследование способа оценки интенсивности предпочтений с
помощью функции полезности. Исследование разделено на несколько этапов.
В первой главе проведено исследование аппарата мер риска, проверка
наиболее часто используемых мер риска на когерентность. Результатом данного
исследования является отбор мер риска для дальнейшего изучения в соотношении
с другими методами ранжирования альтернатив.
Вторая глава посвящена методу ожидаемой полезности и условиям,
накладываемым
мерами
риска на функцию полезности. Кроме того в ней
изучена взаимосвязь степеней рисковости и порядков неприятия риска,
характеризующихся производными функции полезности.
В третьей главе проведено исследование принципа стохастического
доминирования, его взаимосвязи с методом ожидаемой полезности, аппаратом
мер риска и степенями рисковости.
Четвертая глава содержит сравнение трех методов в условиях рисковости
разных степеней на численном примере.
В пятой главе изучены методы оценки интенсивности предпочтения и их
свойства, проведена работа по предложению новых характеристик и критериев
оценивания интенсивности предпочтений на основании функции полезности.
5
Обзор литературы
Теория принятия решений имеет большой спектр различных методов для
принятия решений в условиях риска и неопределенности. Они применяются во
многих областях жизнедеятельности, проводя обработку эмпирических данных.
Сложность принятия решения состоит не только формализации и описании
предпочтений лица принимающего решения (ЛПР), но и в оценке риска.
Пожалуй, одним из основополагающих понятий теории принятия решений
является понятие полезности. Впервые оно было использовано Бернулли [1] в
XVIII, при исследовании им Петербургского парадокса. В дальнейшем, Дж. Фон
Нейман и О. Моргенштерн, которые считаются родоначальниками теории
полезности, сформулировали ряд аксиом, характеризующих рациональное
поведение ЛПР [2]. Кроме того, предложили функцию полезности, которая
помогает отразить индивидуальные предпочтения лица принимающего решения.
Фишберн позднее описал системы аксиом, необходимые для применения метода
ожидаемой полезности [3]. А Дж. Пратт и К. Эрроу охарактеризовали с помощью
функции полезности отношение ЛПР к риску. Развитие теории полезности
продолжается и по сей день, но ее методы не всегда удобны из-за необходимости
задавать функцию полезности индивидуума, так как данная задача бывает весьма
непростой для эмпирических данных.
Альтернативным методом для сравнения решений является принцип
стохастического доминирования [4, 5, 6]. При его использовании нет
необходимости в обращении к явной функции полезности, хотя на нее
накладываются некоторые ограничения [7]. Таким образом, данный метод может
быть полезным в тех случаях, когда ЛПР сложно определить функцию
полезности. Взаимосвязь различных порядков стохастического доминирования и
порядков риска рассматривается в статьях [8, 9]. А в работе [4] описывается связь
степени рисковости с принципом стохастического доминирования.
В теории принятия решений ЛПР может опираться не только на полезность
альтернатив, но и на величину риска потерь или доходов при выполнении
6
альтернативы. Удобным способом количественного описания риска является
построение меры риска. Одни из первых результатов в этом направлении
принадлежат Фишберну [10]. Позднее Ф. Арцер, Ф. Делбаен, Ж.-М. Эбер и Д.
Хиф ввели понятие когерентной меры риска [11] ввиду отсутствия какой-либо
универсальной меры. За рубежом множество работ посвящено исследованию
свойств мер риска [12, 13, 14, 15, 16] и возможности их применения в различных
сферах жизнедеятельности. В некоторых статьях приводятся исследования,
посвященные взаимосвязи порядков стохастического доминирования и мер риска
[9, 10].
Все описанные выше методы позволяют ЛПР упорядочить альтернативы по
величине их полезности или риска. Но в некоторых ситуациях бывает
недостаточно просто упорядочить альтернативы. Необходимо еще и узнать
насколько предпочтительнее для ЛПР одна альтернатива по сравнению с другой.
Поэтому
исследователями
был
предложен
метод
оценки
интенсивности
предпочтения [17, 18, 19], позволяющий выбрать решение в задачах, где нет
единственной наилучшей альтернативы.
7
Глава 1. Анализ, оценка и управление риском
1.1. Понятие риска
Происхождение
понятия
«риск»
неизвестно,
существует
множество
предположений о его возникновении. Некоторые ученые считают, что слово
имеет арабское происхождение, другие предполагают, что оно происходит от
греческих слов ridsikon , ridsa – утес, скала. В Европе слово использовалось еще в
средневековье, но в различных предметных областях, кроме того, было мало
распространено. Благодаря появлению книгопечатания слово получает широкое
распространение,
особенно
в
Испании
и
Италии.
Поэтому
некоторые
исследователи предполагают, что термин «риск» имеет испанское (risco –
отвесная скала) или итальянское (risiko – опасность, угроза; risicare – лавировать
между скал) происхождение [20]. В русский язык слово пришло в XVIII – XIX
веке из французского (risqué – опасность).
Существует большое количество интерпретаций понятия «риск», условно
можно выделить две ветви: экономическая и математическая.
К экономической ветви можно отнести
определение риска, как
возможности получения неблагоприятного результат (потерь) или прибыли. В
данной ветви риск может рассматриваться как негативный фактор, влияние
которого ЛПР стремится снизить. Но при этом риск может быть интерпретирован
и как положительный фактор – возможность получения награды за риск,
прибыли.
Понятие риска как события, которое может произойти в условиях
неопределенности с некоторой вероятностью и приводит к различным исходам,
можно отнести к математической ветви интерпретации понятия «риск». В таком
смысле риск определяется как вероятность отклонения фактического значения от
ожидаемого результата. Соответственно, чем больше это отклонение, тем выше
риск [21].
Благодаря американскому экономисту Фрэнку Найту [22] фактор риска
получил наиболее интенсивное развитие как один из важнейших факторов,
8
влияющих на принятие решений. Найт был первым, кто предложил понимать
риск как объективную вероятность события, которая может быть выражена
количественно, в частности, в виде математически вероятностного распределения
доходов. Тогда с ростом вероятности стандартного отклонения от ожидаемой
величины риск уменьшается, иначе - увеличивается.
На практике возникновение рисковых ситуаций объясняется рядом
факторов: отсутствием полной информации о том или ином событии,
возникновением непредвиденных ранее ситуаций, которые значительно влияют
на решение проблемы и т. д. Для того, чтобы сделать правильный и наиболее
выгодный
выбор
необходимо
правильно
оценить
риск
и
суметь
его
минимизировать. Инструментом для оценки и управления риском служит такой
математический аппарат, как теория риска. Она неразрывно связана с теорией
вероятности, истории развития этих наук тесно связаны друг с другом.
В данной работе будем понимать риск как случайное событие, возникающее
с некоторой вероятностью, тогда для количественного описания риска могут быть
использованы различные меры риска, которые мы рассмотрим далее.
1.2. Методы управления риском
Защита от негативного влияния рисков заключается в прогнозировании
уровня
дополнительных
затрат,
оценке
тяжести
возможного
ущерба,
использовании всего финансового механизма для ликвидации угрозы риска или
его последствий.
Кроме того, общеизвестны четыре общих метода управления риском:
избегание, предотвращение потерь и контроль, снижение степени риска,
поглощение.
1. Избегание заключается в отказе от совершения рискового мероприятия.
Но для финансового предпринимательства исключение риска обычно исключает и
прибыль.
2. Предотвращение потерь и контроль как метод управления финансовым
риском означает определенный набор превентивных и последующих действий,
9
которые обусловлены необходимостью предотвратить негативные последствия,
уберечься от случайностей, контролировать их размер, если потери уже имеют
место или неизбежны.
3. Снижение степени риска. Если невозможно или нежелательно избежать
риска, то можно уменьшить влияние его негативных последствий на работу
компании.
4. Поглощение состоит в признании ущерба и отказе от его страхования. К
поглощению прибегают, когда сумма предполагаемого ущерба незначительна и
ей можно пренебречь.
1.3. Анализ мер риска на когерентность
Введем вероятностное пространство (Ω, А, Р), где Ω – множество
элементарных исходов, Α - σ-алгебра, заданная на Ω, а P – вероятностная мера,
определенная на множествах из Α.
Случайной величиной X (или риском) будем называть произвольное
измеримое отображение из Ω в R.
Обозначим
χ
совокупность
всех
случайных
величин
на
(Ω, А).
Перенумеровав элементы Ω некоторым произвольным образом: Ω = {ω1 ,K , ω n } ,
обозначим X (ω i ) = x i , i = 1,…n и будем отождествлять случайные величины
X ∈ χ с векторами X = ( x1 ,K , x n ) из R n .
Напомним [23], что отношением предпочтения ≼ на множестве
называется произвольное полное транзитивное бинарное отношение, а мерой
риска называется произвольный функционал :
представляет отношение предпочтения ≼, если
( )≤ ( )⟹
→ . Говорят, что мера риска
≼ , для ,
∈ .
Не существует универсальной меры риска, которая могла бы одинаково
хорошо охарактеризовать риск для различных ситуаций. Поэтому в 1999 году
группа ученых выделила ряд свойств, которыми должна обладать мера риска, для
проведения качественного исследования рисков и управления ими [11]. Опишем
10
выделенные учеными свойства и, т.к. на практике меры риска чаще используются
для оценки финансовых рисков, дадим пояснение значения свойств на их
примере.
1. Монотонность
≤
⟹ ( ) ≤ ( ), для ,
∈
Свойство монотонности позволяет описать рациональные предпочтения
ЛПР. Например, при описании ситуации, в которой ЛПР стремится получить
доход, данное свойство определяет правило «чем больше, тем лучше».
Соответственно, если доходность X не меньше доходности Y, то согласно
свойству риск альтернативы с доходностью X
не больше риска, при выборе
альтернативы с доходностью Y. Данное свойство описывает отношения
предпочтения и критерии выбора в условиях риска для рационального ЛПР.
2. Супераддитивность
( + ) ≤ ( ) + ( ), для ,
∈ .
Свойство супераддитивности характеризует поведение меры риска в
ситуациях дробления или, наоборот, агрегирования (объединения) рисков и
говорит о том, что совокупный риск нескольких альтернатив не превышает сумму
рисков
альтернатив.
Обладающая
данным
свойством
мера
отражает
преимущества диверсификации (инвестирование в разные активы, с целью
снижения рисков) и подходит для описания условия: при объединении рисков
ресурсы должны высвобождаться.
3. Положительная однородность
(
)=
( ),
> 0,
∈ .
Данное свойство описывает необходимость независимости меры риска от
выбора денежных единиц, в которых измеряется доходность (инвестиции), т.е.
была инвариантна по отношению к ним. Величина риска не должна зависеть от
единиц измерения — будь то рубли или тысячи рублей. Таким образом, свойство
однородности отражает рост риска пропорционально росту инвестиций.
11
4. Инвариантность относительно сдвига
( +
)= ( )+ ,
∈ ,
∈ .
здесь I = (1,1,K ,1) – вектор с единичными координатами длины n.
Описанное выше свойство отражает следующую ситуацию: при увеличении
доходности X на известную ЛПР константу a, риск для данной альтернативы
уменьшается.
Когерентной мерой риска называется функционал :
→ , обладающий
свойствами монотонности, супераддитивности, положительной однородности и
инвариантности относительно сдвигов [24].
Рассмотрим наиболее часто используемые меры рисков, и проведем анализ
данных мер рисков на когерентность.
1. Математическое ожидание
Непрерывное распределение µ =
∞
∫ x f ( x)dx,
−∞
n
Дискретное распределение µ = ∑ x i p i .
i =1
• Трансляционная инвариантность
Так как математическое ожидание постоянной a равно этой постоянной, т.е.
постоянную a можно рассматривать как случайную величину x, которая может
принимать
только
одно
значение a c
вероятностью
равной
единице.
Поэтому µ = a ⋅ 1 = a . Поэтому
n
µ ( x + a ) = ∑ xi p i + a ⋅ 1 = µ ( x ) + a
i =1
То
есть
свойство
трансляционной
математического ожидания.
12
инвариантности
выполняется
для
• Супераддитивность.
Рассмотрим две случайные величины X и Y (докажем для дискретного
случая)
µ ( X + Y ) = ∑ ( xi + y j ) p ij = ∑ xi pij + ∑ y j p ij = ∑ xi ∑ p ij + ∑ y j ∑ p ij =
ij
ij
ij
ij
ij
ij
ij
= ∑ xi pi• + ∑ y j p• j = ∑ xi p ( X = xi ) + ∑ y j p (Y = y j ) = µ ( X ) + µ (Y ) .
i
j
ij
ij
Таким образом, математическое ожидание – это супераддитивная функция.
• Положительная однородность.
Докажем, что это условие выполняется для математического ожидания,
рассмотрим случайную величину X и постоянную с ≥ 0 :
n
n
i =1
i =1
µ (cx) = ∑ cxi p i = c ∑ xi pi = cµ ( x) .
• Монотонность.
Рассмотрим две случайные величины: X и Y, причем Y = X + ε , где ε ≥ 0 случайная величина, то есть X ≤ Y
µ (Y ) = µ ( X + ε ) = ∑ ( xi + ε j ) p ij = ∑ xi p ij + ∑ ε j p ij = µ ( X ) + µ (ε ) ,
ij
ij
ij
т.е. математическое ожидание обладает свойством монотонности.
2. Дисперсия, стандартное отклонение
Непрерывное распределение σ
2
=
∞
∫ (x − µ)
−∞
n
2
2
Дискретное распределение σ = ∑ xi pi
i =1
• Трансляционная инвариантность
13
2
f ( x)dx
Так как математическое ожидание постоянной a равно этой постоянной, т.е.
постоянную a можно рассматривать как случайную величину x, которая может
принимать только одно значение a с вероятностью равной единице, то:
σ ( x + a) 2 = µ [( x + a − µ ( x + a )) 2 ] = µ [( x − µ ( x) + a − µ (a )) 2 ] =
= µ [( x − µ ( x)) 2 ] = σ ( x) 2
Рассуждения для дисперсии справедливы и для стандартного отклонения:
σ ( x + a) = σ ( x + a) 2 = σ ( x)
То есть свойство трансляционной инвариантности не выполняется для дисперсии
и стандартного отклонения.
• Супераддитивность.
Рассмотрим дисперсию суммы двух случайных величин X и Y (докажем
для дискретного случая):
σ ( X + Y ) 2 = σ ( X ) 2 + σ (Y ) 2 + 2 Cov ( X , Y ) ,
где ковариация случайных величин имеет следующий вид:
Cov ( X , Y ) = σ ( X ) ⋅ σ (Y ) ⋅ r ,
где r – коэффициент корреляции. Так как r ≤ 1 , то верно неравенство
Cov ( X , Y ) ≤ σ ( X ) ⋅ σ (Y )
воспользовавшись которым, получаем
σ ( X + Y ) 2 ≤ σ ( X ) 2 + σ (Y ) 2 + 2 σ ( X ) σ (Y ) ⇒
⇒ σ ( X + Y ) 2 ≤ (σ ( X ) + σ (Y )) 2
Последнее неравенство эквивалентно неравенству
σ ( X + Y ) ≤ σ ( X ) + σ (Y ) .
Таким образом, стандартное отклонение – это супераддитивная функция. Но
дисперсия таковой не является, так как для нее свойство супераддитивности
выполняется тогда и только тогда, когда r ≤ 0 , что неверно в общем случае.
14
• Положительная однородность.
Рассмотрим для дисперсии и стандартного отклонения случайной величины
X и постоянной с ≥ 0 :
σ (сx) 2 = µ[(cx − µ (cx)) 2 ] = µ[(c ( x − µ ( x)) 2 ] = c 2 µ[ x − µ ( x)] = c 2σ ( x) 2
σ (сx) = c 2σ ( x) 2 = cσ ( x)
Следовательно, стандартное отклонение обладает свойством положительной
однородности. Но для дисперсии это свойство выполняется только в случае, если
с = 0и с = 1.
• Монотонность.
Рассмотрим две случайные величины: X и Y, причем Y = X + ε , где ε ≥ 0 –
случайная величина, то есть X ≤ Y , тогда
σ ( y ) 2 = σ ( x + ε ) 2 = σ ( x ) 2 + σ (ε ) 2 + 2 Cov ( x + ε )
Величина σ ( x) 2 меньше σ ( y ) 2 , тогда и только тогда, когда
1
2
σ (ε ) 2 + 2 Cov ( x + ε ) < 0 ⇒ Cov ( x + ε ) < − σ (ε ) 2 .
Однако, так как свойство монотонности должно выполняться для произвольных
положительных случайных величин ε , то дисперсия не удовлетворяет этому
свойству, так же как и стандартное отклонение.
3. Мера риска Market-risk
β x = rX , R
σX
σR
Приведем выражение к более удобному для исследования виду:
β x = rX , R
σ X Cov ( X , R) σ X Cov ( X , R)
=
=
σR
σ Xσ R σ R
σ R2
• Трансляционная инвариантность.
Рассмотрим случайную величину X и постоянную c > 0 :
15
(1.3.1)
β ( X + c) =
Cov ( X + c, R )
σ
2
R
µ[( X + c − µ ( X + c)) ⋅ ( R − µ ( R))] Cov ( X , R)
=
= β x ( x)
σ R2
σ R2
=
Таким образом, для данной меры не выполняется свойство инвариантности.
• Супераддитивность.
Докажем, что выполняется условие супераддитивности:
Cov ( X + Y , R )
β (X + Y ) =
σ
=
=
2
R
=
µ [( X + Y − µ ( X + Y )) ⋅ ( R − µ ( R))]
=
σ R2
µ [( X + Y − µ ( X ) − µ (Y )) ⋅ ( R − µ ( R))]
=
σ R2
µ [( X − µ ( X ))( R − µ ( R) + (Y − µ (Y )) ⋅ ( R − µ ( R))]
=
σ R2
=
µ [( X − µ ( X )) ⋅ ( R − µ ( R)] + µ [(Y − µ (Y )) ⋅ ( R − µ ( R))]
=
σ R2
=
µ [( X − µ ( X )) ⋅ ( R − µ ( R)] µ [(Y − µ (Y )) ⋅ ( R − µ ( R))]
+
=
σ R2
σ R2
=
Cov ( X , R)
σ
2
R
+
Cov (Y , R )
σ R2
= β ( X ) + β (Y ) .
• Положительная однородность.
Для случайной величины X и постоянной c > 0 выполняется
β (cX ) =
Cov (cX , R)
σ
2
R
=
µ [(cX − µ (cX )) ⋅ ( R − µ ( R))] c ⋅ Cov ( X , R)
=
= c ⋅ β (X )
σ R2
σ R2
• Монотонность.
Рассмотрим две случайные величины: X и Y, причем Y = X + ε , где ε ≥ 0 –
случайная величина, то есть X ≤ Y
β(X + ε) =
=
Cov ( X + ε , R)
σ
2
R
=
µ [( X + ε − µ ( X + ε )) ⋅ ( R − µ ( R))]
=
σ R2
µ [( X − µ ( X )) ⋅ ( R − µ ( R)] µ [(ε − µ (ε )) ⋅ ( R − µ ( R))]
+
=
σ R2
σ R2
16
=
Cov ( X , R)
σ
2
R
+
Cov (ε , R)
σ R2
= β ( X ) + β (ε ) .
4. Мера риска Lower Partial Moment
Непрерывное распределение
LPM n,t ( X ) =
t
∫ (t − x)
n
f ( x)dx ,
(1.3.2)
−∞
Дискретное распределение
LPM n,t ( X ) = ∑ max[(t − xi ),0]n pi ,
(1.3.3)
i
В статье [12] приводятся следующие результаты проверки меры на
когерентность.
• Положительная однородность
Мера LPM n ,t ( X ) удовлетворяет условию положительной однородности
тогда и только тогда, когда t = 0 и n = 1 . Более того, если t = 0 , то порядок
однородности LPM n ,t ( X ) равен n, то есть
LPM n,0 (λX ) = λn LPM n,0 ( X ) .
Доказательство. Покажем положительную однородность меры:
n
t
LPM n,t ( X ) = λ ∑ max − xi I t − x≥0 pi .
i
λ
Мера LPM n ,t ( X ) положительно однородна тогда и только тогда, когда t = 0
. Иначе,
если t > 0 ∧ λ > 1 или t < 0 ∧ 0 < λ < 1 , то
LPM n,t (λX ) < λn LPM n,t ( X ) ;
если t > 0 ∧ 0 < λ < 1 или t < 0 ∧ λ > 1 , то
LPM n,t (λX ) > λn LPM n,t ( X ) .
17
Поэтому LPM n ,t ( X ) положительно однородна тогда и только тогда, когда
t = 0 и n = 1.
• Трансляционная инвариантность
Мера
не
LPM n ,t ( X )
инвариантности.
Она
удовлетворяет
уменьшается
условию
(увеличивается),
трансляционной
если
постоянная
(безрисковая) составляющая добавляется к (вычитается из) случайной величине.
Доказательство. Рассмотрим случайную величину Y с yi = ε , ∀i ∈ Ω
ε ∈ R \ {0} и величину Z = X + Y , тогда
LPM n,t ( Z ) ≤ ∑ [(t − xi ) I t − xi ≥0 − ε I t − xi −ε ≥0 ]n pi .
i∈Ω
• Монотонность
Мера LPM n ,t ( X ) монотонна при всех t ∈ R и n > 0 .
Доказательство. Рассмотрим случайную величину Y с yi ≥ 0 , ∀i ∈ Ω и
∃i : yi > 0 . Пусть случайная величина Z определяется как X + Y. Тогда
LPM n,t ( Z ) ≤ ∑ [(t − xi ) I t − xi ≥0 − yi I t − xi − yi ≥0 ]n pi <
i∈Ω
< ∑ [(t − xi ) I t − xi ≥0 ]n pi = LPM n ,t ( X ) .
i∈Ω
• Супераддитивность
Мера LPM n ,t ( X ) супераддитивна при всех t ∈ R+ и 0 < n ≤ 1 .
Доказательство.
Перепишем
формулу
(1.3.2)
с
использованием
индикаторной функции, тогда
LPM n,t ( X ) = ∑ max[(t − xi ) I t − x≥0 ]n pi .
i
Рассмотрим две случайные величины X и Y. Рассмотрим единичную
ситуацию:
18
[(t − x − y ) I t − x− y≥0 ]n p ≤ [(2t − x − y ) I t − x − y≥0 ]n p =
= [(t − x ) I t − x ≥ y + (t − y ) I t − y ≥ x ] n p ≤ [(t − x ) I t − x ≥ 0 + (t − y ) I t − y ≥ 0 ] n p .
(1.3.4.)
Первое неравенство следует из того, что величина t неотрицательна.
Последнее неравенство выполняется, так как верны следующие отношения
(t − x ) I t − x ≥ 0 ≥ (t − x ) I t − x ≥ y и (t − y ) I t − y ≥ 0 ≥ (t − y ) I t − y ≥ x . В первом отношении правая
и левая части неравенства равны всегда, кроме случая, когда I t − x ≥ 0 ≠ I t − x ≥ y . Когда
0 < t − x < y , то (t − x ) I t − x ≥ 0 > 0 = (t − x ) I t − x ≥ y . В случае выполнения неравенства
y ≤ t − x < 0 , получаем (t − x ) I t − x ≥ 0 = 0 > (t − x ) I t − x ≥ y .
Мы рассмотрели одно состояние, для того, чтобы обобщить полученный
вывод на множество состояний, используется лемма 1.3.1:
n > 1,
Лемма 1.3.1. Для
n = 1 , [0 < n < 1] и
ai , bi ∈ R+ ,
∀i ∈ Ω , и
∃ i ∈ Ω : ai > 0 ∧ bi > 0 , выполняется неравенство (равенство):
∑ (ai + bi ) n > (=)[<]∑ (ai ) n + ∑ (bi ) n .
i∈Ω
i∈Ω
i∈Ω
Тогда для n ≥ 1
∑ [(t − xi − yi ) I t − x − y ≥0 ]n pi ≤
i
i∈Ω
i
≤ ∑ [(t − xi ) I t − xi ≥0 + (t − yi ) I t − yi ≥0 ]n pi ≤
(1.3.5)
i∈Ω
≤ ∑ [(t − xi ) I t − xi ≥0 ] n pi + ∑ [(t − yi ) I t − yi ≥0 ]n pi .
i∈Ω
Когда
t ≥ 1,
а
i∈Ω
n < 1,
то
знак
неравенства
(1.3.5)
изменится
на
противоположный. А при условии t < 0 в неравенстве (1.3.4) знак изменится, так
как мы суммируем отрицательные числа. Поэтому, в общем случае, мера не
является супераддитивной.
5. Мера рискаValue-at- risk
VaRα ( X ) = inf{l ∈ R : p ( X ≤ l ) ≥ α } или VaRα ( X ) = − F −1 (α ) .
19
• Монотонность:
Если X и Y – случайные величины, для которых выполняется неравенство:
X ≤ Y, тогда p ( X ≥ l ) ≤ p (Y ≥ l ) для любого l ∈ R . Таким образом, если l
принадлежит множеству {l ∈ R : p (Y ≤ l ) ≥ α } , тогда l принадлежит множеству
{l ∈ R : p ( X ≤ l ) ≥ α} . Отсюда, в частности
inf{l ∈ R : p ( X ≤ l ) ≥ α } ≤ inf{l ∈ R : p (Y ≤ l ) ≥ α } ,
поэтому
VaRα ( X ) ≤ VaRα (Y ) .
• Положительная однородность:
Выполнение данного условия рассмотрено в [13], приведем его ниже.
Для случайной величины X и постоянной c > 0 выполняется
VaRα (ñX ) = ñVaRα ( X ) .
Доказательство:
VaRα (ñX ) = inf{ x ∈ R : p (ñX ≤ x) ≥ α } =
= c ⋅ inf{x c ∈ R : p ( X ≤ x c) ≥ α } = c VaRα ( X ) .
• Трансляционная инвариантность:
Для случайной величины X и постоянной b > 0 выполняется
VaRα ( X + b) = VaRα ( X ) + b .
Доказательство:
VaRα ( X + b) = inf{ x ∈ R : p ( X + b ≤ x) ≥ α } = inf{ x ∈ R : p ( X ≤ x − b) ≥ α } =
= inf{s ∈ R : p ( X ≤ s ) ≥ α } + b = VaRα ( X ) + b ,
где s = x − b .
• Супераддитивность:
В общем случае VaRα (X ) не является супераддитивной мерой, данное
обстоятельство описано в [14, 15]. Приведем для доказательства этого
утверждения контрпример.
20
1
Представим VaRα (X ) в следующем виде: VaRα ( X ) = (c (1 − α )) z . Это
простейшее определение VaRα (X ) , оно верно, если квантиль xα
задается
значением
1 − α = p ( X > xα ) = cxα− z .
Тогда:
VaRα ( X + Y ) ≈ (2c (1 − p))
1
z
1
= 2 z (c (1 − p))
1
z
1
VaRα ( X ) + VaRα (Y ) ≈ (2c (1 − p)) z .
Если z < 1 , то 2
1
z
> 2 и, следовательно, VaRα ( X ) не удовлетворяет
условию супераддитивности. Например, для с = 1 и z = 1 2 мы получаем
VaRα ( X + Y ) ≈ 2(VaRα ( X ) + VaRα (Y ) ,
а для ñ = 1 и z = 1 4 мы имеем
VaRα ( X + Y ) ≈ 8(VaRα ( X ) + VaRα (Y ) .
6. Мера риска Expected Shortfall
X с уровнем доверия α ∈ (0 ;1) мера риска
Для случайной величины
ESα ( X ) будет иметь следующий вид:
Непрерывное распределение
ESα ( X ) =
1
α
α
∫ VaRβ ( X , β ) dβ ,
(1.3.6)
0
Дискретное распределение
ESα ( X ) = −α −1 ( E [ X ⋅ 1{ X ≤ xα } ] + xα (α − p ( X ≤ xα ))) ,
(1.3.7)
где xα = inf{ x ∈ R : p ( X ≤ x ) ≥ α } ,
1, если α удовлетворяет условию X ≤ xα
1{ X ≤ x } =
– индикаторная
α
0, если α не удовлетворяет условию X ≤ xα
функция, причем E [ 1{ X ≤ xα } ] = α . Такое представление описано в статье [25].
21
• Положительная однородность.
Для случайной величины X и постоянной c > 0 выполняется
ES α (ñX ) = ñESα ( X ) .
Доказательство следует из положительной однородности VaRα (X ) по формуле
(1.3.6):
1
1
1
c
ESα (ñX ) =
VaRβ (cX , β ) dβ =
VaRβ ( X , β ) dβ = cESα ( X ) .
∫
1−α α
1 − α α∫
• Трансляционная инвариантность.
Для случайной величины X и постоянной b > 0 выполняется
ESα ( X + b) = ESα ( X ) + b .
Доказательство следует из трансляционной инвариантности VaRα ( X ) .
• Монотонность.
Рассмотрим две случайные величины X и Y. Если X ≤ Y , то
ES α ( X ) ≤ ES α (Y ) .
Доказательство следует из монотонности меры риска VaRα ( X ) .
• Супераддитивность меры.
Доказательство выполнения этого свойства приведено в [25].
Рассмотрим представление меры риска ES α ( X ) в виде (1.3.7):
Введем обозначение:
1{ X
1{(αX) ≤ x} =
1{ X
≤ x}
≤ x}
, если p ( X = x ) = 0
+
α − p( X ≤ x)
p ( X ≤ x)
⋅1
{ X = x}
можем заметить, что
1{(αX)≤ x} ∈ [0 ;1] ,
22
, если p ( X = x ) > 0 ,
E [1{(αX)≤ x} ] = α ,
α −1 E [ X ⋅ 1{(αX)≤ x} ] = − ESα ( X ) .
Тогда
ESα ( X + Y ) = −α −1 E [( X + Y ) ⋅ 1{(αX)+Y ≤( x+ y) } ] =
α
= −α −1 E [ X ⋅ 1{(αX)+Y ≤( x + y )
α}
+ Y ⋅ 1{(αX)+Y ≤( x+ y )
α}
].
Составим разность ESα ( X ) + ESα (Y ) − ESα ( X + Y ) и определим ее знак:
α ( ES α ( X ) + ES α (Y ) − ES α ( X + Y )) = − E [ X ⋅ 1{(αX)≤ x } ] − E [Y ⋅ 1{(αY )≤ y
α
+ E [ X ⋅ 1{(αX)+Y ≤( x+ y )
α}
α}
+ Y ⋅ 1{(αX)+Y ≤( x+ y )
α}
]+
]=
= E [ X (1{(αX)+Y ≤( x+ y )
(α )
(α )
(α )
− 1{ X ≤ x } ) + Y (1{ X +Y ≤( x + y ) } − 1{Y ≤ y } )] ≥
}
α
α
α
α
≥ xα E [1{(αX)+Y ≤( x+ y )
α}
(α )
(α )
(α )
− 1{ X ≤ x } ] + yα E [1{ X +Y ≤( x + y ) } − 1{Y ≤ y } ] =
α
α
α
≥ xα (α − α ) + yα (α − α ) = 0 .
В неравенстве мы использовали следующий факт:
1{(αX) +Y ≤( x+ y ) } − 1{(αX) ≤ x } ≥ 0, если X > xα
α
α
(α )
(α )
1{ X +Y ≤( x+ y )α } − 1{ X ≤ xα } ≤ 0, если X < xα
Таким образом, получаем, что мера ES α ( X ) супераддитивна.
7. Мера риска Worst Conditional Expectation
Для случайной величины
X с уровнем доверия α ∈ (0 ;1) мера риска
WCEα ( X ) будет иметь следующий вид:
WCEα ( X ) = − inf{E[ X | A] : A ∈ Ω, p( A) > α } .
Впервые эта мера была введена в статье [11] как когерентная мера риска. По
своей сути она является наихудшим математическим ожиданием, поэтому
удовлетворяет всем условиям когерентности [16, 25].
Результаты проведенного анализа мер риска были опубликованы в статьях
[26, 27] и в книге [28].
23
1.4. Практическое применение
Области применения мер риска разнообразны, с помощью мер риска можно
оценивать не только финансовые или материальные потери или прибыли, но и
использовать меры рисков для оценки риска смертельных случаев при
возникновении чрезвычайных ситуаций, а именно стихийных бедствий. В данном
параграфе будут найдены оценки рисков по количеству погибших при
наводнениях, землетрясениях, извержениях вулканов и цунами.
Постановка задачи
По статистическим данным, приведенным в табл. 1. Приложения, о
количестве погибших при наводнениях, землетрясениях, извержениях вулканов и
цунами
• провести оценку рисков в количестве погибших с помощью
рассмотренных ранее мер риска;
• составить шкалы рисковости по каждой мере рисков;
• сравнить полученные результаты.
Решение
Так как в расчетах использованы эмпирические данные [29, 30], то для
расчета мер использовались их дискретные формы записи. Существуют
некоторые особенности вычисления мер для дискретной функции распределения,
которые будут описаны ниже. Расчеты проводились с помощью программы
Microsoft Excel 2007.
1. Расчет мер
Проведем расчеты для мер, проанализированных в данной главе.
• дисперсия и стандартное отклонение
Стихийное бедствие
Дисперсия
Стандартное отклонение
наводнение
землетрясение
извержение вулкана
цунами
1915184872
3898477104
268796024
144777643
43763
62438
16395
12032
24
• Мера β-риск
Для расчета меры риска β-риск используем формулу (1.3.1). Для начала
необходимо рассчитать величину σ R – стандартное отклонение для случайной
величины, характеризующей общее количество жертв для всех стихийных
бедствий, которое получилось следующим σ R = 40810,35.
Промежуточные результаты, которые были получены для нахождения меры
β-риск, приведены в табл. 2 Приложения. В таблице расположенной ниже
приведены значения меры β-риск для разных видов стихийных бедствий.
Стихийное бедствие
наводнение
землетрясение
извержение вулкана
цунами
β-риск
1,03288
1,299354
0,317872
0,292144
• Мера Value-at-risk
Мера риска VaRα (X )
для эмпирических данных представляет собой
выборочный квантиль соответствующего уровня α, результаты для данной меры
риска приведены ниже в таблице.
Стихийное бедствие
наводнение
землетрясение
извержение вулкана
цунами
VaRα (X )
α = 0,05
1525
1000
9
119
α = 0,3
5624
1621
700
2182
α = 0,5
7625
2500
1335
3620
α = 0,95
107485
87587
29025
30000
• Мера Lower Partial Moment
При расчете меры LPM n ,t по формуле (1.3.3) необходимо задать уровень t –
уровень
потерь,
будем
принимать
его
равным
различным
распределения. Тогда для удобства обозначим LPM n ,α –
квантилям
мера LPM n ,t , для
которой величина t равна квантилю с уровнем вероятности α. Уровень α = 0,05
25
нецелесообразно, так как этому уровню соответствует только одно наблюдение, и
при расчете меры получаем значение 0.
LPM n ,α
Стихийное бедствие
α = 0,3
540
306
72
123
наводнение
землетрясение
извержение вулкана
цунами
α = 0,5
1298
895
360
622
α = 0,95
87432
99341
29249
24700
• Мера Expected Shortfall
Для расчета меры ESα ( X ) можно использовать несколько способов, работе
были использованы два из них:
‒ Первый способ, представлен в статье [16],
ESα ( X ) = −
1 [ Nα ]
∑ Xi ,
[ Nα ] i =1
(1.4.1)
‒ второй способ, представлен в статье [17],
ES α ( X ) = −
1
Nα
[ Nα ]
∑ X i + ( Nα − [ Nα ]) X [ Nα ] ,
i =1
(1.4.2)
здесь N – количество исходов, [ Nα ] – целая часть Nα . Для удобства обозначим
ES α(1) ( X ) – меру Expected Shortfall, рассчитанную по первому способу, а ES α( 2 ) ( X ) –
меру, рассчитанную по первому способу. Результаты вычислений приведены
ниже в виде таблицы.
Стихийное бедствие
наводнение
землетрясение
извержение вулкана
цунами
Стихийное бедствие
наводнение
ES α(1) ( X )
α = 0,05
1525
2000
40
119
α = 0,5
4929
1802
664
2273
α = 0,95
13059
17920
4642
9480
α = 0,3
3802
1421
340
2064
ES α( 2 ) ( X )
α = 0,05
1916
26
α = 0,5
5036
α = 0,95
16533
α = 0,3
3829
Стихийное бедствие
землетрясение
извержение вулкана
цунами
ES α( 2 ) ( X )
α = 0,05
1000
19
325
α = 0,5
1658
609
2340
α = 0,95
12189
3741
8861
α = 0,3
1267
280
1762
• Мера Worst Conditional Expectation
Мера
риска
WCEα ( X )
дает
грубые
результаты
в
случае,
когда
рассматриваемые исходы близки к равновероятным, это обусловлено ее
природой. Использование данной меры распределений с равновероятными
событиями нецелесообразно.
2. Составление шкал
Теперь составим шкалы рисковости для рассмотренных мер, на них сверху
вниз в порядке убывания рисковости расположены стихийные бедствия.
Шкала №1 изображенная на рис. 1.1, получена для мер: VaR0,05 ( X ) ,
)
)
)
)
VaR0,95 ( X ) , LPM 1,0, 03 ( X ) , LPM 1,0, 05 ( X ) , ES 0(1, 05
( X ) , ES 0(1,95
( X ) , ES 0( 2, 05
( X ) , ES 0( ,295
(X ) .
Наводнения
Землетрясения
Цунами
Извержения вулканов
Рис. 1.1. Шкала рисковости №1
Шкала №2 изображенная на рис. 1.2, получена для мер: VaR0,3 ( X ) , VaR0,5 ( X )
, ES 0(1,3) ( X ) , ES 0(1,5) ( X ) , ES 0( ,23) ( X ) , ES 0( ,25) ( X ) .
Наводнения
Цунами
Землетрясения
Извержения вулканов
Рис. 1.2. Шкала рисковости №2
27
Шкала №3 изображенная на рис. 1.3, получена для мер LPM 1,0,95 ( X ) и βриск.
Землетрясения
Наводнения
Извержения вулканов
Цунами
Рис. 1.3. Шкала рисковости №3
3. Сравнение результатов
Можно сделать вывод, что самая большая дисперсия (разброс) жертв
характерна для землетрясений, за ними следуют наводнения, извержения и
цунами.
По шкале рисковости VaRα (X ) с уровнем α = 0,05 и α = 0,95 видно, что
самые
большие
потери
ожидаются
при
наводнениях,
затем
следуют
землетрясения, цунами и извержения. Но для шкалы рисковости VaRα (X ) с
уровнем α = 0,3 и α = 0,5 ситуация меняется, в порядке убывания следуют:
наводнения, цунами, землетрясения, извержения. То есть шкала рисковости
VaRα (X ) зависит от того, для какой части распределения рассчитывается мера: для
хвостовой или центральной. Аналогичная ситуация складывается и для меры
ESα (X ) . Так же при сравнении результатов ES α(1) ( X ) и ES α( 2 ) ( X ) видно, что
ES α( 2 ) ( X ) имеет больший запас, чем ES α(1) ( X ) , но все же не такой большой, как
VaRα (X ) . Таким образом, наиболее предпочтительной мерой при расчете рисков
для эмпирических распределений будет ES α( 2 ) ( X ) .
Шкала меры LPM n ,α ( X ) дает абсолютно другие результаты. Для уровня
вероятности α = 0,95 получена следующая шкала: землетрясения, наводнения,
извержения, цунами, а α = 0,3 и α = 0,5 – наводнения, землетрясения, цунами,
извержения, что совпадает со шкалой рисковости ESα (X ) с уровнем α = 0,05 и α =
28
0,95. Мера LPM n ,α ( X ) не является когерентной, как и VaRα (X ) , но дает результаты
менее грубые, чем VaRα (X ) .
Шкала для β-риска совпадает со шкалой LPM n ,α ( X ) для уровня α = 0,95.
Она не имеет аналогов среди когерентных мер риска.
Меру WCEα ( X ) возможно вычислить только для распределений, в которых
события не являются равновероятными. Поэтому использование WCEα ( X ) в
распределениях
близких
к
равновероятным
опубликованы в книге [28].
29
нецелесообразно.
Результаты
Глава 2. Принцип ожидаемой полезности
2.1. Понятие ожидаемой полезности
Впервые математическая формулировка понятия ожидаемой полезности
появилась еще в XVIII веке. В попытке объяснить Санкт-Петербургский парадокс
ученые Габриэль Крамер и Даниэль Бернулли [1] пришли к выводу, что люди при
принятии решения стремятся увеличить ожидаемую ценность (полезность), а не
ожидаемый максимальный выигрыш. В ходе своих рассуждений Бернулли
предложил функцию полезности и показал, что функция ожидаемой полезности
для нее конечна. При этом ученый не пытался измерить полезность с помощью
той функции и не оценивал рациональность поведения людей. Только спустя
почти два века Джон фон Нейман и Оскар Моргенштерн «практически
определили
численную
полезность
как
объект,
для
которого
подсчет
математических ожиданий является законным» (Нейман и Моргенштерн, 1970,
с. 55 в [2]). Они описали рациональное поведение ЛПР с помощью системы из
пяти аксиом, из которых можно вывести принцип максимума ожидаемой
полезности. Этот принцип используется и в настоящее время при принятии
решений в условиях риска и неопределенности.
2.2. Функция полезности фон Неймана – Моргенштерна
Принцип рационального принятия решений в условиях неопределенности с
использованием функции полезности опирается на пять аксиом, определенных
фон Нейманом и Моргенштерном. Они отражают необходимые условия для
непротиворечивого поведения индивидуума. Данные аксиомы описаны в книге
[31] и приведены ниже. Для начала введем определение.
Определение 2.2.1. Обозначим G( x, z : α ) ситуацию, в которой индивид с
вероятностью α получает денежную сумму х, и с вероятностью (1 − α ) – сумму z.
Кроме того, введем следующее обозначение: пусть S – множество всех
возможных исходов.
30
Аксиома 2.2.1.
Аксиома сравнимости (полноты). Для всего множества
возможных исходов S ЛПР может либо предпочесть исход х исходу у ( x f y) ,
либо ( y f x) , либо он безразличен в отношении к выбору между х и у ( x ~ y) .
Запись (x f y) означает, что исход х предпочтительнее исхода у либо индивид
безразличен в отношении к выбору между х и у.
Аксиома 2.2.2. Аксиома транзитивности (состоятельности). Если ( x f y)
и ( y f z) , то ( x f z) . Если ( x ~ y) и ( y ~ z) , то ( x ~ z) .
Аксиома 2.2.3.
Аксиома сильной независимости. Предположим, что
сложилась ситуация, в которой ЛПР с вероятностью α получает денежную сумму
х, и с вероятностью (1 − α ) – сумму z, т. е. G( x, z : α ) . Сильная независимость
означает, что если индивид безразличен в отношении к выбору между х и у
( x ~ y) , то он также будет безразличен в отношении к выбору между ситуациями
G( x, z : α ) и G( y, z : α ) , т. е. из ( x ~ y) следует G( x, z : α ) ~ G( y, z : α ) .
Аксиома 2.2.4. Аксиома измеримости. Если ( x f y ~ z) или ( x ~ y f z) , то
существует единственная вероятность α, такая, что y ~ G( x, z : α ) .
Поясним смысл этой аксиомы на примере. Пусть существуют три исхода:
х = 1000; у = 0; величина z означает смерть игрока. Исходя из здравого смысла,
смерть нельзя сравнивать ни с каким выигрышем, и соответствующего этому
исходу значения вероятности α существовать не может. Однако в жизни бывают
ситуации, когда некий проигрыш равнозначен смерти. Тогда утверждение
y ~ G( x, z : α ) можно считать справедливым для некоторого значения 0 ≤ α ≤ 1 .
Аксиома 2.2.5. Аксиома ранжирования. Если альтернативы у и и находятся
по предпочтительности между альтернативами х и z и можно построить такие
ситуации, в которых ЛПР безразличен в отношении выбора между у и G ( x, z : α1 ) ,
a также при выборе между и и G( x, z : α 2 ) , то при α1 > α 2 ( y f u) .
Утверждение аксиомы вполне соответствует здравому смыслу: чем больше
вероятность крупного выигрыша, тем больше затрат требует игра, т. е. тем
большая плата потребуется за приобретение права участвовать в этой игре.
31
Учитывая то, что люди выбирают по принципу «чем больше, тем лучше», и
выполняются приведенные аксиомы, можно сделать вывод, что ЛПР при
принятии решения будет стремиться к максимизации ожидаемой полезности.
Другими словами, из всех возможных решений он выберет то, которое
обеспечивает наибольшую ожидаемую полезность. Сформулируем определение
полезности по Нейману-Моргенштерну, так же описанное в книге [31].
Определение 2.2.2. Полезность – это некоторое число, приписываемое
лицом, принимающим решение, каждому возможному исходу. Функция
полезности Неймана-Моргенштерна для ЛПР показывает полезность, которую он
приписывает каждому возможному исходу. У каждого ЛПР своя функция
полезности, которая показывает его предпочтение к тем или иным исходам в
зависимости от его отношения к риску.
Определение 2.2.3. Пусть на некотором множестве альтернатив X
определено отношение предпочтения ≽. Тогда вещественнозначная функция ( )
определенная на X будет называться функцией полезности, если выполняется
условие: ≽ ⟺ ( ) ≥ ( ) для всех ,
Функция
полезности
не
только
∈ .
помогает
количественно
описать
предпочтения ЛПР, но и оценить его отношение к риску. ЛПР по отношению к
риску условно делятся на три группы: склонные к риску, несклонные и
безразличные. В общем случае функция полезности может быть трех типов,
графики которых можно охарактеризовать так:
•
несклонность к риску – строго вогнутая функция, у которой каждая дуга
кривой лежит выше своей хорды;
•
безразличие к риску – прямая линия;
•
склонность к риску – строго выпуклая функция, у которой каждая дуга
кривой лежит ниже своей хорды.
Независимо друг от друга ученые Джон Пратт (1964) [32] и Кеннет Эрроу
(1971) [33] разработали меру неприятия риска, которая представляет собой
отношение второй производной функции полезности к первой, взято с
32
отрицательным знаком. Эту меру принято называть коэффициентом ЭрроуПратта неприятия риска:
Rα = −
u ′′( x)
.
u ′( x)
Мера удобна тем, что является инвариантной относительно линейных
преобразований и, кроме того, принимает постоянное значение для линейных и
экспоненциальных функций полезности.
Прежде, чем переходить к вероятностному описанию принципа ожидаемой
полезности, дадим общее определение ожидаемой полезности события.
Определение 2.2.4. Ожидаемая полезность события равна сумме
произведений вероятностей исходов на значения полезностей этих исходов.
2.3. Модель ожидаемой полезности
Принятие решений в условиях риска, где риск выражается в вероятностной
неопределенности
исходов,
подразумевает
выбор
альтернатив,
которые
описываются с помощью вероятностной меры на множестве исходов, например, с
помощью
функции
распределения.
Опишем
соотношение
нескольких
альтернатив, используя соответствующие функции распределения, а также
функции полезности для данных альтернатив.
Введем вероятностное пространство (Ω, А, Р), где Ω – множество
элементарных исходов, Α – σ-алгебра, заданная на Ω, а P – вероятностная мера,
определенная на множествах из Α.
Случайной величиной X (или риском) будем называть произвольное
измеримое отображение из Ω в R.
Обозначим
χ
совокупность
всех
случайных
величин
на
Перенумеровав элементы Ω некоторым произвольным образом: Ω = {
обозначим
(
#)
=
X ∈ χ с векторами
#,
(Ω, А).
,…,
! },
$ = 1, … , & и будем отождествлять случайные величины
= ( ,…,
!)
из
!
.
33
Будем
рассматривать
отношение
альтернатив, которые являются
предпочтения
вероятностными
≻
на
множестве
мерами. В работе
[3],
разработано несколько систем аксиом, таких, что для любого удовлетворяющего
им отношения ≻ на P существует такая функция ∈ (, что
) ≺ + ⟺ ,( , )) < ,( , +),
(2.3.1)
где для любых ), + ∈ . через ,( , )) и ,( , +) обозначено математическое
ожидание функции u относительно вероятностных мер ), +.
В теории ожидаемой полезности могут быть установлены аксиоматические
характеристики отношений ≻ и ≥ на P, которые обладают односторонними
представлениями вида:
или
) ≺ + ⟹ ,( , )) < ,( , +),
) ≈ + ⇒ ,( , )) = ,( , +),
) ≼ + ⟹ ,( , )) ≤ ,( , +),
) ≺ + ⟹ ,( , )) < ,( , +).
В книге [3] модели ожидаемой полезности посвящено несколько глав. В
данной работе мы не будем приводить подробное описание всех моделей, а
воспользуемся моделью (2.3.1) и рассмотрим ее взаимосвязь с мерами риска.
2.4. Метод ожидаемой полезности и меры риска
Совместное использование мер риска и ожидаемой полезности накладывает
некоторые ограничения на производные функции полезности. В предыдущих
разделах мы рассмотрели несколько различных мер риска. Выберем из них
наиболее часто используемые на практике, таковыми являются VaRαF ( X ) ,
LPM n ,t ( X ) ,
ES αF ( X ) , и проанализируем их связь с методом ожидаемой
полезности.
Для начала ознакомимся с теоремой, опубликованной в статье [34].
Теорема 2.4.1. Пусть F и G являются кумулятивными распределениями
двух различных случайных величин с квантилями QF и
34
QG соответственно.
Тогда F доминирует G в терминах стохастического доминирования второго
порядка тогда и только тогда, когда
α
α
0
0
∫ QF (t )dt ≥ ∫ QG (t )dt
для всех значений α в диапазоне 0 ≤ 1 ≤ 1, и имеет место строгое неравенство
при некотором значении α.
Так же воспользуемся леммой 2.4.1 приведенной А. А. Новоселовым в его
книге [23].
Лемма 2.4.1. Математическое ожидание от функции E (h( X )) от случайной
величины, имеющей функцию распределения F, может быть вычислено по
формуле
∞
1
−∞
0
−1
∫ h ( x) dF ( x) ≥∫ h ( F (v))dv .
Теперь докажем предложение.
Теорема 2.4.1. Пусть F (x) и G(x) являются функциями распределения двух
различных случайных величин. И пусть ES αF ( X ) – мера Expected Shortfall, а U ( x) F
– неубывающая (невозрастающая) функция полезности, соответствующие
случайной величине с функцией распределения F (x) . Тогда для них выполняется:
ES αF ( X ) ≤ ES αG ( X )
⇔
E[U ( x )] F ≥ E[U ( x )]G .
Доказательство. Используем представление Expected Shortfall (1.3.7) и
распишем представление меры через квантиль распределения величины X
ES α ( X ) =
1
α
α
∫ VaR β ( X , β ) dβ = −
0
1
α
F
α∫
−1
0
Тогда ES αF ( X ) ≤ ES αG ( X ) можно записать так
−
1
α
F
α∫
−1
( β ) dβ ≤ −
0
−
1
1
α
G
α∫
−1
( β ) dβ
0
α
(F
α∫
−1
( β ) − G −1 ( β ))dβ ≤ 0
0
35
⇒
( β ) dβ .
отсюда следует, что F −1 ( β ) − G −1 ( β ) ≥ 0 или F −1 ( β ) ≥ G −1 ( β ) , следовательно,
выполняется соотношение F ( x) ≤ G ( x) . Тогда, используя лемму 2.4.1 и учитывая,
что функция полезности монотонная, получим
∞
1
−∞
0
E[U ( x)] F = ∫ U ( x)dF ( x) = ∫ U ( F −1 ( β ))dβ ≥
1
∞
−1
≥ ∫ U (G ( β ))dβ = ∫ U ( x)dG ( x) = E[U ( x)]G .
−∞
0
Таким образом, получили E[U ( x)] F ≥ E[U ( x)]G . Обратное утверждение
доказывается аналогично. ■
Теорема 2.4.2. Пусть F (x) и G (x) являются функциями распределения двух
различных случайных величин. И пусть LPM nF,t ( X ) – мера Lower Partial Moment, а
U ( x ) F – неубывающая (невозрастающая) функция полезности, соответствующие
случайной величине с функцией распределения
F (x) .
Тогда для них
выполняется:
LPM nF,t ( X ) ≤ LPM nG,t ( X ) ⇔ E[U ( x)] F ≥ E[U ( x)]G .
Доказательство. Используем формулу записи (1.3.2) для меры LPM n,t ( X ) и
распишем ее представление через квантиль распределения величины X
n
LPM n ,t ( X ) = E{max[(t − x) ,0]} =
t
∫ (t − x)
n
f ( x)dx .
−∞
Тогда LPM nF,t ( X ) ≤ LPM nG,t ( X ) можно записать так
t
∫ (t − x )
t
n
f ( x )dx ≤
−∞
t
∫ (t − x)
∫ (t − x)
n
g ( x )dx ⇒
−∞
n
( f ( x) − g ( x))dx ≤ 0
−∞
тогда верно выражение
справедливо
неравенство
f ( x) − g ( x) ≤ 0
(t − x ) n ≥ 0 .
36
или F ( x) ≤ G ( x) , так как при ∀t
Отсюда
следует,
что
выполняется
соотношение F −1 ( β ) ≥ G −1 ( β ) . Тогда, используя лемму 2.4.1 и учитывая, что
функция полезности монотонная, получим следующее выражение
∞
1
−∞
0
E[U ( x)] F = ∫ U ( x)dF ( x) = ∫ U ( F −1 ( β ))dβ ≥
1
∞
0
−∞
≥ ∫ U (G −1 ( β ))dβ = ∫ U ( x)dG ( x) = E[U ( x)]G .
Таким образом, получили E[U ( x)] F ≥ E[U ( x)]G . Обратное утверждение
доказывается аналогично. ■
Теорема 2.4.3. Пусть F (x) и G (x) являются функциями распределения двух
различных случайных величин. И пусть VaRαF ( X ) – мера Value-at-risk, а U ( x) F –
неубывающая
(невозрастающая)
функция
полезности,
соответствующие
случайной величине с функцией распределения F (x) . Тогда для них выполняется:
VaRαF ( X ) ≤ VaRαG ( X )
⇔
E[U ( x)] F ≥ E[U ( x )]G .
Доказательство. Используем представление VaRα (X ) для непрерывного
распределения и запишем представление меры через квантиль распределения
величины X
VaRα ( X ) = − F −1 (α ) .
Тогда VaRαF ( X ) ≤ VaRαG ( X ) можно записать так F −1 ( β ) ≥ G −1 ( β ) . Отсюда
следует, что F ( x) ≤ G ( x) . Тогда, используя лемму 2.4.1 и учитывая, что функция
полезности монотонная, получим
∞
1
−∞
0
E[U ( x)] F = ∫ U ( x)dF ( x) = ∫ U ( F −1 ( β ))dβ ≥
1
∞
0
−∞
≥ ∫ U (G −1 ( β ))dβ = ∫ U ( x)dG ( x) = E[U ( x)]G .
Таким образом, получили
E[U ( x)] F ≥ E[U ( x)]G .
доказывается аналогично. ■
37
Обратное утверждение
Из проведенного исследования можем сделать вывод, что использование
мер риска накладывает лишь одно ограничение на функцию полезности – это ее
неубывание (невозрастание), т. е. U ′( x) ≥ 0 ( U ′( x) ≤ 0 ).
2.5. Связь степеней рисковости и производных функции полезности
Лицо, принимающее решение, стремится получить наиболее выгодное
решение. Если ЛПР не склонен к риску, то для него наиболее выгодным будет
наименее рискованное решение, но если ЛПР предпочитает риск, то выберет
наиболее рискованную альтернативу. Чаще всего на практике ЛПР стремится
минимизировать риск, для достижения этой цели ЛПР старается принять решение
в условиях малорисковости или рисковости меньшей степени. Понятие степени
рисковости введено в работе [4].
Пусть Y – некоторая случайная величина, описывающая состояние
неблагоприятной внешней среды, ее функция распределения имеет вид
Φ( y) = P(Y ≤ y) . Вероятность принимает значения в интервале от 0 до 1, разобьем
данный
промежуток
на
n
интервалов
следующим
образом:
[k1 , k 2 ] ∪ (k 2 , k 3 ] ∪ ... ∪ (k n−1 , k n ] , где k1 = 0, k n = 1, а ki – значения вероятности,
заданные экспертом.
Определение 2.5.1. Будем полагать, что имеет место рисковость степени
ноль, если p( y) ∈ [k1 , k 2 ] .
Определение 2.5.2. Будем полагать, что имеет место рисковость степени
(n − 2) , если p( y ) ∈ [k n−1 , k n ] .
Как мы знаем, функция полезности характеризует отношение к риску ЛПР.
Например, положительная первая производная функции полезности описывает
случай ненасыщаемости потребностей, то есть для ЛПР с такой функцией
полезности больше лучше, чем меньше. А отрицательная вторая производная
функции полезности ЛПР показывает, что он строго не расположен к риску. Но
кроме склонности или несклонности индивидуума к риску, функция также
38
позволяет определить порядок неприятия риска, впервые данное понятие было
введено в работе [8].
Сначала дадим определение неприятия риска порядка n. Для любой
функции U (x) : [ a , b ] → R 1 предположим, что U ∈ Ñ ∞ , тогда U ( k ) ( x ) – производная
функции U (x) степени k, k = 2,3,...
Определение 2.5.3. Говорят, что ЛПР с функцией полезности U (x) имеет
неприятие риска порядка n, если
( −1) n −1U ( n ) ( x ) > 0 для всех x∈[a, b] , n ≥ 1 .
(2.5.1)
Теперь определим понятие порядка риска, предложенное в работе [8].
Пусть F (x) и G(x) – функции распределения случайных величин X 1 и X 2
соответственно, x∈[a, b] . Функция распределения принимает значения от 0 до 1,
т. е. F (a) = G(a) = 0 и F (b) = G(b) = 1. Введем следующие обозначения: F 1 ( x ) = F ( x )
x
, F ( x) = ∫ F k −1 (t ) dt , k = 2,3,... Эти обозначения применимы и к G(x) .
k
a
Определение 2.5.4. G(x) имеет большую величину риска порядка n, чем
F (x) , если
G k (b ) = F k (b ) , для k = 1,..., n − 1 ,
(2.5.2)
G n ( x ) ≥ F n ( x ) , для всех x∈[a, b] ,
(2.5.3)
со строгим неравенством > в (2.5.3) для некоторого x ∈ (a, b) .
В терминах ожидаемой полезности соотношение функций распределений,
приведенное выше, можно записать следующим образом:
b
b
n
a
a
k =1
1
1
k −1 ( k −1)
k
k
∫ U ( x)dF ( x) − ∫ U ( x)dG ( x) = ∑ (−1) U (b)[ F (b) − G (b)] +
b
+ ∫ (−1) n U ( n ) ( x)[ F n (b) − G n (b)]dx .
a
Рассмотрим
лемму,
которая
отражает
определений [8].
39
взаимосвязь
данных
выше
Лемма 2.5.1. G(x) имеет большую степень риска порядка n, чем F (x) , тогда
и только тогда, когда каждое лицо принимающее решение, имеющее неприятие
риска порядка n, предпочитает величину с функцией распределения F (x) величине
с функцией распределения G(x) .
В
статье
было
[9]
описано
соотношение
между
стохастическим
доминированием, порядками риска и неприятием риска.
Определение
2.5.5.
Любое
ЛПР,
имеющее
функцию
полезности,
удовлетворяющую условию
( −1) k −1U ( k ) ( x ) > 0 для всех x∈[a, b] и k = 1...n
предпочтет величину с функцией распределения F (x) величине с функцией
распределения G(x) тогда и только тогда, когда
G k (b ) ≥ F k (b ) , для k = 1,...n −1 ,
G n ( x ) ≥ F n ( x ) , для всех x∈[a, b] ,
со строгим неравенством > в (2.5.3) для некоторого x ∈ (a, b) .
Выше были изложены сведения необходимые для определения связи
степеней рисковости и производных функции полезности.
Теорема 2.5.1. Если ЛПР, характеризующееся неприятием риска порядка m,
предпочитает величину с функцией распределения F (x) величине с функцией
распределения G(x) в условиях рисковости степени (n − 1) , то ЛПР предпочитает
F (x) альтернативе G(x) и в условиях рисковости степени n.
Доказательство.
Пусть
имеет
место
ситуация,
когда
ЛПР,
характеризующееся неприятием риска порядка m, предпочитает F (x) , а не G(x) ,
т. е. для функции полезности ЛПР U (x) : [ a , b ] → R 1 предположим, что U ∈ Ñ ∞ ,
тогда U ( k ) ( x )
– производная функции U (x) степени k, k = 2,3,... выполняется
условие ( −1) m −1U ( m ) ( x ) > 0 для всех x∈[a, b] , m ≥ 1 .
Тогда по условию леммы 2.5.1 выполняется следующее:
большую величину риска порядка m, чем F (x) , если
G k (b ) = F k (b ) , для k = 1,..., m − 1 ,
40
G(x) имеет
G m ( x ) ≥ F m ( x ) , для всех x∈[a, b] ,
со строгим неравенством > в (2.5.3) для некоторого x ∈ (a, b) .
Пусть p1 ( y) – вероятность возникновения состояния y внешней среды,
которое соответствует условию рисковости степени (n − 1) , а p2 ( y ) соответствует
условию рисковости степени n.
С учетом независимости случайных событий G(x) имеет большую величину
риска порядка m, чем F (x) в условиях рисковости степени (n − 1) , обозначим это
как G m ( x) ≥ ( n−1) F m ( x) , если
G k (b ) = F k (b ) , для k = 1,..., m − 1 ,
G m ( x ) ⋅ p1 ( y ) ≥ F m ( x ) ⋅ p1 ( y ) , для всех x∈[a, b] .
Разделим обе части неравенства на p1 ( y) и умножим на p2 ( y ) . Так как
p1 ( y ) > 0 и p2 ( y) > 0 , то знак неравенства не изменится, получим
G m ( x ) ⋅ p 2 ( y ) ≥ F m ( x ) ⋅ p 2 ( y ) , для всех x∈[a, b] ,
следовательно G m ( x) ≥ (n) F m ( x) , что и требовалось доказать.
Теорема 2.5.2. Если ЛПР, характеризующееся неприятием риска порядка n,
предпочитает величину с функцией распределения F (x) величине с функцией
распределения G(x) , то ЛПР предпочитает F (x) альтернативе G(x) в случае
неприятия риска порядка (n + 1) .
Доказательство.
Пусть
имеет
место
ситуация,
когда
ЛПР,
характеризующееся неприятием риска порядка n , предпочитает F (x) , а не G(x) ,
т. е. для функции полезности ЛПР U (x) : [ a , b ] → R 1 предположим, что U ∈ Ñ ∞ ,
тогда U ( k ) ( x )
– производная функции U (x) степени k, k = 2,3,... выполняется
условие
( −1) n −1U ( n ) ( x ) > 0 для всех x∈[a, b] , n ≥ 1 .
Тогда по условию леммы 2.5.1 выполняется следующее:
(2.5.4)
G(x) имеет
большую величину риска порядка n , чем F (x) , если
G k (b ) = F k (b ) , для k = 1,..., n − 1 ,
41
(2.5.5)
G n ( x ) ≥ F n ( x ) , для всех x∈[a, b] ,
(2.5.6)
со строгим неравенством > в (2.5.3) для некоторого x ∈ (a, b) .
Знаем, что при выполнении условий (2.5.5, 2.5.6, 2.5.7), по определению
2.5.5 F (x) доминирует G(x) в условиях стохастического доминирования порядка
n, т. е. F ( x) ≥ n G ( x) .
Приведем здесь вспомогательную теорему, она доказана в работе [4].
Теорема 2.5.3. Если имеет место стохастическое доминирование порядка
(k − 1) : X 1 ≥ ( k −1) X 2 , то имеет место и стохастическое доминирование порядка k:
X 1 ≥ (k ) X 2 .
Тогда по теореме 2.5.3, так как выполняется F ( x) ≥ n G ( x) , то выполняется и
F ( x ) ≥ ( n+1) G ( x ) . Таким образом, для стохастического доминирования порядка
(n + 1) следует выполнение неравенства ( −1) n U ( n +1) ( x ) > 0 . Что и требовалось
доказать.
Теорема 2.5.4. Пусть α , β такие, что 1 ≤ α ≤ k , 0 ≤ β ≤ n − 3 , где k , n ∈ N . Если
найдутся такие i = α и j = β , что ЛПР, характеризующееся неприятием риска
порядка i, предпочитает величину с функцией распределения F (x) величине с
функцией распределения G(x) в условиях рисковости степени j, то имеет место
ситуация, когда ЛПР
с порядком неприятия риска (i + 1) предпочитает F (x)
альтернативе G(x) и в условиях рисковости степени ( j + 1) .
Введем обозначение предпочтения (i)RA( j) для ситуации, когда ЛПР,
характеризующееся неприятием риска порядка i, предпочитает величину с
функцией распределения F (x) величине с функцией распределения G(x) в
условиях рисковости степени j, т. е. F ( x) (i ) RA( j ) G ( x) .
(RA – risk aversion – неприятие риска).
Доказательство данной теоремы следует из доказательства теорем 2.5.1 и
2.5.2.
Данную взаимосвязь можно представить в виде схемы:
42
(1) RA(n − 2) ⇒ (2) RA(n − 2) ⇒ K ⇒ (k ) RA(n − 2),
↓
↓
↓
(1) RA( n − 3) ⇒ (2) RA(n − 3) ⇒ K ⇒ (k ) RA(n − 3),
↓
↓
↓
M
M
M
↓
(1) RA(0)
↓
(2) RA(0)
↓
(k ) RA(0).
⇒
⇒ K ⇒
Рис. 2.2. Схема, отражающая связь степеней рисковости и порядка неприятия риска
(порядка производной функции полезности)
43
Глава 3. Концепция стохастического доминирования
3.1. Принцип стохастического доминирования
Стохастическое доминирование является одной из концепций оценки в
условиях неопределенности. Стохастическое доминирование примечательно тем,
что для определения выгодной альтернативы достаточно знать соответствующие
функции распределения негарантированных результатов и отношение к риску
ЛПР. В обращении к явной функции полезности нет необходимости, хотя на нее
накладываются некоторые ограничения [5, 35].
Пусть задано вероятностное пространство (R, В, Р), где R – множество
вещественных чисел, В – борелевская сигма-алгебра, Р – вероятностная мера на В.
Пусть Р – множество вероятностных распределений на измеримом пространстве
(R, B). Тогда будем обозначать F ( x ) = Ρ (( −∞ , x ]) , функцию распределения,
соответствующую
X1 ∈ Ρ ,
G ( x ) = Ρ (( −∞ , x ]) ,
а
функцию
распределения,
соответствующую X 2 ∈ Ρ .
Определение 3.1.1. Пусть
X 1 , X 2 ∈ Ρ . Имеет место стохастическое
доминирование первого порядка: X 1 ≤ 1 X 2 , если F ( x) ≥ G ( x), x ∈ R .
Определение 3.1.2. Имеет место стохастическое доминирование k-го
порядка, X 1 ≤
k
k
k
X 2 , если F ( x ) ≥ G ( x ), x ∈ R , где
k
F ( x) =
x
∫F
k −1
(t ) dt ,
k = 2,3...
−∞
Стохастическое
доминирование
можно
характеризовать
различными
способами; используем описание в терминах ожидаемой полезности [6, 7].
1. Стохастическое доминирование первого порядка
Стохастическое доминирование первого порядка характеризует ситуацию,
когда ЛПР предпочитает большее количество полезности меньшему количеству,
т.е. функция полезности должна быть возрастающей.
44
Определение 3.1.3. Пусть U (x) – строго возрастающая монотонная функция
полезности (U ' ( x) > 0) , E [U ( x)]F – ожидаемая полезность, соответствующая
функции распределения F (x) . Пусть a ≤ x ≤ b , тогда говорят, что X 2 доминирует
X 1 в смысле стохастического доминирования первого порядка, если
F ( x) ≥ G ( x) ⇔ E[U ( x)] F ≤ E[U ( x )]G .
Стохастическое доминирование первого порядка принято обозначать как
≥ FSD или ≥1 . Графически условие проиллюстрировано на рис. 3.1.
F(x), G(x)
1
G(x)
F(x)
a
x
b
Рис. 3.1. Стохастическое доминирование первого порядка
2. Стохастическое доминирование второго порядка
Стохастическое
доминирование
второго
порядка
накладывает
дополнительное условие на функцию полезности: предполагается, что ЛПР
является строго не расположенным к риску, т. е. должны выполняться условия
U ' ( x) > 0 и U " ( x ) < 0 .
Определение 3.1.4. Пусть U (x) – строго возрастающая монотонная функция
полезности
(U ' ( x) > 0)
и
U " ( x) < 0 ,
E [U ( x)]F
соответствующая функции распределения
–
ожидаемая
полезность,
F (x) . Тогда говорят, что
X2
доминирует X 1 в смысле стохастического доминирования второго порядка, если
выполняется
x
x
−∞
−∞
∫ F (t ) dt ≥ ∫ G (t ) dt ⇔ E[U ( x)] F ≤ E[U ( x)]G .
45
Стохастическое доминирование второго порядка принято обозначать как
≥ SSD или ≥ 2 . Графически условие проиллюстрировано на рис. 3.2.
F(x), G(x)
1
G(x)
F(x)
x
a
b
Рис.3. 2. Стохастическое доминирование второго порядка
3. Стохастическое доминирование третьего порядка
Стохастическое доминирование третьего порядка накладывает еще одно
условие на функцию полезности: ЛПР обладает убывающей абсолютной
несклонностью к риску, это значит, что должно выполняться условие U ' ' ' ( x) > 0 .
Соответственно, для выполнения условий стохастического доминирования
третьего порядка для функции полезности должны выполняться условия U ' ( x) > 0,
U " ( x ) < 0 и U ' ' ' ( x) > 0 .
Определение 3.1.5. Пусть U (x) – строго возрастающая монотонная функция
полезности (U ' ( x) > 0) и U " ( x) < 0 , U ' ' ' ( x) > 0 , E [U ( x)]F – ожидаемая полезность,
соответствующая функции распределения
F (x) . Тогда говорят, что
X2
доминирует X 1 в смысле стохастического доминирования третьего порядка, если
выполняется
x ϑ
x ϑ
−∞ −∞
−∞ −∞
∫ ∫ F (v, t ) dvdt ≥ ∫ ∫ G (v, t ) dvdt ⇔ E[U ( x)] F
≤ E[U ( x)]G .
Стохастическое доминирование третьего порядка принято обозначать как
≥ TSD или ≥ 3 .
Можно развить концепцию стохастического доминирования до более
высоких порядков, при этом в условиях, накладываемых на функцию полезности,
присутствует следующая закономерность:
46
U
(k )
d k U ( x) > 0, если k нечетное
( x) =
=
.
dx k
< 0, если k четное
Но интерпретация данных условий является неоднозначной, а иногда и
невозможной [36].
3.2. Стохастическое доминирование и меры риска
К сожалению, зачастую различные меры риска дают противоречивые
результаты для одного и того же набора альтернатив. В свою очередь, принцип
максимизации
ожидаемой
полезности
и
критерий
стохастического
доминирования для функции полезности позволяют определить предпочтения
ЛПР и дают соответствующие правила принятия решений.
В данном разделе мы используем критерий стохастического доминирования
совместно с ранее выбранными мерами риска. Изучим поведение мер риска в
условиях
первого
и
второго
порядка
стохастического
доминирования
соответственно и узнаем, как меры риска способны дать согласованные
результаты.
Как мы знаем, различные меры риска дают отличные друг от друга
результаты. Аналогично, от выбора функции полезности для описания
предпочтений ЛПР зависит и результат принятия решения. Но если мы сможем
оценить варианты с помощью критерия стохастического доминирования, то для
широкого класса функций полезности предпочтения будут одни и те же.
Проследим взаимосвязь общей меры риска и ожидаемой полезности при условии
стохастического доминирования.
Рассмотрим две случайные величины X и Y с функциями распределения
F (x) и G ( y) соответственно. Будем обозначать X >1 Y , если X стохастически
доминирует Y в смысле первого порядка доминирования, и X > 2 Y , если
X
стохастически доминирует Y в смысле второго порядка доминирования. В общем
случае, обозначим X > SD Y ситуацию, когда X стохастически доминирует Y.
Дадим определение стохастического доминирования в терминах мер риска.
47
Определение 3.2.1. Говорят, что X стохастически доминирует Y – X > SD Y ,
если ρ ( X ) ≤ ρ (Y ) , где ρ – мера риска.
1. Меры риска и стохастическое доминирование первого порядка
Теорема 3.2.1.
1. Независимо от распределения X и Y , следующие выражения эквивалентны:
a. X >1 Y ;
b. VaRα ( X ) ≤ VaRα (Y ) ;
c. LPM n,t ( X ) ≤ LPM n,t (Y ) , где n = 0 .
2. Если X >1 Y , тогда:
d. ESα ( X ) ≤ ESα (Y ) ;
e. LPM n,t ( X ) ≤ LPM n,t (Y ) , где n = 1 .
Но обратное неверно.
Доказательство:
2 ⟺ 3: Если X и Y – случайные величины, для которых выполняется X >1 Y ,
следовательно, выполняется соотношение F ( x ) ≤ G ( y ) . Тогда для любого l ∈ R
справедливо следующее p ( X ≥ l ) ≤ p (Y ≥ l ) .
Из неравенства следует: если l принадлежит множеству {l ∈ R : p (Y ≤ l ) ≥ α } ,
тогда l принадлежит множеству {l ∈ R : p ( X ≤ l ) ≥ α} . Тогда выполняется
inf{l ∈ R : p ( X ≤ l ) ≥ α } ≤ inf{l ∈ R : p (Y ≤ l ) ≥ α } .
А так как VaRα ( X ) = inf{l ∈ R : p ( X ≤ l ) ≥ α} , то VaRα ( X ) ≤ VaRα (Y ) . Обратное
утверждение доказывается аналогично.
2 ⟺ 4:. При & = 0 мера принимает вид:
t
LPM 0,t ( X ) = ∫ (t − x) f ( x)dx =
0
−∞
t
∫ f ( x)dx = p( X ≤ t ) ,
−∞
48
аналогично она описывается и для величины Y. А так как X >1 Y , то F ( x ) ≤ G ( y ) ,
p ( X ≤ t ) ≤ p (Y ≤ t ) ,
следовательно
то
LPM n,t ( X ) ≤ LPM n,t (Y ) .
Обратное
доказывается аналогично.
d: Докажем, что если выполняется условие
X >1 Y , то выполняется и
ESα ( X ) ≤ ESα (Y ) .
Если
X >1 Y , то F ( x ) ≤ G ( y ) , следовательно, справедливо соотношение
{x : F ( x ) ≥ β } ⊆ { y : G ( y ) ≥ β } , тогда F −1 ( β ) ≥ G −1 ( β ) , где β ∈ (0 ;1] .
Напомним, что
ES α ( X ) =
1
α
α
∫ VaR β ( X , β ) dβ = −
0
1
α
F
α∫
−1
( β ) dβ .
Рассмотрим
0
разность мер для случайных величин X и Y ES α ( X ) − ES α (Y ) и определим знак
выражения.
1 α −1
ES α ( X ) − ES α (Y ) = − ∫ F ( β )dβ − − ∫ G ( β )dβ ) =
α0
α 0
1
−
1
α
(F
α∫
−1
α
−1
( β ) − G −1 ( β ))dβ
0
Так как F −1 ( β ) ≥ G −1 ( β ) , то получим, что F −1 ( β ) − G −1 ( β ) ≥ 0 , тогда
−
1
α
(F
α∫
−1
( β ) − G −1 ( β ))dβ ≤ 0 .
0
Отсюда следует, что ESα ( X ) ≤ ESα (Y ) .
e: Докажем, что если выполняется условие X >1 Y , то LPM 1,t ( X ) ≤ LPM 1,t (Y ) , где
n =1.
Если
X >1 Y , то
F ( x ) ≤ G ( y ) , а значит, выполняется неравенство
f ( x) ≤ g ( x) .
Запишем представление меры для непрерывного распределения:
t
LPM 1,t ( X ) =
∫ (t − x)
1
f ( x)dx .
−∞
Тогда LPM 1,t ( X ) − LPM 1,t (Y ) можно записать так
49
t
t
∫ (t − x) f ( x)dx − ∫ (t − x) g ( x)dx
−∞
=
−∞
t
= ∫ (t − x )( f ( x ) − g ( x ))dx
−∞
Тогда, так как выполняется f ( x) ≤ g ( x) , то верно и f ( x) − g ( x) ≤ 0 . Так как
при ∀t справедливо неравенство (t − x ) ≥ 0 . Отсюда следует, что выполняется
LPM 1,t ( X ) ≤ LPM 1,t (Y ) ■.
2. Меры риска и стохастическое доминирование второго порядка
Теорема 3.2.2. Независимо от распределения X и Y , следующие выражения
эквивалентны:
a. X > 2 Y ;
b. ESα ( X ) ≤ ESα (Y ) ;
c. LPM n,t ( X ) ≤ LPM n,t (Y ) , где n = 1 .
Если X > 2 Y , тогда выполняется LPM n,t ( X ) ≤ LPM n,t (Y ) , где n = 2 , но
обратное неверно.
Доказательство:
2 ⟺ 3: следует из результата, полученного Г. Фёльмером и А. Шидом в книге
[37]. В этой книге есть теорема 2.58,
которая утверждает, что X > 2 Y
эквивалентно неравенству:
α
∫F
−1
α
( β )dβ ≥ ∫ G −1 ( β ))dβ для всех α ∈ (0,1] .
0
0
А из данного неравенства следует:
−
1
α
F
α∫
−1
( β ) dβ ≤ −
0
Зная, что ES α ( X ) = −
1
α
G
α∫
−1
( β ))dβ для всех α ∈ (0,1] .
0
1
α
F
α∫
−1
( β ) dβ
, получаем ESα ( X ) ≤ ESα (Y ) .
0
50
3 ⟺ 4: по определению, так как X > 2 Y , то есть
t
t
∫ F ( x) dx ≥ ∫ G ( x) dx , а
−∞
−∞
t
LPM 1,t ( X ) = ∫ (t − x)1 f ( x) dx = (t − x ) F ( x) −∞ + ∫ F ( x )dx = ∫ F ( x) dx при условии,
−∞
−∞
−∞
t
t
t
t
xF ( x) = 0 , что следует из ∫ xf ( x)dx > −∞ .
что xlim
→ −∞
−∞
Докажем теперь: X > 2 Y ⇒ LPM2,t ( X ) ≤ LPM2,t (Y ) .
t
t
−∞
−∞
Так как X > 2 Y , то есть выполняется ∫ F ( z ) dz ≥ ∫ G ( z ) dz .
В работе Портера [38] доказано (данное доказательство можно провести с
помощью теоремы Фубини), что
t
t
1
F
(
z
)
dz
=
(t − z ) 2 f ( z )dz .
∫−∞
∫
2 −∞
(3.2.1)
Мы знаем, что
t
LPM 2,t ( X ) = ∫ (t − z ) 2 f ( z )dz ,
−∞
тогда из равенства 3.2.1 следует, что LPM 2,t ( X ) ≤ LPM 2,t (Y ) ■.
3.3. Стохастическое доминирование и степени рисковости
Во второй главе было дано определение 2.5.1 степени рисковости,
введенное в работе [4]. Так же в ней была доказана теорема, связывающая
понятие степени рисковости и порядок стохастического доминирования.
Теорема 3.3.1. Пусть α , β такие, что 1 ≤ α ≤ k − 1, 0 ≤ β ≤ n − 3 , где k , n ∈ N .
Если
найдутся
такие
i =α
и
j=β,
что имеет место стохастическое
доминирование порядка i в условиях рисковости степени j: A ≥ i , j B (или иначе
(i)SDR ( j) ),
то имеет место стохастическое доминирование порядка (i + 1) в
условиях рисковости степени ( j + 1) : A ≥ i +1, j +1 B (т. е. (i + 1)SDR( j + 1) ).
51
(1) SDR ( n − 2) ⇒ ( 2) SDR ( n − 2) ⇒ K ⇒ ( k ) SDR ( n − 2),
↓
↓
↓
(1) SDR ( n − 3) ⇒ ( 2) SDR ( n − 3) ⇒ K ⇒ ( k ) SDR ( n − 3),
↓
↓
↓
M
M
M
↓
(1) SDR (0)
↓
( 2) SDR (0)
↓
(k ) SDR (0).
⇒
⇒ K ⇒
Рис. 3.1. Схема, отражающая связь степеней рисковости и порядка стохастического
доминирования
52
Глава 4. Сравнение методов
Методы оценивания рисков, а так же ранжирования оценок рисков,
рассмотренные в работе, чаще всего применяются при решении финансовых
задач, но это не единственная область их применения. Так же их можно
использовать для оценки количества людей с различными заболеваниями,
основываясь на статистических данных Министерства Здравоохранения. В
частности, проведем оценку риска заражения людей ВИЧ и туберкулезом.
Постановка
задачи.
С
помощью
изученных
методов:
принцип
стохастического доминирования, метод ожидаемой полезности, оценка с
помощью мер риска, сравнить случайные величины в условиях различных
степеней рисковости, реализации которых приведены в Приложении, табл. 3, это
статистические данные по количеству людей с диагнозами ВИЧ и туберкулез [39].
Решение. Обозначим X 1 – случайную величину, которая соответствует
данным о количестве людей с диагнозом ВИЧ, а X 2 – случайную величину,
отражающую количество людей с диагнозом туберкулез. И запишем функции
распределения данных величин
F (x) и G (x) соответственно, их графики
представлены на рис. 4.1.
F(x) и G(x)
1,00
G(x)
ВИЧ
0,80
F(x)
0,60
Туберкулез
0,40
0,20
0,00
0
100 000 200 000 300 000 400 000 500 000 600 000 700 000
x
Рис. 4.1. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез).
53
0 при x < 887
0,05 при 887 ≤ x < 1090
0,11 при 1090 ≤ x < 2063
0,16 при 2063 ≤ x < 6918
0,21 при 6918 ≤ x < 10889
0,26 при 10889 ≤ x < 30647
0,32 при 30647 ≤ x < 89808
0,37 при 89808 ≤ x < 177579
0,42 при 177579 ≤ x < 227502
0,47 при 227502 ≤ x < 263898
F ( x) =
0,53 при 263898 ≤ x < 296045
0,58 при 296045 ≤ x < 300252
0,63 при 300252 ≤ x < 373718
0,68 при 373718 ≤ x < 403100
0,74 при 403100 ≤ x < 471676
0,79 при 471676 ≤ x < 530185
0,84 при 530185 ≤ x < 589581
0,89 при 589581 ≤ x < 650100
0,95 при 650100 ≤ x < 703781
1 при x ≥ 703781
0 при x < 94214
0,05 при 94214 ≤ x < 101570
0,11 при 101570 ≤ x < 115792
0,16 при 115792 ≤ x < 142411
0,21 при 142411 ≤ x < 142871
0,26 при 142871 ≤ x < 143293
0,32 при 143293 ≤ x < 148762
0,37 при 148762 ≤ x < 150879
0,42 при 150879 ≤ x < 150904
0,47 при 150904 ≤ x < 154344
G ( x) =
0,53 при 154344 ≤ x < 154964
0,58 при 154964 ≤ x < 157110
0,63 при 157110 ≤ x < 157649
0,68 при 157649 ≤ x < 157859
0,74 при 157859 ≤ x < 159110
0,79 при 159110 ≤ x < 160130
0,84 при 160130 ≤ x < 163618
0,89 при 163618 ≤ x < 165401
0,95 при 165401 ≤ x < 171352
1 при x ≥ 171352
Зададим значения вероятности, которые разбивают промежуток от 0 до 1 на
некоторое количество интервалов и воспользуемся определением 2.4.2. Разделим
промежуток на пять интервалов, которые будут характеризовать риск следующим
образом:
•
p( y ) ∈ [0; 0,2] – рисковость степени ноль (низкая степень риска);
•
p( y ) ∈ (0,2; 0,4] – рисковость первой степени (степень риска ниже средней)
•
p( y ) ∈ (0,4; 0,6] – рисковость второй степени (средняя степень риска)
•
p( y ) ∈ (0,6; 0,8] – рисковость третьей степени (степень риска выше средней)
•
p( y ) ∈ (0,8;1] – рисковость четвертой степени (высокая степень риска)
Тогда сравним случайные величины на каждом из промежутков с помощью
принципа стохастического доминирования, метода ожидаемой полезности и
аппарата мер риска, для этого воспользуемся теоремой 2.4.3 и теоремой 3.3.1.
54
1.
Принцип стохастического доминирования
Проверим условие стохастического доминирования первого и второго
порядка. По определениям 3.1.1 и 3.1.2, справедливо соотношение X 1 ≤ 1 X 2 если
функции удовлетворяют неравенству F ( x) ≥ G ( x ) для любых x ∈ R , а X 1 ≤ 2 X 2
x
справедливо, если
x
∫ F (t ) dt ≥ ∫ G (t ) dt
−∞
для любых x ∈ R .
−∞
Рассмотрим промежуток [0; 0,2] , проверим, как себя ведут функции
распределения на данном промежутке (рис. 4.2).
F(x) и G(x)
0,18
0,16
F(x)
0,14
0,12
0,10
ВИЧ
G(x)
0,08
Туберкулез
0,06
0,04
0,02
0,00
x
0
20 000
40 000
60 000
80 000
100 000
120 000
Рис. 4.2. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке [0; 0,2] (низкая степень рисковости).
По графику видно, что F ( x) ≤ G ( x ) на всем промежутке p( y ) ∈ [0; 0,2].
Таким образом, условие для первого порядка стохастического доминирования
выполнено, т. е. X 1 ≥ 1 X 2 на данном промежутке.
Теперь проверим выполнение условия стохастического доминирования
второго порядка. В статье [40] приведена формула для вычисления интеграла
дискретной функции распределения:
F 2 ( xi ) = F ( xi )( xi − xi −1 ) + F ( xi −1 ) ,
55
(4.1.1)
здесь F 2 ( xi ) интеграл функции распределения F (x ) для наблюдения xi .
используем данную формулу для поиска интегралов функций распределения, и
затем сравним их. Запишем результаты расчетов в виде таблицы
F 2 (887) = 0,05 ⋅ (887 − 0) + 0 = 46,68
F 2 (1090) = 0,11 ⋅ (1090 − 887 ) + F 1 (887 ) = 68,05
F 2 (2063) = 0,16 ⋅ ( 2063 − 1090) + F 1 (1090) = 306,95
G 2 (94214 ) = 0,05 ⋅ (94214 − 0) + 0 = 4958,63
G 2 (101570 ) = 0,11 ⋅ (101570 − 94214 ) + G 1 (94214 ) = 5732,91
G 2 (115792 ) = 0,16 ⋅ (115792 − 101570 ) + G 1 (101570 ) = 7978,48
Для удобства сравнения запишем вычисления в табличном виде:
Соответствующее
значение вероятности
0,05
0,11
0,16
Исходя
из
полученных
F 2 ( x)
G 2 ( x)
F 2 ( x) − G 2 ( x)
46,68
68,05
306,95
4958,63
5732,91
7978,48
-4 911,95
-5 664,86
-7 671,53
результатов,
можем
сделать
вывод,
что
F 2 ( x ) ≤ G 2 ( x) , значит, в условиях рисковости степени ноль X 1 доминирует X 2 в
0
смысле стохастического доминирования второго порядка, т. е. X 1 ≥ 2 X 2 .
Рассмотрим промежуток (0,2; 0,4] , характеризующий первую степень
рисковости, на рис. 4.3 представлен график функций распределения на этом
промежутке. По графику можем увидеть, что и на рассматриваемом промежутке
F ( x) ≤ G ( x ) , т. е. X 1 ≥ 1 X 2 . То есть в заданных условиях при уровне вероятности
p( y ) ∈ (0,2; 0,4] ВИЧ менее опасен, чем туберкулез.
56
F(x) и G(x)
0,40
0,38
F(x)
0,36
0,34
0,32
G(x)
0,30
ВИЧ
0,28
Туберкулез
0,26
0,24
0,22
0,20
0
50 000
100 000
150 000
200 000 x
Рис. 4.3. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке (0,2; 0,4] (степень рисковости ниже средней).
Приведем значения интегралов, соответствующих данному промежутку для
того, чтобы определить, выполняется ли условие стохастического доминирования
второго порядка. Интегралы были рассчитаны по формуле (4.1.1).
Соответствующее
значение вероятности
0,21
0,26
0,32
0,37
F 2 ( x)
G 2 ( x)
F 2 ( x) − G 2 ( x)
1 215,37
2 260,37
8 499,74
30295,89
8530,03
11282,59
13626,33
17483,16
-7 314,66
-9 022,22
-5 126,60
12 812,74
Можно увидеть, что отрицательная разность F 2 ( x ) − G 2 ( x) в сумме
больше, чем положительная, т. е. F 2 ( x) ≤ G 2 ( x) . К такому же выводу можно
прийти и исходя из графика функций распределения на рис. 4.4. Можно заметить,
что площадь ограниченная функциями распределения обозначенная знаком плюс,
меньше, чем площадь, обозначенная знаком минус. Таким образом,
доминирует X 2 в смысле стохастического доминирования второго порядка.
57
X1
Рис. 4.4. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке (0,2; 0,4] (степень рисковости ниже средней).
Для второй степени рисковости, на рис. 4.5 представлен график функций
распределения на этом промежутке.
F(x) и G(x)
0,60
F(x)
0,55
G(x)
ВИЧ
0,50
Туберкулез
0,45
0,40
0
100 000
200 000
300 000
400 000 x
Рис. 4.5. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке (0,4; 0,6] (средняя степень рисковости).
По графику можем увидеть, что на рассматриваемом промежутке
F ( x) ≥ G ( x ) , т. е.
X 1 ≤ 1 X 2 . В данном интервале вероятности туберкулез
является менее опасным, чем ВИЧ.
58
Проверим, выполняется ли условие стохастического доминирования
второго порядка, для этого вычислим интегралы функций распределения:
Соответствующее
значение вероятности
0,42
0,47
0,53
0,58
F 2 ( x)
G 2 ( x)
F 2 ( x) − G 2 ( x)
67 252,11
90 899,84
110 055,63
128 667,05
19 216,98
19 228,55
21 039,13
21 398,32
48 035,13
71 671,29
89 016,50
107 268,73
По полученным значениям можно сделать вывод, что F 2 ( x ) ≥ G 2 ( x) , то
есть X 2 доминирует X 1 в смысле стохастического доминирования второго
порядка.
На рис. 4.6 и 4.7 представлены графики функций на промежутках (0,6; 0,8]
и (0,8;1] , относящихся к третьей и четвертой степени рисковости соответственно.
F(x) и G(x)
0,80
0,78
0,76
0,74
0,72
0,70
0,68
0,66
0,64
0,62
0,60
F(x)
G(x)
ВИЧ
Туберкулез
0
100 000
200 000
300 000
400 000
500 000 x
Рис. 4.6. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке (0,6; 0,8] (степень рисковости выше средней).
Поведение функций на данных промежутках подобно, здесь величина X 2
доминирует X 1 в смысле стохастического доминирования первого порядка, это
можно увидеть на графиках.
59
F(x) и G(x)
1,00
F(x)
0,95
G(x)
ВИЧ
0,90
Туберкулез
0,85
0,80
0
200 000
400 000
600 000
800 000
x
Рис. 4.7. Функции распределения случайных величин X 1 (ВИЧ), X 2 (Туберкулез) на
промежутке (0,8; 1] (высокая степень рисковости).
Приведем ниже значения интегралов для каждого из интервалов.
Соответствующее
F 2 ( x)
G 2 ( x)
значение вероятности
Интервал (0,6; 0,8]
0,63
131 324,11 22753,23
0,68
181 590,32 23122,27
0,74
203 240,21 23277,21
0,79
257 379,16 24264,81
Интервал (0,8;1]
0,84
306 649,89 25123,41
0,89
359 793,68 28244,04
0,95
417 127,47 29933,34
1
470 808,47 35884,89
F 2 ( x) − G 2 ( x)
108 570,88
158 468,05
179 963,00
233 114,35
281 526,49
331 549,65
387 194,14
434 923,59
Исходя из расчетов, можем сделать вывод, что выполняется следующее
неравенство F 2 ( x ) ≥ G 2 ( x) , следовательно, величина X 2 доминирует X 1 в
смысле стохастического доминирования
вероятности,
попадающих
в
второго
рассмотренные
порядка. Для
промежутки
ВИЧ
значений
опаснее
туберкулеза.
Если же не использовать понятие степени рисковости, то результат для
принципа стохастического доминирования будет следующим.
60
Как можно увидеть по графику на рис.1 для функций F (x) и G (x) условие
принципа стохастического доминирования первого порядка не выполняется. Так
как функции не удовлетворяют неравенству F ( x) ≥ G ( x ) для любых x ∈ R ,
которое приведено в определении 3.1.1.
Проверим,
выполняется
ли
условие
принципа
стохастического
доминирования второго
рого порядка. Исходя из графика функций распределения на
рис. 8, можно заметить, что площадь ограниченная функциями распределения
обозначенная знаком плюс, значительно меньше, чем площадь, обозначенная
знаком минус, справа от точки пересечения функций распределения
расп
. Поэтому
можно сделать вывод, что при вычислении разности функций F ( x ) − G ( x ) для
всех x ∈ R , получим функцию, которая принимает отрицательные значения, а это
2
2
означает F ( x) ≥ G ( x) . Это говорит о том, что для ЛПР не склонного к риску
предпочтительнее будет альтернатива
G ( x ) , соответственно, по данному
принципу ВИЧ опаснее туберкулеза.
F(x) и G(x)
1,00
G(x)
0,80
‒
ВИЧ
F(x)
0,60
Туберкулез
0,40
+
0,20
0,00
x
0
100 000 200 000 300 000 400 000 500 000 600 000 700 000
Рис. 4.8. Функции распределения случайных величин X 1 (ВИЧ),
ВИЧ), X 2 (Туберкулез).
2.
Метод ожидаемой полезности
Выберем функцию ожидаемой полезности, она должна неограниченно
увеличиваться при возрастании аргумента, причем достаточно быстро, так как,
61
чем больше людей заболевает, тем хуже. Выбирать будем, соответственно
наименьшее из полученных значений функции полезности для величин с
функциями распределения F (x) и G (x) .
Тогда пусть функция полезности имеет следующий вид:
U ( x ) = x 0 ,9
Данная функция дифференцируема U ∈ Ñ ∞ , причем U ′( x ) = 0,9 ⋅ x −0,1 > 0,
U ′′( x ) = − 0 , 09 ⋅ x − 1,1 < 0 , U ′′′( x ) = 0,099 ⋅ x −2,1 > 0 , … ( −1) m −1U ( m ) ( x ) > 0 .
Ниже изображен график данной функции на рис. 9.
Рис. 4.9. График функции полезности в общем виде.
Тогда рассчитаем для наблюдений случайных величин X 1 и X 2 значения
функции полезности для каждой степени рисковости. Расчеты приведены в
табл. 4 Приложения.
После этого мы можем найти ожидаемую полезность
каждой альтернативы по формуле:
n
E[U ( x )] = ∑ U ( x i ) ⋅ p i ,
i =1
где n – количество наблюдений, pi – вероятность наблюдения xi .
Для общего случая, без разделения по степеням рисковости, получены
следующие результаты для статистик:
E [U F ( x)] = 74043,76 и E [U G ( x)] = 43991,5 ,
E [U F ( x)] ≥ E [U G ( x)] .
62
Для ЛПР в данном случае предпочтительна ожидаемая полезность с
минимальным значением, т. е. E [U G ( x)] . Соответственно, по методу ожидаемой
полезности туберкулез является менее опасным заболеванием, чем ВИЧ.
Теперь рассчитаем значения ожидаемой полезности для каждой степени
рисковости:
Степень
рисковости
Нулевая степень
E [U F ( x)]
E [U G ( x)]
Соотношение полезностей
114,58
5164,41
E [U F ( x)] ≤ E [U G ( x)]
Первая степень
2461,61
8578,86
E [U F ( x)] ≤ E [U G ( x)]
Вторая степень
14684,74
9748,66
E [U F ( x)] ≥ E [U G ( x)]
Третья степень
22487,17
10044,52
E [U F ( x)] ≥ E [U G ( x)]
Четвертая степень
34295,66
10455,05
E [U F ( x)] ≥ E [U G ( x)]
Таким образом, для нулевой и первой степеней рисковости ВИЧ является
менее
опасным,
чем
туберкулез.
Для
остальных
степеней
ситуация
противоположная. Подводя итог, можно сказать, что при малых значениях
вероятности возникновения болезни нужно опасаться туберкулеза, а при больших
значениях – ВИЧ.
3.
Оценивание альтернатив с помощью мер риска
По результатам, приведенным в главе 1, мы сделали вывод, что наилучшей
когерентной мерой для расчета риска эмпирических наблюдений является
ES a( 2 ) ( X ) , которая вычисляется по формуле (1.4.2), из некогерентных мер
рассмотрим меру LPM n,t ( X ) , которую будем рассчитывать по формуле (1.3.3).
Расчет мер проводим так же, как и для задачи, рассмотренной в главе 1.
Возьмем уровень вероятности α как среднее значение каждого промежутка,
характеризующего степени рисковости. После проведенных в Microsoft Excel
расчетов, получены следующие результаты:
63
Болезнь
ВИЧ
Туберкулез
Болезнь
ВИЧ
Туберкулез
ES a( 2 ) ( X )
α = 0,1
983
97698
α = 0,3
7691
114783
α = 0,5
71565
128512
α = 0,7
143065
136475
α = 0,9
224643
141814
α = 0,7
161921
14825
α = 0,9
328763
21228
LPM n,t ( X )
α = 0,1
0
0
α = 0,3
1687
4500
α = 0,5
78926
11286
По ним можно сделать следующий вывод. Для нулевой и первой степени
рисковости наиболее опасным является туберкулез. Для второй степени
рисковости
ситуация
складывается
неоднозначная,
так
как
меры
дают
противоположное упорядочение. Это объясняется тем, что на промежутке
(0,4; 0,6] функции распределения пересекаются. Для третьей и четвертой
степеней рисковости наиболее опасным является ВИЧ, это подтверждается и
графиком.
Результаты. Итак, нами была проведена оценка риска тремя методами,
сравним полученные результаты для каждой степени рисковости.
Степень рисковости ноль: все методы дают единый
результат – для
значения вероятности принадлежащего промежутку [0; 0,2] ВИЧ будет менее
опасен, чем туберкулез.
Первая степень рисковости: туберкулез опаснее, чем ВИЧ для уровня
вероятности из интервала
(0,2; 0,4] . Такой результат получен для всех
используемых методов. Причем, для стохастического доминирования второго
порядка получены разности интегралов с разными знаками, это говорит о том, что
при вероятности близкой к 0,4 функции пересекаются.
Вторая степень рисковости: для принципа стохастического доминирования
первого и второго порядка, а так же для метода ожидаемой полезности при
вероятности (0,4; 0,6] туберкулез менее опасен, чем ВИЧ. Но по результатам
применения аппарата мер риска однозначно определить, какая из болезней
64
опаснее, нельзя. Так как меры дали противоположные результаты, это произошло
потому, что при вероятности близкой к 0,4 функции пересекаются и не позволяют
ответить на вопрос при расчете мер только для середины отрезка (0,4; 0,6] . Чтобы
уточнить результат, нужно было бы рассчитать меры риска для других уровней
вероятности данного промежутка.
Третья и четвертая степень рисковости: для данных степеней результаты
одинаковы для всех методов. Во всех случаях, когда вероятность высока ВИЧ
опаснее туберкулеза.
Если бы мы проводили наше исследование без использования понятия
степеней рисковости, результаты получились бы совершенно другими. Так, по
принципу стохастического доминирования первого порядка нельзя сказать, какая
из болезней опаснее, так как условия этого принципа не выполнены для всех
x ∈ R . По принципу стохастического доминирования второго порядка можно
сделать вывод, что для ЛПР не склонного к риску ВИЧ опаснее туберкулеза. И
метод ожидаемой полезности дает такой же результат. А результат, получаемый с
помощью аппарата мер риска зависит от того, какой уровень вероятности выберет
ЛПР. В этом смысле, меры риска являются более гибким инструментом для
оценки риска и упорядочения альтернатив, даже без использования понятия
степеней рисковости. Так же можно сделать вывод, что применение степеней
рисковости помогает более качественно и подробно изучить соотношение
альтернатив, особенно для методов стохастического доминирования и ожидаемой
полезности, что очень важно в условиях реальной жизни. Результаты
исследования опубликованы в книге [28].
65
Глава 5. Оценка интенсивности предпочтения
5.1. Некоторые свойства интенсивности предпочтения
Мы знаем, что функция полезности является возрастающей (или
неубывающей), если с ростом количества ресурсов ЛПР растет их полезность. И
наоборот, функция будет убывающей (или невозрастающей), если при росте
количества ресурсов полезность будет уменьшаться.
Рассмотрим функцию
полезности u ( x) : Ω → R1 предположим, что u ∈ С ∞ , тогда u ( k ) ( x ) – производная
функции u ( x) степени k, k = 2,3,... Введем определение степени интенсивности
предпочтения.
Определение 5.1.1. Когда для отношения f и {xi } = Ω
принятия
решений
существует
функция
–
множества исходов
полезности
такая,
что
x i f x j ⇒ u ( x i ) > u ( x j ) , будем называть разность полезностей u ( x i ) − u ( x j ) > 0
интенсивностью предпочтения [17].
Обозначим u ( x i , x j ) = u ( x i ) − u ( x j ) – функция интенсивности предпочтения,
тогда условие будет выглядеть следующим образом:
xi f x j ⇒ u ( xi , x j ) > 0 для всех xi , x j ∈ Ω .
Бинарное отношение ≺ может быть описано с помощью функции
интенсивности предпочтения, когда оно ациклично.
Оценка интенсивности предпочтения играет решающую роль в ситуациях,
когда нет единственной наилучшей альтернативы. Впервые ситуацию (парадокс),
где отсутствует максимальная альтернатива, описал Кондорсе в конце XVIII века.
В качестве примера такой ситуации является случай, когда множество Ω состоит
из трех элементов , , 5 и имеется три таких строгих упорядочения {≺# }7#6
что
≺
≺ 5,
≺8 5 ≺ 8 ,
5 ≺7
≺7 .
66
на Ω,
В этом случае отношение ≺, представляющее собой:
где & ( , ) – разность между числом таких i, для которых
i, для которых
9 # , имеет вид:
≺
≺
≺5≺ .
≺
⟺ & ( , )< 0,
9 # , и числом таких
≺
Получается, что при очередном принятии решения в данных условиях
может быть выбрана любая из альтернатив, и, в этом смысле, устойчивости
решение этой задачи не имеет. Но оценить наиболее вероятное к принятию
решение можно с использованием интенсивности предпочтения. Причем,
существуют ситуации, когда оценка интенсивности предпочтения позволяет
принять оптимальное решение.
Функция
интенсивности
предпочтения
позволяет
упорядочивать
альтернативы по величине интенсивности. Она задает на множестве Ω×Ω слабое
упорядочение:
( , ) ≺∗ (5, ;) ⟺ ( , ) < (5, ;);
при этом условие о кососимметричности функции u означает, что
( , ) ≺∗ (5, ;) ⟺ (;, 5) ≺∗ ( , ).
Из определения 5.1.1. получим:
≺
⟺ ( , ) ≺∗ ( , ) .
В книге [17] предложен такой подход к оценке интенсивности, при котором на
отношение предпочтения и множество альтернатив накладываются определенные
условия, при которых функция u, обладающая свойством:
( , ) ≺∗ (5, ;) ⟺ ( ) − ( ) < (5) − (;)
единственна с точностью до положительного линейного преобразования.
Формулировки алгебраических аксиом, обладающих указанными свойствами,
приведены в [3].
Приведем в данном разделе некоторые алгебраические и порядковые свойства
интенсивности предпочтения, опубликованные в [17].
Определение 5.1.2. Пусть на Ω задано отношение предпочтения ≺, тогда
функция интенсивности предпочтения u, представляющая данное отношение
будем называть
67
аддитивной на транзитивных компонентах, если
≺
≺5∧
≺ 5 ⟹ ( , 5) = ( , ) + ( , 5);
аддитивной на отрицательно транзитивных компонентах, если
9
≺
9 5∧
≺
9 5 ⟹ ( , 5) = ( , ) + ( , 5);
≺
супераддитивной на транзитивных компонентах, если
≺
≺5∧
≺ 5 ⟹ (5, ) ≥ (5, ) + ( , );
супераддитивной на отрицательно транзитивных компонентах, если
9
≺
9 5∧
≺
9 5 ⟹ ( , 5) ≥ ( , ) + ( , 5).
≺
Если для данного отношения ≺ на Ω задана односторонняя функция
полезности >:
≺
⟹ >( ) < >( ),
то функция интенсивности предпочтения имеет вид:
где A( , ) = 0 при
≈
( , ) = ?>( ) − >( )@A( , ),
и A( , ) = 1 при
≈
9 , и является аддитивной на
транзитивных компонентах. Кроме того она аддитивна на отрицательно
транзитивных компонентах при выполнении условия
≺
⟺ >( ) < >( ).
Для отношения ≺ на Ω, задающего интервальное упорядочение
≺
⟺ >( ) + B( ) < >( ),
где > - функция полезности, а B - неотрицательная пороговая функция, функция
интенсивности принимает следующий вид
и
задает
>( ) + B( ) − >( ) при ≺ ,
( , ) = C0 при ≈ , G
>( ) − >( ) − B( ) при ≻
супераддитивное
отношения ≺, так как
≺
на
транзитивных
компонентах
представление
≺ 5 ⟹ ( , 5) = ( , ) + ( , 5) + A( ).
Благодаря свойству аддитивности функции интенсивности предпочтения,
ЛПР может восстановить функцию полезности. Приведем описание ниже.
68
Лемма 5.1.1. На Ω существует такая функция полезности > единственная с
точностью до прибавления константы, что для любых ,
( , ) = >( ) − >( ),
∈ Ω выполняется
если функция интенсивности предпочтения u аддитивна на отрицательно
транзитивных компонентах и представляет слабое упорядочение ≺ на Ω.
Лемма 5.1.2. На Ω существует функция полезности > такая, что
выполняется
т.е. для любых ,
∈Ω
где A( , ) = 0 при
⟹ >( ) − >( ) = ( , ),
≺
≈
( , ) = ?>( ) − >( )@A( , ),
и A( , ) = 1 при
≈
9 , если функция интенсивности
предпочтения u является аддитивной на транзитивных компонентах и задает
строгое упорядочение ≺ на Ω.
Доказательство этих лемм можно найти, например, в [17].
Определение
Пусть
5.1.3.
на
Ω
задано
отношение
≺,
тогда
представляющую его функцию интенсивности предпочтения u будем называть
монотонной на транзитивных компонентах, если
≺
≺5∧
≺ 5 ⟹ (5, ) ≥ max { (5, ), ( , )};
если же неравенство в правой части данного условия является строгим, то
функцию интенсивности предпочтения u назовем строго монотонной на
транзитивных компонентах.
≺
Определение 5.1.4. Пусть на Ω задано отношение ≺, которое имеет цикл
8
≺⋯≺
!
≺
, тогда если существует такой элемент 5 ∈ Ω, что
при всех i, или существует такой элемент ; ∈ Ω, что ; ≺
≺
8
≺⋯≺
!
≺
называется f-циклом.
#
#
≺5
при всех i, то цикл
Теорема 5.1.1. Пусть Ω есть не более чем счетное множество альтернатив,
на котором задано асимметричное отношение ≺. Тогда
1) ≺ не имеет f-циклов, если оно задается строго монотонной на транзитивных
компонентах функцией интенсивности предпочтения,
69
2) ≺ представимо с помощью аддитивной на транзитивных компонентах
функции интенсивности предпочтения, если ≺ не имеет f-циклов и является
слабо полным.
На множестве альтернатив Ω, являющемся выпуклым подмножеством в
некотором векторном пространстве E, интенсивность предпочтения обладает
рядом особых свойств, которые будут описаны ниже. [17].
Определение 5.1.4. Функцию интенсивности предпочтения u на Ω×Ω будем
называть биаффинной, если выполняются следующие равенства
( ,
при любом
(
+ (1 − )5) =
+ (1 − ) , 5) =
∈ L0,1M и при любых , , 5 ∈ Ω
( , ) + (1 − ) ( , 5),
( , 5) + (1 − ) ( , 5),
Так как функция u кососимметрична, то условия этого определения
эквивалентны.
Теорема 5.1.2 (о единственности). Пусть существует биаффинная функция
интенсивности предпочтения u, представляющая на выпуклом множестве
альтернатив Ω отношение ≺, тогда существует такая константа
> 0, что для
любой другой биаффинной функции интенсивности предпочтения > на Ω×Ω
выполняется
( , ) = >( , ), для всех ,
∈Ω
Так как в теории полезности основным предположением является
монотонность функции интенсивности предпочтения, то рассмотрим связь между
биаффинностью и монотонностью на транзитивных компонентах.
Теорема 5.1.3 (о монотонности). Пусть на Ω задано отношение ≺, которое
представляется биаффинной функцией интенсивности предпочтения u, которая
будет монотонной на транзитивных компонентах тогда и только тогда, когда
отношение ≺ – слабое упорядочение и на Ω существует такая аффинная функция
>, что для всех ,
∈ Ω выполняется
( , ) = >( ) − >( ).
Доказательства теорем можно найти в [4].
70
Определение 5.1.5. Функцию интенсивности предпочтения u будем
называть
выпукло-вогнутой, если для любых , , 5 ∈ Ω и любого
( ,
+ (1 − )5) ≥
строго выпукло-вогнутой, если при
( ,
+ (1 − )5) >
∈ L0,1M
( , ) + (1 − ) ( , 5);
∈ L0,1M и
≠5
( , ) + (1 − ) ( , 5);
квазивыпукло-квазивогнутой, если для любых , , 5 ∈ Ω и любого
( ,
+ (1 − )5) ≥ min { ( , ), ( , 5)}.
∈ L0,1M
Функция интенсивности предпочтения может быть использована для поиска
максимальных элементов отношения ≺, данное свойство функции будет описано
ниже.
Теорема 5.1.4. Пусть множество альтернатив Ω является выпуклым
компактом, на котором задано отношение ≺, обладающее представлением с
помощью функции интенсивности предпочтения u. Отношение ≺ имеет
максимальный элемент, если функция u является квазивыпукло-квазивогнутой и
полунепрерывна снизу по второй переменной.
Отношение
≾ = ≺∪≈ имеет
строго наибольший элемент, если функция u является строго выпукло-вогнутой.
5.2. Степени и меры интенсивности предпочтения
Концепция интенсивности предпочтения развивается и критикуется уже
долгое время. Критика данной концепции касается способов измерения
интенсивности. Различные способы оценки интенсивности приведены в работах
[3, 18, 19, 28, 41]. В настоящей работе мы предложили и исследовали новые
критерии оценки интенсивности предпочтения: степени и меры интенсивности.
•
Производные функции полезности и интенсивность предпочтения
Будем рассматривать три группы ЛПР, сформированные в зависимости от
отношения к риску: склонные к риску – рискофилы, несклонные – рискофобы и
безразличные – рисконейтралы.
71
В статье [42] приводится интерпретация знаков пяти первых производных
функции полезности. А так же обосновывается теория о том, что рискофилы и
рискофобы имеют одинаковое отношение к риску для нечетных порядков. А для
четных порядков их отношение противоположны: для рискофобов u ( 2n) ( x) < 0 , а
для рискофилов u ( 2 n ) ( x ) > 0 .
Рассмотрим
зависимость
предпочтений
от
производных
функции
полезности, то есть от отношения ЛПР к риску.
‒ Нечетные порядки производных функции полезности
Известно, что первая производная для любого ЛПР должна быть больше
нуля, причем, чем больше значение производной, тем больше полезность. Это
обусловлено свойствами функции полезности, аксиомами фон НейманаМоргенштерна о рациональном выборе ЛПР [31].
Аналогично, для нечетных порядков функции полезности действует
следующая интерпретация: чем больше значение производной, тем выше
полезность. Следовательно, интенсивность предпочтения для нечетных порядков
производной функции полезности, будет выражаться так:
xi f x j ⇒ u ( 2n −1) ( xi ) − u ( 2n −1) ( x j ) > 0 .
‒ Четные порядки производных функции полезности
Как известно, чем больше абсолютное значение второй производной, тем
больше несклонность / склонность к риску ЛПР [8]. Кроме того, знак второй
производной функции полезности определяет отношение ЛПР к риску: если
u" ( x) < 0 , то ЛПР является рискофобом, если u" ( x) > 0 , то ЛПР – рискофил.
Тогда интенсивность предпочтения для разных типов отношения к риску
будет выглядеть следующим образом:
‒ для рискофила: x i f x j ⇒ u" ( x i ) − u" ( x j ) > 0 .
‒ для рискофоба: x i f x j ⇒ u" ( x i ) − u" ( x j ) < 0 .
Можно показать, что приведенные выше рассуждения, распространяются и
на более высокие четные порядки производных.
Введем определение степени интенсивности предпочтения.
72
Определение 5.2.1. Имеет место интенсивность предпочтения степени 0
по модулю µ 0 , если
xi f x j ⇒ u ( xi ) − u ( x j ) = µ 0 для всех xi , x j ∈ Ω , n ≥ 1 , µ 0 > 0 .
Определение 5.2.2. Имеет место интенсивность предпочтения степени n
по модулю µ n , если
xi f x j ⇒ u ( n) ( xi ) − u ( n) ( x j ) = µ n для всех xi , x j ∈ Ω , n = 1,2...
(5.2.1)
µ n > 0 n = 2m − 1, m = 1,2,... ,
µ n > 0 n = 2m, m = 1,2,... для рискофилов,
µ n < 0 n = 2m, m = 1,2,... для рискофобов.
Обозначим
u ( n ) ( xi , x j ) = u ( n ) ( xi ) − u ( n ) ( x j )
–
функция
интенсивности
предпочтения степени n по модулю µ n , тогда условие (5.2.1) будет выглядеть
следующим образом:
x i f x j ⇒ u ( n ) ( x i , x j ) = µ n для всех xi , x j ∈ Ω , n = 1,2... .
•
Относительная интенсивность предпочтения
Приведенные выше определения 5.2.1 и 5.2.2 были сформулированы для
абсолютной интенсивности предпочтения. Далее мы рассмотрим относительную
интенсивность предпочтения.
Используя метод естественной нормализации, предложенный в книге [43],
введем определение относительной интенсивности предпочтения, проводя
аналогию с первой интервальной мерой ценности, представленной в работе [19]
Определение 5.2.3. Для множества исходов принятия решений {x j } = Ω и
заданной на них функции полезности u ( x) : Ω → R1 пусть выполняются условия:
‒ Существует исход x1 такой, что для всех j выполняется условие x1 f x j , и
существует исход x r такой, что для всех j выполняется x j f x r .
‒ Неверно соотношение: x r f x1 ,
‒ Выполняется условие x j f x r ⇒ u ( x j ) > u ( x r ) для всех x j ∈ Ω .
73
Тогда относительной интенсивностью предпочтения j–го исхода назовем
величину
R( x j ) =
u( x j ) − u( xr )
u ( x1 ) − u ( x r )
> 0 , j = 2K r − 1 .
Легко заметить, что 0 ≤ R( x j ) ≤ 1 . То есть относительная интенсивность
предпочтения служит интервальной мерой.
Сформулируем
предпочтения,
определение
опираясь
на
степени
результаты
относительной
исследования,
интенсивности
полученные
для
формулировки определения 5.2.1.
Определение
5.2.4.
Имеет
место
относительная
интенсивность
предпочтения j–го исхода степени n, если
R( n) ( x j ) =
u ( n) ( x j ) − u ( n) ( xr )
u ( n) ( x1 ) − u ( n) ( xr )
для всех x j ∈ Ω , j = 2 K r − 1, n = 1,2... (5.2.2)
Зная, что u ( n ) ( xi , x j ) – функция интенсивности предпочтения степени n,
перепишем условие (5.2.2) следующим образом:
R( n ) ( x j ) =
Необходимо
u ( n) ( xi , x r )
u ( n ) ( x1 , x r )
понимать,
что
для всех xi , x j ∈ Ω , n = 1,2...
значение
относительной
интенсивности
предпочтения степени n во многом зависит от вида функции полезности и ее
производной. Поэтому для использования данной характеристики интенсивности
предпочтения, необходимо исследовать поведение функции на рассматриваемом
промежутке и определить наилучшую и наихудшую альтернативу для каждой
степени интенсивности предпочтения.
Тогда мы обеспечим выполнение условия: 0 ≤ R( n ) ( x j ) ≤ 1 , которое, в свою
очередь позволяет определить меру интенсивности предпочтения.
•
Меры интенсивности предпочтения
Функция относительной интенсивности предпочтения степени n принимает
значения в интервале от 0 до 1, разобьем данный промежуток на m интервалов
74
следующим образом: [k1 , k 2 ] ∪ (k 2 , k 3 ] ∪ ... ∪ (k m−1 , k m ] , где k1 = 0, k m = 1, а ki –
значения, заданные экспертом.
Определение 5.2.5. Имеет место мера интенсивности предпочтения степени
n порядка ноль, если R ( n ) ( x j ) ∈ [ k1 , k 2 ] .
Определение 5.2.6. Имеет место мера интенсивности предпочтения степени
n порядка (m − 2) , если R ( n ) ( x j ) ∈ ( k m −1 , k m ] .
•
Ожидаемая интенсивность предпочтения
Один из видов неопределенности при принятии решений характеризуется
вероятностной неопределенностью исходов, подразумевает выбор альтернатив,
которые описываются с помощью вероятностной меры на множестве исходов. В
работе [3] разработано несколько систем аксиом для отношения предпочтения f ,
которое рассматривается на множестве альтернатив, являющихся вероятностными
мерами.
Введем вероятностное пространство (Ω, А, Р), где Ω – множество
элементарных исходов, Α – σ-алгебра, заданная на Ω, а P – вероятностная мера,
определенная на множествах из Α. Случайной величиной X будем называть
произвольное измеримое отображение из Ω в R.
Обозначим
χ
совокупность
всех
случайных
величин
на
(Ω, А).
Перенумеровав элементы Ω некоторым произвольным образом: Ω = {ω1 , K, ω n } ,
обозначим X (ωi ) = xi , i = 1,…n и будем отождествлять случайные величины
X ∈ χ с векторами X = ( x1 , K, xn ) из R n .
Рассмотрим модель ожидаемой интенсивности предпочтения:
ν p κ ⇔ E(u,ν × κ ) < 0,
u( x, y) = −u( y, x) для всех x, y ∈ Ω ,
которая
дает
возможность
практического
шкалирования
сравнительной
полезности и интенсивности предпочтения на основе оценок предпочтения между
вероятностными мерами на множестве исходов [3].
75
Определение 5.2.7. Имеет место ожидаемая интенсивность предпочтения
степени n по модулю µ n , если
ν f κ ⇒ E (u ( n) ,ν × κ ) = E (u (n) ( x) − u ( n) ( y)) = µ n для всех n = 1,2...
µ n > 0 n = 2m − 1, m = 1,2,... ,
µ n > 0 n = 2m, m = 1,2,... для рискофилов,
µ n < 0 n = 2m, m = 1,2,... для рискофобов,
где ν , κ ∈ P – вероятностные меры соответствующие исходам x, y ∈ Ω .
•
Пример применения новых определений
Задача.
ожидаемую
Оценить
интенсивность
интенсивность
предпочтения
предпочтения,
относительную
однокомпонентных
и
медицинских
препаратов на основе ибупрофена и лоперамида в зависимости от их стоимости
для потребителя.
Решение. Существует широкий спектр медицинских препаратов на основе
ибупрофена и лоперамида, производимый различными фармацевтическими
компаниями.
Соответственно
их
цена
варьируется,
в
зависимости
от
производителя и аптеки. Диапазоны цен препаратов с учетом формы выпуска
оказались следующими:
• активное вещество ибупрофен: от 0,25 € до 1,90 €,
• активное вещество лоперамид: от 0,45 € до 4,70 €.
Тогда
множество
исходов
(альтернатив)
для
ибупрофена
есть
Ω I = [0,25; 1,90] , а для лоперамида – Ω L = [0,45; 4,70] .
В качестве лица, принимающего решение, выступает группа студентов.
Ясно, что с ростом цены препарата, его полезность снижается, поэтому функция
полезности будет монотонно убывающей в указанном диапазоне. С помощью
метода, предложенного в книге [41], строятся функции полезности для данных
препаратов. Сначала для каждого препарата находится детерминированный
эквивалент лотереи 50-50 вида <0,25;1,90> и <0,45;4,70>, который равен
соответственно 1,00 и 3,30. Так как по смыслу задачи функция полезности
76
убывающая, то полагаем, что u I (0,25) = u L (0,45) = 1 и u I (1,90) = u L (4,70) = 0 , где
u I ( x) – функция полезности для препарата с ибупрофеном,
u L (x) – функция
полезности для препарата с лоперамидом. Далее с использованием полученных
значений, находятся еще три детерминированных эквивалента. Полученные
данные приведены в таблицах ниже:
u I ( x)
x
0
2,00
0,25
1,50
0,325
1,25
0,5
0,80
0,75
0,50
1
0,25
Табл. 5.2.1. Таблица значений функции полезности для ибупрофена.
u L (x)
x
0
4,70
0,25
4,30
0,5
3,30
0,75
0,85
0,875
0,60
1
0,45
Табл. 5.2.2. Таблица значений функции полезности для лоперамида.
Затем проводится проверка согласованности полученных результатов, и
строятся кривые функций полезности для ибупрофена и лоперамида. Далее с
помощью программы «Advanced Grapher» были вычислены аппроксимирующие
функции, проходящие через полученные точки:
• u I ( x) = −0,4628 ⋅ ln(x) + 0,3952 , с коэффициентом детерминации R 2 = 0,98 ;
• u L ( x) = 4,6880 ⋅10−4 ⋅ x 4 - 0,6332 ⋅ x3 + 0,4671 ⋅ x 2 - 1,1477 ⋅ x + 1,4210 ,
с коэффициентом детерминации R 2 = 0,99 .
1. Найдем значения функции полезности для набора альтернатив:
xi
u ( x)
u I ( x)
u L (x)
x1 = 0,80
x2 = 1,00
x3 = 1,20
x4 = 1,40
x5 = 1,60
0,4985
0,7696
0,3952
0,6776
0,3108
0,6080
0,2395
0,5579
0,1777
0,5243
Табл. 5.2.3. Значения функций полезности для набора альтернатив
{ , 8 , 7 , S , T }={0,80; 1,00;1,20;1,40;1,60}.
2. Используя определение 5.2.1, найдем интенсивность предпочтения для пар
альтернатив:
77
xi f x j
u ( xi , x j )
u I ( xi , x j )
u L ( xi , x j )
x1 f x2
x 2 f x3
x3 f x4
x4 f x5
x1 f x3
x2 f x4
x3 f x5
0,1033
0,0920
0,0844
0,0696
0,0713
0,0501
0,0618
0,0336
0,1876
0,1616
0,1557
0,1197
0,1331
0,0837
Табл. 5.2.4. Значения функций интенсивности предпочтения для набора альтернатив
{ , 8 , 7 , S , T }={0,80; 1,00;1,20;1,40;1,60}.
По
полученным
данным
можем
сделать
вывод,
что
наибольшая
интенсивность предпочтения для обоих препаратов приходится на альтернативы
x1 и x 2 , а наименьшая – между альтернативами x 4 и x5 . Это говорит о том, что
разница в между двумя низкими ценами для потребителя более важна, чем между
двумя высокими ценами. Причем, для лоперамида интенсивность между
альтернативами x 2 и x3 почти в 2 раза больше, чем для x 4 и x5 . Можно сделать
следующий вывод: чем выше цена, тем безразличнее относится к ее увеличению
ЛПР.
3. Функции полезности u I ( x) и u L (x) являются убывающими на промежутках
Ω I и Ω L соответственно, поэтому производная функции полезности меньше
нуля. Используя определение 5.2.2, найдем первую и вторую степени
интенсивности предпочтения:
xi f x j
u (1) ( x i , x j )
u I (1) ( x i , x j )
u L (1) ( x i , x j )
x2 f x4
x3 f x5
-0,1157 -0,0771 -0,0551 -0,0413 -0,1928 -0,1322
-0,1194 -0,1046 -0,0900 -0,0754 -0,2240 -0,1946
-0,0964
-0,1654
x1 f x2
x 2 f x3
x3 f x4
x4 f x5
x1 f x3
Табл. 5.2.5. Значения функций интенсивности предпочтения первой степени для набора
альтернатив { , 8, 7 , S , T }={0,80; 1,00;1,20;1,40;1,60}.
В нашем примере, для ЛПР, чем меньше цена, тем выше полезность.
Соответственно,
интенсивность
предпочтения
первой
степени
можно
интерпретировать следующим образом: увеличение цены на 0,2€ в диапазоне
низких цен снижает полезность для ЛПР значительно больше, чем увеличение на
ту же сумму в диапазоне высоких цен. То есть, чем ниже цена, тем сильнее
реагирует покупатель на ее повышение.
78
xi f x j
u ( 2) ( xi , x j )
u I ( 2) ( xi , x j )
u L ( 2) ( xi , x j )
x1 f x2
x 2 f x3
x3 f x4
x4 f x5
x1 f x3
x2 f x4
x3 f x5
0,2603
0,0740
0,1414
0,0735
0,0853
0,0731
0,0553
0,0726
0,4017
0,1475
0,2267
0,1466
0,1406
0,1457
Табл. 5.2.6. Значения интенсивности предпочтения второй степени для набора альтернатив
{ , 8 , 7 , S , T }={0,80; 1,00;1,20;1,40;1,60}.
Как известно, знак второй производной функции полезности определяет
отношение к риску. В нашем примере ЛПР по отношению к ценам на оба
препарата является рискофилом. Кроме того чем больше абсолютное значение
второй производной, тем больше ЛПР рискует.
Для двух рассматриваемых активных веществ склонность к риску у ЛПР
уменьшается по мере роста цены препаратов, их содержащих. Но в случае, с
препаратами, содержащими ибупрофен, значения функции интенсивности второй
степени указывают нам на то, что для наиболее предпочтительных альтернатив
ЛПР склонен к риску в большей степени. Тогда как для препаратов на основе
лоперамида, интенсивность предпочтения второй степени показывает, что ЛПР
рискует в одинаковой степени, вне зависимости от того, является ли альтернатива
наиболее предпочтительной или наименее предпочтительной.
4. Найдем относительную интенсивность предпочтения для рассматриваемых
препаратов, используя определение 5.2.3.
xj
R( x j )
RI ( x j )
RL ( x j )
x2 = 1,00
x3 = 1,20
x4 = 1,40
0,3333
0,5828
0,2457
0,5229
0,1715
0,4797
Табл. 5.2.7. Значения функции относительной интенсивности предпочтения для набора
альтернатив { 8 , 7 , S }={1,00;1,20;1,40}.
Величина
относительной
интенсивности
предпочтения
характеризует
насколько текущая альтернатива лучше наихудшей альтернативы относительно
разницы между наилучшей и наихудшей альтернативой. Такая характеристика
позволяет сравнить интенсивность каждой отдельно взятой альтернативы. В
нашем случае значения относительной интенсивности альтернатив указывают на
79
следующее: для ибупрофена альтернатива x 2 почти в два раза предпочтительнее
по интенсивности, чем альтернатива x 4 , для лоперамида разница в интенсивности
предпочтения одиночных альтернатив примерно одинакова.
5. Найдем относительную интенсивность предпочтения первой и второй
степени для рассматриваемых препаратов, используя определение 5.2.4.
xj
R (1) ( x j )
R I (1) ( x j )
R L (1) ( x j )
x2 = 1,00
x3 = 1,20
x4 = 1,40
0,1429
0,5248
0,0952
0,3883
0,0612
0,2708
Табл. 5.2.8. Значения функции относительной интенсивности предпочтения первой степени
для набора альтернатив { 8 , 7 , S }={1,00;1,20;1,40}.
При расчете относительной интенсивности R L (1) ( x j ) мы должны учесть, что
функция u ′L (x) на промежутке
[0,45; 4,70]
вогнутая. Максимальное значение
функции в точке (2,97; 0,51). Тогда, будем исследовать R L (1) ( x j ) на промежутке
[0,45; 2,97], в который попадают рассматриваемые альтернативы.
Исходя из полученных значений R (1) ( x j ) , с учетом свойств первой
производной функции полезности можно сделать вывод, что увеличение цены
почти не влияет на общую тенденцию интенсивности предпочтения альтернатив
для ибупрофена. Для препаратов, имеющих в составе лоперамид, следует
заметить, что увеличение цены повышает разницу в интенсивности предпочтения
альтернатив.
xj
R( 2) ( x j )
R I (2) ( x j )
R L ( 2) ( x j )
x2 = 1,00
x3 = 1,20
x4 = 1,40
0,0476
0,8629
0,0282
0,8136
0,0165
0,7646
Табл. 5.2.9. Таблица значений относительной интенсивности второй степени для набора
альтернатив { 8 , 7 , S }={1,00;1,20;1,40}.
Выводы
в
случае
R( 2) ( x j )
аналогичны
интенсивности предпочтения второй степени.
80
полученным
ранее
для
6. Используя определение 5.2.6, найдем меру интенсивности предпочтения.
Зададим значения ki , которые разбивают промежуток [0;1]
на некоторое
количество интервалов. Разделим промежуток на пять интервалов и найдем
соответствующие меры интенсивности предпочтения ЛПР. Получим следующие
интервалы: [0; 0,2] , (0,2; 0,4] , (0,4; 0,6] , (0,6; 0,8] , (0,8;1] . Далее определим меру
интенсивности предпочтения различных степеней.
Для препаратов с ибупрофеном альтернативы x 2 и x3 имеют меру
интенсивности предпочтения (нулевой степени) первого порядка, а альтернатива
x 4 – порядка ноль. Мера интенсивности предпочтения первой и второй степени
порядка ноль характерна для всех рассматриваемых альтернатив.
Исследование меры для препаратов с лоперамидом дают следующий
результат: все альтернативы имеют меру интенсивности предпочтения третьего
порядка. Интересно, что альтернативы x3 и x 4 имеют меру интенсивности
предпочтения первой степени второго порядка, а альтернатива
x 2 – третьего
порядка.
Чем выше порядок меры интенсивности предпочтения, тем больше
относительная интенсивность предпочтения альтернативы. Меры дают нам
возможность группировать альтернативы по величине интенсивности. Что
упрощает для ЛПР принятие решения, оно может рассматривать не все
возможные альтернативы, а только те, которые имеют меру определенного
порядка, удовлетворяющего ЛПР.
7. Найдем ожидаемую интенсивность предпочтения различных степеней.
Для начала определим интервалы, на которых хотим исследовать ожидаемую
полезность. Для этого воспользуемся предыдущими рассуждениями и выберем
наиболее интересные для исследования альтернативы.
Для исследования ожидаемой интенсивности предпочтения первой степени
возьмем для рассмотрения пары альтернатив: 0,80 и 1,00; 1,00 и 1,20.
81
Воспользуемся определением 5.2.7 для нахождения значения ожидаемой
интенсивности предпочтения первой и второй степени, получим следующие
результаты:
[ xi ; x j ]
E (u (1) ,ν × κ )
E (u I (1) ,ν × κ )
E (u L ( n ) ,ν × κ )
[0,80; 1,00]
[1,00; 1,20]
-0,0926
-0,0762
-0,0926
-0,0824
Табл. 5.2.10. Таблица значений ожидаемой интенсивности предпочтения первой степени
[ xi ; x j ]
E (u ( 2 ) ,ν × κ )
E (u I ( 2) ,ν × κ )
E (u L ( 2 ) ,ν × κ )
[0,80; 1,00]
[1,00; 1,20]
0,0844
0,1149
0,1033
0,1072
Табл. 5.2.11. Таблица значений ожидаемой интенсивности предпочтения второй степени
Результаты: введенные в данной работе понятия позволяют наиболее
тщательно и глубоко изучить интенсивность предпочтения ЛПР. Степени
интенсивности предпочтения позволяют оценить соотношение интенсивности
предпочтения и степень реакции ЛПР на увеличение / уменьшение благ, а также
влияние отношения к риску ЛПР на интенсивность. Понятие меры интенсивности
позволяют сгруппировать альтернативы по уровню интенсивности и отсеять
альтернативы с уровнем, не удовлетворяющим интересы исследователя.
Важно заметить, что использование введенных определений зависит от вида
функции полезности, не всякая функция полезности будет удобна для проведения
глубокого исследования предпочтений. В частности, полиномиальная функция
полезности нечетного порядка для исследований степеней интенсивности не
годится, так как на промежутке исследования ее производные ведут себя
непредсказуемо.
Результаты данного исследования приведены в статье [44].
82
5.3. Порог принятия решения
Функция интенсивности предпочтения (определение 5.1.1), предполагает
существование монотонной функции полезности. На практике не всегда удается
задать монотонную функцию интенсивности предпочтения на всем множестве
альтернатив. Поэтому далее мы рассмотрим несколько немонотонных функций
полезности и оценку интенсивности предпочтения для них.
Примером случая, когда функция полезности бывает немонотонной, может
служить задача из области медицины. Например, исследование результатов
анализа крови на содержание различных компонентов: гормонов, сахара,
лейкоцитов и т.д. Для каждого компонента существует своя «норма», и
отклонение в большую или меньшую сторону, означает отрицательный результат.
В этом случае, функция полезности при содержании компонента меньше нормы
возрастать, а больше нормы – убывать.
По смыслу своего применения функция полезности непрерывная и кусочномонотонная (немонотонна) на множестве альтернатив, то есть она является
возрастающей или убывающей на интервалах промежутка ее определения.
Рассмотрим функцию полезности ( ): L ; VM → , где промежуток L ; VM можно
разбить на k интервалов, на которых функция или убывает, или возрастает. В
зависимости от способа построения функции полезности ЛПР может изначально
знать границы промежутков монотонности функции на L ; VM, а может для ранее
выбранной функции полезности определить промежутки ее монотонности.
Теорема
5.3.1.
Непрерывная
кусочно-монотонная
функция
( )
определенная на промежутке L ; VM имеет абсолютный максимум и абсолютный
минимум на этом промежутке.
Доказательство:
Так как функция непрерывна на промежутке La; bM, значит, по теореме
Вейерштрасса (о непрерывности функций), она ограничена на нем и достигает
своих минимального и максимального значений.
Если функция непрерывна и монотонна на отрезке, то она достигает своих
максимального и минимального значений в концах отрезка.
83
Функция ( ) кусочно-монотонна и непрерывна на L ; VM и, следовательно,
,
существует конечное число таких точек
a<
8 … X,
< x8 < ⋯ <
и ( ) монотонна в каждом из промежутков L ;
Y
что
<V
M,L ;
8M
, …, L
X ; VM.
Следовательно, максимальным значением функции ( ) будет наибольшее
из значений
( ), ( ), ( 8 ) , … , (
Минимальное значение
значений.■
X ),
(V).
( ) будет соответствовать наименьшему из этих
Мы выяснили, что у функции полезности может быть несколько
максимумов и минимумов. Так как данные значения являются порогом, при
переходе через который функция меняет свой характер: из возрастающей
становится убывающей и наоборот, они представляют особый интерес для ЛПР.
Исследуя интенсивность предпочтения, ЛПР при достижении экстремума
функции полезности должен производить замену функции интенсивности
предпочтения.
Определение 5.3.1. Порогом принятия решения для функции полезности
( ) определенной на L ; VM будем называть такую альтернативу
max ( ) или ( ∗ ) = min ( ) для любых Z L ; VM.
∗
, что ( ∗ ) =
Определение 5.3.2. Абсолютным (максимальным или минимальным)
порогом принятия решения для функции полезности ( ) будем называть такую
альтернативу
∗
Z L ; VM.
[ , что
( [) = max ( ∗ ) или
( [) = min ( ∗ ) для любых
5.4. Порог принятия решения на практике
Постановка задачи: владелец
аптеки, зная статистические данные о
заболеваемости пневмонией, опубликованные на сайте Роспотребнадзора РФ (см.
Приложение, табл. 5) [45], хочет выяснить, в какое время нужно увеличить
84
закупаемое количество препаратов для лечения пневмонии с целью повышения
прибыли.
Для решения данной задачи построим функцию полезности, определим
порог принятия решения и вычислим интенсивность предпочтения.
Решение: ясно, что потребность в препаратах для лечения пневмонии
напрямую зависит от количества заболевших. Следовательно, при построении
функции полезности ЛПР поставит максимальное значение полезности в
соответствие месяцу, когда заболевших больше всего. По индексу сезонности
заболеваемости понятно, что ноябрь и февраль – месяцы с повышенным
количеством случаев заболевания, а август – месяц, с наименьшим показателем
заболеваемости. На эти месяцы следует обратить особое внимание, но
необходимо также понять: какова потребность в медикаментах, а соответственно,
возможность получить наибольшую прибыль вблизи пикового значения.
Функция полезности строилась, исходя из среднего количества больных в
каждом месяце, с помощью метода, уже использованного в пункте 5.2 данной
работы. Полезность в нашем примере выражается не в денежном эквиваленте, а в
баллах (условных единицах), заданных ЛПР как оценка прибыли на конец месяца,
величина которой составляет от 15 до 150 баллов.
#
1
2
3
4 5 6 7 8 9 10 11 12
( ( # ) 125 150 120 90 50 30 20 15 40 85 100 75
∗
Табл. 5.4.1. Значения функции полезности.
С помощью программы Advanced Grapher 2.2 (методом регрессионного
анализа) построена аппроксимирующая полиномиальная функция полезности,
которая имеет вид:
(( ) = 5,1724 ∗ 10`S
a
− 107,6842
Коэффициент корреляции:
− 0,0256
7
8
c
+ 0,5190
+ 159,3187
= 0,9973
85
8
e
− 5,5287
T
+ 33,0775
− 54,8592 + 100,3408
S
Рис. 5.4.1. График функции полезности
Полученная нами функция является кусочно-монотонной и имеет четыре
промежутка монотонности.
С помощью вышеназванной программы проведено исследование функции
полезности на экстремумы, на промежутке
1,91,
∈ L1; 12M их три:
Экстремум
x
(( )
max
1,91
149,03
min
7,53
13,55
max
10,86 102,19
Таким образом, мы нашли три порога принятия решения для ЛПР:
∗
8
= 7,53,
∗
7
∗
=
= 10,86, один из этих порогов является абсолютным
максимальным: [hij = 1,91, а другой – абсолютным минимальным: [h#! = 7,53.
Рассмотрим интенсивность предпочтения альтернатив, близких к порогам
принятия решения, которые интересуют ЛПР, а именно пороги с максимальными
значениями функции полезности.
1. Абсолютный
максимальный
порог
принятия
решения
является
альтернативой [hij = 1,91. Так как множество альтернатив – временной
промежуток – 12 месяцев, то [hij = 1,91 обозначает, в данном контексте, 25
февраля (мы знаем, что промежуток от 1 до 2 – это февраль, в этом месяце 28
дней, следовательно: 0,91 часть от 28 дней это 25,48 ≈ 25 февраля). Теперь
рассмотрим альтернативы, соответствующие 22 февраля:
86
∈ L1,78; 1,81M, 23
февраля:
∈ L1,83; 1,85M, 24 февраля:
рассмотрим альтернативы со значениями
88
∈ L1,86; 1,89M. Для определенности
= 1,81,
Интенсивность предпочтения определяется как:
Тогда
(
функция
#`
,
#)
= (( # ) − ((
интенсивности
87
#`
предпочтения
= 1,85,
).
для
8S
= 1,89.
рассматриваемых
альтернатив будет иметь следующие значения:
xi f x j
U ( xi , x j )
x 25 f x 24
x 24 f x 23
x 23 f x22
0,0096
0,0938
0,1867
U ( xi , x j )
Табл. 5.4.2. Значения интенсивности предпочтения для альтернатив, соответствующих 22-25
февраля
Легко заметить, что интенсивность предпочтения за 3 дня до получения
максимальной прибыли (22 февраля), значительно выше, чем за день до этого
события. То есть ЛПР следует позаботиться о том, чтобы сделать запасы
медикаментов не только заранее, но и не позднее, чем за 3 дня до наступления
максимального спроса на препараты.
Теперь рассмотрим период времени после пикового значения полезности и
возьмем для рассмотрения 26-28 февраля, по аналогии с предыдущими
вычислениями, 26 февраля:
8a
8e
= 1,96, 27 февраля:
8c
= 1,99, 28 февраля:
= 2,03. При этом учитывая, что функция полезности убывает после перехода
абсолютного
порога
предпочтения
для
принятия
данного
решения,
промежутка
т.е.
функция
интенсивности
будет
другой.
Интенсивность
предпочтения для альтернатив февраля будет определяться так: (
((
#`
) − (( # ).
#`
,
#)
=
Величина функции интенсивности предпочтения для рассматриваемых
альтернатив равна:
xi f x j
U ( xi , x j )
U ( xi , x j )
x25 f x26
x 26 f x 27
x 27 f x28
0,0642
0,1067
0,2181
Табл. 5.4.3. Значения интенсивности предпочтения для альтернатив, соответствующих 25-28
февраля
87
Получаем, что интенсивность предпочтения для 28 февраля значительно
выше, чем для 26 февраля, то есть прибыль через 3 дня начнет убывать с большей
скоростью, чем через день-два после максимального значения. Кроме того,
интенсивность предпочтения в период 26-28 февраля выше, чем в период 22-24
февраля, из этого можно сделать вывод скорость уменьшения прибыли после
пикового значения, будет больше, чем ее рост до пикового значения. Вероятно, в
период после 3 дней следует провести акцию по снижению стоимости
медикаментов, с целью увеличения количества продаваемого товара.
∗
7
2. Второй максимальный порог принятия решения является альтернативой
= 10,86, соответствующей 25 ноября. Аналогично предыдущим рассуждениям
исследуем альтернативы, соответствующие 22-24 ноября:
10,79,
8S
= 10,83 и 26-28 ноября:
8e
= 10,89,
Получим следующие интенсивности предпочтения:
xi f x j
U ( xi , x j )
U ( xi , x j )
8c
88
= 10,76,
= 10,93,
x 25 f x 24
x 24 f x 23
x 23 f x22
0,0256
0,1148
0,1462
8a
87
=
= 10,96.
Табл. 5.4.4. Значения интенсивности предпочтения для альтернатив, соответствующих 22-25
ноября
xi f x j
U ( xi , x j )
U ( xi , x j )
x25 f x26
x 26 f x 27
x 27 f x28
0,0265
0,1166
0,1484
Табл. 5.4.5. Значения интенсивности предпочтения для альтернатив, соответствующих 25-28
ноября
Легко заметить, что в обоих случаях интенсивность предпочтения соседних
альтернатив увеличивается при отдалении от точки порога принятия решения,
причем почти симметрично. В данном временном периоде рост и уменьшение
прибыли происходит без резких скачков, следовательно, не стоит принимать
специальных мер (например, скидочных акций) для привлечения клиентов.
88
Результаты
В первой главе проведено исследование различных мер риска на
когерентность. Исходя из условий когерентности, можно сделать вывод, что
создание
меры
риска
универсальной
для
порядков
стохастического
доминирования невозможно. Потому что мера должна содержать в себе
реализацию случайной величины в соответствующей степени, как, например,
LPM n,t ( X ) , чтобы удовлетворять разным порядкам доминирования.
После построения мер риска для статистических данных, можно сказать,
что наиболее удобной когерентной мерой является ESα ( X ) , рассчитанная по
формуле (1.4.2). Из некогерентных мер риска можно выделить две меры:
LPM n ,t ( X ) и VARα ( X ) . Мера VARα ( X ) дает результаты, схожие с ESα ( X ) , но с
большим запасом, когда применяется к эмпирическим данным. А мера LPM n ,t ( X )
интересна тем, что в зависимости от степени n мера соответствует разным
порядкам стохастического доминирования, разным порядкам риска.
Во второй главе оказана теорема, связывающая понятия степени рисковости
и порядок неприятия риска. Это позволяет использовать метод ожидаемой
полезности в терминах степеней рисковости. В третьей главе доказаны теоремы,
описывающие взаимосвязь мер риска и порядка стохастического доминирования.
Решена задача по сравнению результатов методов оценки риска в условиях
рисковости разных степеней. Можно сделать вывод о том, что использование
степеней рисковости дает более гибкие и подробные результаты, нежели
исследование данных без применения этого понятия. А так же можно заметить,
что аппарат мер риска является более гибким, чем методы ожидаемой полезности
и стохастического доминирования.
Таким образом, итоги сравнения альтернатив зависят и от метода,
применяемого для их ранжирования, и от степени рисковости, выбранной
экспертом.
В пятой главе введены дополнительные измерения интенсивности
предпочтения: степени и меры, применение которых рассмотрено на практике.
89
Также предложен порог принятия решения для кусочно-монотонных функций
полезности, доказана теорема о существовании порога и приведен пример его
применения. Оценка интенсивности предпочтений позволяет проанализировать
близкие по величине полезности альтернативы, сделать выводы о скорости
увеличения или уменьшения полезности при выборе из множества наиболее
предпочтительных альтернатив.
В работе исследованы методы, решающие задачу оптимизации и
принадлежащие
различным
математическим
аппаратам.
Продолжить
это
исследование можно в области мер риска, для того, чтобы создать меру, которая
позволяла бы сравнивать результаты метода ожидаемой полезности и принципа
стохастического доминирования между собой. Так как на практике лица
принимающие
решения
сталкиваются
с
ситуациями,
когда
необходимо
использовать оба метода и иметь возможность сравнить результаты упорядочения
обоих методов.
90
Заключение
Современная теория решений постоянно развивается, это подтверждают
многочисленные исследования различных подходов к оценке риска. В литературе
подробно описана связь метода ожидаемой полезности и некоторых мер риска с
принципом стохастического доминирования. В данной работе исследованы три
метода количественной оценки риска: аппарат мер риска, метод ожидаемой
полезности и принцип стохастического доминирования.
В первой главе проведен анализ мер риска, используемых в теории
принятия решений на когерентность. Исследовано их практическое применение к
эмпирическим данным, выявлены наиболее предпочтительные меры рисков для
проведения количественной оценки риска для статистических данных.
Во второй главе проведено исследование условий, которые накладываются
на функцию полезности при совместном ее использовании с мерами риска,
доказаны теоремы 2.4.1-2.4.3. Так же изучена и доказана связь степеней
рисковости с методом ожидаемой полезности, а именно какое соотношение
связывает степени рисковости с производными функции полезности, см. теоремы
2.5.1,
2.5.2,
2.5.4.
В
третьей
главе
описан
принцип
стохастического
доминирования и его связь с мерами риска, доказаны теоремы 3.2.1 и 3.2.2 о
соотношении порядков стохастического доминирования и мер риска.
В четвертой главе решается задача сравнения результатов трех методов
оценки риска с использованием понятия рисковости различных степеней и без
использования данного определения.
В пятой главе исследован метод оценки интенсивности предпочтения,
предложены новые инструменты для измерения интенсивности (введены
определения 5.2.1 – 5.2.7), результат применения которых проверен на примере.
Также предложено понятие порога принятия решений (введены определения
5.3.1, 5.3.2) для кусочно-монотонных функций полезности, доказана теорема
5.3.1, проведена интерпретация результатов на примере.
91
Литература
1.
Бернулли Д. Опыт новой теории измерения жребия [1738] // Вехи
экономической мысли. Т.1. Теория потребительского поведения и спроса.
СПб.: Экономическая школа, 1999, С.11–27.
2.
Дж. фон Нейман, О. Моргенштерн. Теория игр и экономическое поведение. /
Перев. с англ. под ред. и с доб. Н. Н. Воробьева. М.: Наука, 1970. 708 с.
3.
Фишберн П. Теория полезности для принятия решений. М.: Наука, 1978.
352 с.
4.
Колбин В.В., Белоносова И.Ю. Анализ рисков в страховании, СПб.:
Издательство СПбГУ, 2006. 68 с.
5.
Крувшиц Л., Шефер Д., Шваке М. Финансирование и инвестиции. Сборник
задач и решений / Пер. с нем. под общей редакцией Сабова З. А., Дмитриева
А. Л. СПб.: Питер, 2001. 320 с.
6.
Новоселов А. А. Стохастическое доминирование и его приложения в
моделировании риска // Записки ФАМ Семинара, Красноярск, 2002. Т. 7.
С. 37–44.
7.
Новоселов А. А., Варочкина Т. С. Стохастическое доминирование I и II
рода // Вестник Красноярского государственного университета, 2004. Т. 5.
Вып. 2. С. 15–21.
8.
Ekern S. Increasing Nth degree risk // Economics Letters 6, 1980. P. 329-333.
9.
Jean W. H. The geometric mean and stochastic dominance // Journal of Finance
35, 1980. P. 151-158.
10. Fishburn P. C. Mean-risk analysis with risk associated with below-target returns //
The American Economic Review, 1977. Vol. 67, 2. P. 116–126.
11. Artzner P., Delbaen F., Eber J. M. Coherent measures of risk //Math. Fin., 1999.
Vol. 9, 3. P. 203–228.
12. Barbosa A., Ferreira M. A. Beyond coherence and extreme losses: Root Lower
Partial
Moment
as
a
risk
measure,
2004.
Available
at
http://ssrn.com/abstract=609221 or http://dx.doi.org/10.2139/ssrn.609221.
92
SSRN:
13. Peng J. Value at Risk and Tail Value at Risk in uncertain environment //
Proceedings of the Eighth International Conference on Information and
Management Sciences, 2009. P. 787–793.
14. Rootzén H., Kluppelberg C. A single number can’t hedge against economic
catastrophes // Ambio, 1999. Vol. 28, No. 6. P. 550–555.
15. Danielsson J., Jorgensen B. N., Sarma M., Samorodnitsky G., de Vries C. G. Fat
tails, VaR and subadditivity // Journal of Econometrics, 2013. Vol. 172, 2. P. 283–
291.
16. Benati S. The optimal portfolio problem with coherent risk measure constraints //
Europian Journal of Operational Research, 2003. Vol. 150, No. 3. P. 572–584.
17. Кирута
А.Я.,
Рубинов
А.М.,
Яновская
Е.Б.
Оптимальный
выбор
распределений в сложных социально-экономических задачах. Л.: Наука,
1980. 167 c.
18. Peter H. Farquhar, L. Robin Keller Preference intensity measurement //Annals of
Operations Research, 1989, Vol. 19, 1, P. 205–217.
19. Фишберн П. К. Измерение относительных ценностей / Статистическое
измерение качественных характеристик. – М.: Статистика (Перевод с
английского под редакцией Е. М. Четыркина), 1972. С. 34–94.
20. Алехин Е.И. Теория риска и моделирование рисковых ситуаций. Конспект
лекций. Орел: ОГУ, 2006. 90 с.
21. Панягина А.Е. Подходы к пониманию и классификации рисков //
Современная экономика: проблемы, тенденции, перспективы. 2012. Вып. №6.
С. 1–11.
22. Найт Ф. Понятие риска и неопределенности // Thesis, 1994. Вып. 5. C. 12–28
(Перев. с англ. С. А. Афонцева из: Frank H. Knight. The Meaning of Risk and
Uncertainty // F. Knight. Risk, Uncertainty, and Profit. Boston: Houghton Mifflin
Co, 1921. P. 210–235).
23. Новоселов А. А. Математическое моделирование финансовых рисков: теория
измерения. Новосибирск: Наука, 2001. 99 c.
93
24. Новоселов А. А. Обобщенные когерентные меры риска // Труды IV
Всероссийской
конференции
по
финансово-актуарной
математике
и
смежным вопросам. Т. 1. Красноярск, 2005. С. 325–339.
25. Acerbi C., Tasche D. On the coherence of Expected Shortfall // Journal of Banking
& Finance, 2002. Vol. 26(7). P. 1487–1503.
26. Ледовская, В.А. Анализ некогерентных мер риска // Процессы управления и
устойчивость: Труды 44-й международной научной конференции аспирантов
и студентов / СПб.: Издат. Дом СПбГУ, 2013. С. 605–610.
27. Ледовская, В.А. Колбин В.В. Когерентные меры риска // Процессы
управления
и
устойчивость:
Труды
44-й
международной
научной
конференции аспирантов и студентов / СПб.: Издат. Дом СПбГУ, 2013.
С. 611–616.
28. Ледовская В.А, Колбин, В.В. Исследование мер риска в теории принятия
решений / В.В. Колбин, // Оценка и управление риском / Palmarium Academic
Publishing, 2014. С. 191–245.
29. The United States Geological Survey. http://www.usgs.gov/natural_hazards/
30. Statistic Brain. http://www.statisticbrain.com/
31. Дубров А.М., Лагоша Б.А., Хрусталев Е.Ю.
Моделирование рисковых
ситуаций в экономике и бизнесе. / Под ред. Б.А. Лагоши. М.: Финансы и
статистика, 2000. 176 с.
32. Pratt J.W. Risk Aversion in the Small and in the Large // Econometrica, January–
April 1964, Vol. 32, P. 122–136.
33. Arrow K.J. Essays in the Theory of Risk-Bearing. Chicago: Markham Pub. Co.,
1971. 278 p.
34. Levy H., Kroll Y. Ordering uncertain options with borrowing and lending // The
Journal of Finance, 1978. Vol. 33, 2. P. 553–574.
35. Колбин В. В. Анализ и оценка рисков. СПб: Изд. дом СПбГУ, 2005. 70 c.
36. Буренин А.Н. Управление портфелем ценных бумаг. М., Научно-техническое
общество имени академика С.И. Вавилова, 2008. 440 с.
94
37. Follmer. H., Schied, A. Stochastic Finance: An Introduction in discrete Time, 2nd
rev. and extended ed. De Gruyter Studies in Mathematics; 27. 2004. P. 459.
38. Porter R. Semivariance and Stochastic Dominance: A Comparison // The
American Economic Review, 1974. Vol. 64, 1. P. 200–204.
39. Европейское региональное бюро Всемирной организации здравоохранения.
http://data.euro.who.int/hfadb/shell_ru.html
40. Jean W. H., Helms B. P. Stochastic dominance as a decision model // Quarterly
Journal of Business and Economics, 1986. Vol. 25, 1. P. 65–101.
41. Кини Р. Л., Райфа X. Принятие решений при многих критериях:
предпочтения и замещения / Пер. с англ. под ред. И. Ф. Шахнова. М.: Радио и
связь, 1981. 560 с.
42. Deck C., Schlesinger H. Consistency of Higher Order Risk Preferences //
Econometrica, Econometric Society, 2014. Vol. 82. P. 1913–1943.
43. Трухачев Р.И. Модели принятия решений в условиях риска и не
определенности. М.: Наука, 1981. 151 с.
44. Ledovskaya, V.A., Kolbin V.V. Degrees and measures of preference intensity //
International Conference on «Stability and Control Processes» in Memory of V.I.
Zubov, SCP / IEEE, Article number 7342192, 2015. P. 470–474.
http://ieeexplore.ieee.org/document/7342192/.
45. Роспотребнадзор. Статистические материалы.
http://www.rospotrebnadzor.ru/activities/statistical-materials/
95
Приложение
Табл. 1. Исходные статистические данные по количеству погибших во время
природных катастроф.
Количество погибших, чел.
№
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
Землетрясение
1000
1000
1000
1000
1117
1185
1313
1342
1500
1567
1621
1989
2000
2183
2200
2266
2323
2400
2500
2735
2800
3000
4000
5000
5502
5749
9500
9748
17118
20085
Извержение
вулкана
9
31
43
150
152
153
245
353
500
700
700
700
800
1000
1184
1200
1335
1335
1369
1377
1475
1680
2000
2000
2942
2957
3360
3500
4011
5110
96
Наводнение
Цунами
1525
3034
3247
3841
4170
4421
4554
5624
5951
6229
6401
7453
7625
7829
7890
7992
8720
9838
10140
10698
12502
21621
23896
33780
107485
150180
173123
119
1700
1997
2000
2144
2182
2243
3000
3000
3000
3620
3800
5000
5000
5233
8000
13486
15030
25674
26360
30000
36500
40000
Количество погибших, чел.
№
Землетрясение
31
32
33
34
35
36
37
38
20896
25000
31000
50000
86000
87587
227898
316000
Извержение
вулкана
6000
9350
14300
25000
29025
36417
92000
Наводнение
Цунами
Табл. 2. Промежуточные вычисления для меры β-риск.
Катастрофа
Наводнение
Землетрясение
Извержение
Цунами
Стандартное
отклонение
43762,83
62437,79
16395
12032,36
Коэффициент
корреляции
0,96319599
0,84927857
0,79124502
0,99086903
Табл. 3. Статистические данные о количестве больных ВИЧ и туберкулезом в
России за период с 1994 по 2012 г.г..
Год
1 994
1 995
1 996
1 997
1 998
1 999
2 000
2 001
2 002
2 003
2 004
Количество инфицированных
ВИЧ
887
1090
2603
6918
10889
30647
89808
177579
227502
263898
296045
97
Количество больных
туберкулезом
94214
101570
115792
118411
150879
128871
165401
171352
163618
160130
154344
Год
2 005
2 006
2 007
2 008
2 009
2 010
2 011
2 012
Количество инфицированных
ВИЧ
300252
373718
403100
471676
530185
589581
650100
703781
Количество больных
туберкулезом
150904
157649
154964
157859
159110
157110
146762
136293
Табл. 4. Значения функции полезности для случайных величин c функциями
распределения F ( x) и G(x) .
Соответствующее значение
функции распределения
U F ( xi )
U G ( xi )
Интервал [0; 0,2]
0,05
0,11
0,16
449,9158
541,6065
1 185,5688
29 971,1553
32 069,1383
36 083,5436
Интервал (0,2; 0,4]
0,21
0,26
0,32
0,37
2 857,4754
4 298,2255
10 908,1015
28 706,7166
36 817,4889
39 731,9676
41 785,5515
44 663,3123
Интервал (0,4; 0,6]
0,42
0,47
0,53
0,58
53 021,5532
66 265,3807
75 734,2779
83 988,9124
45 789,5918
45 796,2645
46 734,8036
46 903,8446
Интервал (0,6; 0,8]
0,63
0,68
0,74
0,79
85 062,3385
103 583,2694
110 884,6726
127 726,0263
98
47 487,8302
47 634,5317
47 691,7115
48 031,7153
Соответствующее значение
функции распределения
U F ( xi )
U G ( xi )
Интервал (0,8;1]
0,84
0,89
0,95
1,00
141 900,7506
156 130,9964
170 483,4202
183 102,2935
48 308,6386
49 254,5978
49 737,4472
51 345,3024
Табл. 5. Статистические данные по месяцам за 2013-2016 г.г. о количестве
больных пневмонией (внебольничной) в России.
№
Месяц
2013
2014
2015
2016
1
2
3
4
5
6
7
8
9
10
11
12
январь
февраль
март
апрель
май
июнь
июль
август
сентябрь
октябрь
ноябрь
декабрь
68314
70547
68669
53194
35852
29887
27041
26896
32329
46962
52020
45668
52785
52295
53960
53097
41967
33549
29440
26612
32571
46318
44779
39658
51174
53645
52989
47731
37277
31208
27607
26069
33651
43407
45722
42180
64652
88347
59052
45261
39340
37570
32789
33524
40106
55143
61645
53653
99
Среднее
Индекс
значение сезонности
59231,25
1,3113
66208,50
1,4658
58667,50
1,2988
49820,75
1,1030
38609,00
0,8548
33053,50
0,7318
29219,25
0,6469
28275,25
0,6260
34664,25
0,7674
47957,50
1,0617
51041,50
1,1300
45289,75
1,0027
Отзывы:
Авторизуйтесь, чтобы оставить отзыв