Некоторые оценки непрерывности для вероятности разорения и других связанных с разорением величин

Основная информация

• ID статьи: 2511.12218
• Название: Some continuity estimates for ruin probability and other ruin-related quantities
• Автор: Lazaros Kanellopoulos (Университет Эгейского моря, Греция)
• Классификация: math.PR (Теория вероятностей)
• Дата публикации: 18 ноября 2025 г. (препринт arXiv)
• Ссылка на статью: https://arxiv.org/abs/2511.12218

Аннотация

В данной работе исследуются свойства непрерывности вероятности разорения в классической модели риска. Используя свойства сжимающих интегральных операторов, выводятся оценки непрерывности для дефицита в момент разорения. Эти результаты также применяются к непрерывным процессам прибыли с диффузионными возмущениями, в рамках которых вероятность разорения может быть представлена как свёртка составного геометрического распределения K с диффузионным членом. В статье выводятся неравенства непрерывности для K и предлагается итерационный метод приближения для связанных с разорением величин, результаты проверяются на числовых примерах.

Исследовательский контекст и мотивация

Основные проблемы

1. Проблема устойчивости: В актуарной математике фактические параметры риска (такие как интенсивность поступления исков λ и распределение размеров исков F) обычно неизвестны и должны оцениваться приблизительными параметрами. Как измерить влияние возмущений параметров на выходные величины, такие как вероятность разорения?
2. Оценки непрерывности: Поиск неравенств вида ν(ψ_α, ψ_α̃) ≤ w(δ(α, α̃)), где δ измеряет расстояние между входными параметрами, а ν измеряет расстояние между выходными функциями.
3. Модель с диффузионными возмущениями: Для моделей риска, содержащих диффузионный член, в литературе отсутствуют работы, использующие вероятностные метрики для анализа непрерывности.

Значимость исследования

• Теоретическое значение: Совершенствование рамок анализа устойчивости в теории риска, обеспечение теоретических гарантий для распространения ошибок оценки параметров
• Практическая ценность: Помощь актуариям в оценке границ ошибок при расчёте вероятности разорения с использованием приблизительных параметров
• Методологическое новшество: Применение принципа Банаха о сжимающем отображении и теории неподвижных точек к новым областям теории разорения

Ограничения существующих методов

• Методы сильной устойчивости в основном применяются к классическим моделям без диффузии
• Явные решения существуют только для специальных распределений (экспоненциальное, фазового типа)
• Недостаточно инструментов для анализа непрерывности в моделях с диффузионными возмущениями

Основные вклады

1. Неравенства непрерывности для вероятности разорения (теорема 1): Для распределений исков с конечным (γ+1)-м моментом установлены границы непрерывности вероятности разорения на основе метрики ν_γ
2. Границы непрерывности для дефицита в момент разорения (теорема 2): Выведены границы единообразной метрики между двумя распределениями дефицита в момент разорения вида: $\sup_{u≥0}|G(u,y) - \tilde{G}(u,y)| ≤ \frac{1}{c-λμ}[λQ_y(F,\tilde{F}) + |λ-\tilde{λ}|\tilde{μ}]$
3. Непрерывность K(u) в модели с диффузией (теорема 3): Впервые установлены оценки непрерывности для составного геометрического распределения K(u) в модели с диффузионными возмущениями
4. Итерационный метод приближения (лемма 1): На основе теоремы Банаха о неподвижной точке предложена явная итерационная формула приближения K(u) с чёткими гарантиями сходимости
5. Численная проверка: На нескольких примерах проверена эффективность теоретических границ (хотя границы относительно широкие)

Подробное описание методов

Постановка задачи

Даны два процесса риска:

• Процесс 1: параметры (λ, c, F), вероятность разорения ψ(u), дефицит в момент разорения G(u,y)
• Процесс 2: параметры (λ̃, c, F̃), вероятность разорения ψ̃(u), дефицит в момент разорения G̃(u,y)

Цель: Установить неравенства вида d_out(ψ, ψ̃) ≤ f(d_in(λ,F; λ̃,F̃))

где d_in и d_out — надлежащим образом выбранные метрики, f — возрастающая функция, удовлетворяющая f(0)=0.

Теоретическая база

1. Конструкция банахова пространства

• Пространство D_γ: Функции h, удовлетворяющие h_γ(x) = (1+x)^γ|h(x)| ∈ D_∞ (функции с правой непрерывностью и левыми пределами)
• Метрика: ν_γ(x,y) = ∫_0^∞ (1+t)^γ|x(t)-y(t)|dt
• Полнота: (D_γ, ν_γ) — полное метрическое пространство

2. Сжимающий оператор

Определим интегральный оператор: $Tx(u) = \frac{λ}{c}\left[\int_u^∞ F(t)dt + \int_0^u x(u-t)F(t)dt\right]$

Ключевые свойства:

• T: D_γ → D_γ
• Коэффициент сжатия: λM_γ^X/c < 1 (при условии положительной чистой прибыли)
• где M_γ^X = ∫_0^∞(1+t)^γF(t)dt = E(X+1)^(γ+1)-1/(γ+1)

3. Метод неподвижной точки

• Вероятность разорения ψ является единственной неподвижной точкой оператора T: ψ = Tψ
• Использование неравенства треугольника: ν_γ(ψ,ψ̃) ≤ ν_γ(Tψ,Tψ̃) + ν_γ(Tψ̃,T̃ψ̃)
• Первый член контролируется сжимаемостью, второй измеряет разницу операторов

Подробное описание теоремы 1 (непрерывность вероятности разорения)

Формулировка: Предположим EX^(γ+1) < ∞, EX̃^(γ+1) < ∞, тогда: $\int_0^∞(1+t)^γ|ψ(t)-\tilde{ψ}(t)|dt ≤ \frac{c}{c-λM_γ^X}\left[\frac{ν_{γ+1}(F,\tilde{F})}{γ+1} + ν_γ(F,\tilde{F})M_γ^L + \frac{|λ-\tilde{λ}|}{c}M_{γ+1}^{\tilde{X}}(1+M_γ^L)\right]$

Основные этапы доказательства:

1. Установление ν_γ(Tx,Ty) ≤ (λM_γ^X/c)ν_γ(x,y) (сжимаемость)
2. Оценка ν_γ(Tψ̃,T̃ψ̃), разложение на I_1 (разница распределений) + I_2 (разница интенсивностей)
3. Для I_1 использование интегрирования по частям:
   • Члены ∫0^t(1+u)^γdu дают ν{γ+1}(F,F̃)/(γ+1)
   • Свёрточные члены используют неравенство (1+z+t)^γ ≤ (1+z)^γ(1+t)^γ
4. Аналогичная обработка I_2

Частный случай (γ=0): Сводится к границе в метрике Канторовича

Подробное описание теоремы 2 (дефицит в момент разорения)

Предпосылка: Хвост дефицита в момент разорения G(u,y) = ψ(u) - G(u,y) удовлетворяет: $G(u,y) = \frac{λ}{c}\left[\int_0^u G(u-t,y)F(t)dt + \int_{u+y}^∞ F(t)dt\right]$

Формулировка: $\sup_{u≥0}|G(u,y)-\tilde{G}(u,y)| ≤ \frac{1}{c-λμ}[λQ_y(F,\tilde{F}) + |λ-\tilde{λ}|\tilde{μ}]$

где Q_y(F,F̃) = ∫_y^∞|F(t)-F̃(t)|dt

Стратегия доказательства:

• Использование единообразной метрики ν_d(x,y) = sup_t|x(t)-y(t)|
• Фиксирование y, рассмотрение G(·,y) как функции u
• Проверка сжимаемости оператора T^z (коэффициент λμ/c)
• Оценка ν_d(T^zΓ̃^z, T̃^zΓ̃^z)

Подробное описание теоремы 3 (K(u) в модели с диффузией)

Модель: U(t) = u + ct - S(t) + σB(t)

Ключевая величина: K(u) = Pr(L_K > u), где L_K — совокупная потеря от рекордных максимумов, вызванных исками

Структура: K(u) удовлетворяет несобственному интегральному уравнению: $K(u) = \frac{λμ}{c}\left[A(u) + \int_0^u K(u-t)dA(t)\right]$

где A(x) — функция распределения L_o + L_c, L_o ~ H_1 (диффузионный член), L_c ~ F_e (сбалансированное распределение)

Граница непрерывности: $\sup_u|K(u)-\tilde{K}(u)| ≤ \frac{1}{c-λμ}\left[λμ\left(\frac{c}{D}K(H_1,\tilde{H}_1) + \frac{|\tilde{D}-D|}{D} + \frac{K(F,\tilde{F})}{μ} + \frac{|\tilde{μ}-μ|}{μ}\right) + |λμ-\tilde{λ}\tilde{μ}|\right]$

Технические сложности: Требуется оценка ∫|a(t)-ã(t)|dt, где a — плотность свёртки H_1 и H_2

Итерационное приближение (лемма 1)

Итерационная схема: K_n = T_dK_, K_0 ∈ K_d

Явная формула (K_0=k — константа):

• n=1: K_1(u) = φ - (1-k)φA(u)
• n≥2: K_n(u) = φ - (1-k)φ^nA^{*n}(u) - (1-φ)∑_^{n-1}φ^iA^{*i}(u)

где φ = 1/(1+θ)

Сходимость: Гарантируется теоремой Банаха о неподвижной точке ν_d(K_n,K) → 0

Частный случай (экспоненциальные иски, β=c/D): $K_n(u) = φe^{-βu}S_1(βu) + \sum_{m=2}^{n-1}φ^me^{-βu}[S_{2m-1}(βu)-S_{2m-3}(βu)] + ...$

где S_m(z) = ∑_^m z^r/r!

Экспериментальная установка

Конструкция численных примеров

Пример 1 (проверка теоремы 1):

• Распределения: X ~ смешанное экспоненциальное (1/2·Exp(5/4) + 1/2·Exp(5/6)), X̃ ~ Exp(1)
• Параметры: EX = EX̃ = 1, множество комбинаций (λ,c)
• Метрики: γ = 0, 1

Пример 2 (проверка теоремы 2):

• Три распределения: X^(1) ~ Erlang(3,3), X^(2) ~ Exp(1), X^(3) ~ смешанное экспоненциальное
• Параметры: θ = 1 или 4, множество пар (u,y)
• Сравнение: |G^(i)(u,y) - G^(j)(u,y)| с теоретической границей DK2

Пример 3 (проверка теоремы 3):

• Распределения: X ~ Exp(3), X̃ ~ смешанное экспоненциальное
• Параметры диффузии: множество комбинаций (D,D̃)
• Сравнение: sup|K(u)-K̃(u)| с границей DK3

Пример 4 (итерационное приближение):

• Распределение: X ~ Exp(β)
• Начальные значения: k = 0.0, 0.1, ..., 1.0
• Количество итераций: n = 1,...,5
• Два набора параметров

Методы вычисления

• Точные значения: использование известных явных решений (экспоненциальный, Erlang случаи)
• Численное интегрирование: расчёт метрик и границ
• Итерационные вычисления: использование явной формулы из леммы 1

Результаты экспериментов

Результаты примера 1 (таблица 1)

Наблюдения:

• γ=0, λ=5/6:
  • Истинное расстояние: 0.0154 (c=3) → 0.0060 (c=7)
  • Теоретическая граница DK1: 0.1211 → 0.0929
  • Степень консервативности границы: примерно 8-15 раз
• При γ=1 границы ещё более консервативны (примерно 8-13 раз)
• При увеличении λ граница слегка возрастает
• При увеличении c истинное расстояние и граница уменьшаются

Вывод: Границы консервативны, но качественные тенденции верны

Результаты примера 2 (таблица 2)

Четыре группы сравнений:

1. θ=1, X^(1) vs X^(2): истинное расстояние 0.0080-0.0832, граница 0.1547-0.2900
2. θ=4, X^(1) vs X^(2): истинное расстояние 0.0034-0.0309, граница 0.0387-0.0735
3. θ=4, X^(1) vs X^(3): истинное расстояние 0.0035-0.0328, граница 0.0413-0.0779
4. θ=4, X^(3) vs X^(2): истинное расстояние 0.0001-0.0023, граница 0.0027-0.0054

Обнаруженные закономерности:

• При увеличении θ (увеличение коэффициента безопасности) расстояние и граница значительно уменьшаются
• При увеличении y расстояние сначала растёт, затем уменьшается (существует максимум)
• Для похожих распределений (X^(3) и X^(2)) расстояние очень мало (< 0.003)
• Степень консервативности границы: примерно 2-10 раз, более узкая при θ=4

Замечание авторов: Консервативность границ объясняется требованием единообразности по всем u

Результаты примера 3 (таблица 3)

Влияние параметров диффузии:

• D=1, D̃=1/10: истинное расстояние 0.0854, граница 0.4837 (5.66 раз)
• D=3, D̃=1/20: истинное расстояние 0.1334, граница 0.5254 (3.94 раз)

Тенденции:

• Чем больше различие между D и D̃, тем больше расстояние
• При увеличении D истинное расстояние растёт
• Степень консервативности границы: примерно 3-7 раз

Анализ соотношений: Столбец |1-Ratio| показывает относительную ошибку 2.94-6.75

Результаты примера 4 (таблицы 4-5)

Скорость сходимости (таблица 4, β=2, точное значение K(1)=0.3325717):

• n=1: 0.2030-0.5000 (зависит от k)
• n=3: 0.3316-0.3336 (ошибка <0.001)
• n=5: примерно 0.3326 (ошибка <10^-5)

Чувствительность к начальному значению:

• k=0.0, n=5: ошибка 0.0000008
• k=1.0, n=5: ошибка 0.0000007
• Вывод: После 5 итераций результат не зависит от начального значения

Таблица 5 (β=3/2, точное значение K(1)=0.6573777):

• Аналогичная быстрая сходимость
• При n=5 ошибка для всех начальных значений <3×10^-6

Практическая ценность: 3-5 итераций достаточно для получения высокоточного приближения

Связанные работы

Методы устойчивости

• Beirlant & Rachev (1987): Первая формализация проблемы устойчивости в актуарной математике
• Kalashnikov (2000): Методы сильной устойчивости для вероятности разорения
• Enikeeva et al. (2001): Оценки непрерывности вероятности разорения
• Gordienko & Vázquez-Ortega (2016): Простые оценки непрерывности с использованием сжимающих отображений (прямое расширение в данной работе)

Многомерные и марковские модели с переключением режимов

• Benouaret & Aı̈ssani (2010): Сильная устойчивость двумерной классической модели риска
• Rusaityte (2001): Устойчивость модели с марковским переключением и инвестициями
• Gordienko & Vázquez-Ortega (2018): Оценки непрерывности для многомерной модели риска обновления

Применение принципа сжатия

• Gajek & Rudź (2018): Принцип Банаха о сжимающем отображении для моделей с переключением режимов
• Jiang (2021): q-масштабирующие функции и окончательная вероятность разорения
• Sánchez & Baltazar-Larios (2020): Приближение вероятности разорения с использованием теоремы о неподвижной точке
• Gajek & Rudź (2025): Анализ неподвижных точек для проблем разорения

Модели с диффузией

• Gerber (1970): Введение классической модели с диффузионными возмущениями
• Dufresne & Gerber (1991): Систематическое исследование трёх классов вероятностей разорения
• Tsai (2003, 2006): Структура и свойства K(u)
• Veraverbeke (1993): Анализ асимптотического поведения

Новшество данной работы: Первое систематическое применение методов вероятностных метрик к моделям с диффузионными возмущениями

Заключение и обсуждение

Основные выводы

1. Теоретические вклады:
   • Установлены три основных неравенства непрерывности (теоремы 1-3)
   • Предоставлены вычислимые границы ошибок
   • Расширены методы сжимающих операторов на модели с диффузией
2. Методологические вклады:
   • Единая рамка: выбор подходящего банахова пространства и метрики
   • Итерационное приближение: явные формулы + гарантии сходимости
   • Численная проверка: несколько примеров поддерживают теоретические результаты
3. Практическое значение:
   • Количественная оценка влияния ошибок оценки параметров
   • Контроль ошибок приближённых методов
   • Быстрое вычисление K(u) в моделях с диффузией

Ограничения

1. Консервативность границ:
   • Численные примеры показывают, что границы в 3-15 раз больше истинного расстояния
   • Причина: требование единообразности по всем u и многократное применение неравенства треугольника
   • Пространство для улучшения: более тонкий анализ или адаптивные границы
2. Условия на моменты:
   • Теорема 1 требует EX^(γ+1) < ∞
   • Исключает некоторые распределения с тяжёлыми хвостами (например, Парето с α≤γ+1)
   • Для тяжелохвостых случаев требуются другие методы
3. Ограничения численных примеров:
   • Главным образом используются распределения экспоненциального семейства (имеют явные решения)
   • Проверка для более общих распределений ограничена
   • Многомерные случаи не рассмотрены
4. Предположения в модели с диффузией:
   • Теорема 3 предполагает D≥D̃ и μ≥μ̃
   • Необходимость этих предположений недостаточно обсуждена
   • Методы для случаев, когда они не выполняются, не предложены
5. Вычислительная сложность:
   • Вычисление границ включает множество интегралов
   • Практическая осуществимость в реальных приложениях не подробно обсуждена

Направления будущих исследований

1. Теоретические улучшения:
   • Поиск более узких границ
   • Обработка распределений с тяжёлыми хвостами
   • Расширение на многомерные модели
2. Расширение методов:
   • Другие типы моделей риска (например, процессы Леви)
   • Параметры, зависящие от времени
   • Зависимые иски
3. Прикладные исследования:
   • Оценка параметров на реальных данных
   • Анализ робастности
   • Оптимальные стратегии приближения

Глубокая оценка

Достоинства

1. Математическая строгость:
   • Полные детали доказательств
   • Ясная математическая рамка
   • Явные предположения и условия
2. Методологическое новшество:
   • Первое систематическое применение методов вероятностных метрик к моделям с диффузией
   • Единая обработка вероятности разорения и дефицита в момент разорения
   • Явные формулы для итерационного приближения (лемма 1)
3. Ясная структура:
   • Переход от простого к сложному (классический → дефицит → диффузия)
   • Каждая теорема сопровождается численной проверкой
   • Хороший баланс между техническими деталями и интуитивным объяснением
4. Практическая ценность:
   • Предоставлены вычислимые границы ошибок
   • Итерационный метод быстро сходится (3-5 итераций)
   • Применим к анализу чувствительности параметров
5. Обзор литературы:
   • Полное охватывание связанных работ
   • Ясное позиционирование вклада данной работы
   • Надлежащие ссылки (32 источника)

Недостатки

1. Эффективность границ:
   • Численные примеры показывают консервативные границы (3-15 раз)
   • Может ограничить практическую применимость
   • Отсутствует анализ тесноты границ (можно ли улучшить?)
2. Численные эксперименты:
   • Главным образом ограничены экспоненциальным семейством (легко вычислить точные решения)
   • Отсутствует проверка методом Монте-Карло для общих распределений
   • Вычислительные затраты не обсуждены
3. Теоретические ограничения:
   • Условия на моменты исключают тяжелохвостые случаи
   • Односторонние предположения в теореме 3 (D≥D̃) не обсуждены
   • Недостаточно рекомендаций по выбору оптимальной метрики
4. Рекомендации по применению:
   • Отсутствуют советы по выбору γ
   • Интеграция с практической оценкой параметров не обсуждена
   • Сравнение с другими методами (например, сильной устойчивостью) недостаточно
5. Детали изложения:
   • Некоторые обозначения недостаточно объяснены (например, смысл M^L_γ)
   • Замечание в теореме 2 упоминает "единообразность", но не развивает
   • Мотивация выбора параметров в примере 4 неясна

Оценка влияния

Теоретическое влияние:

• Предоставляет новые инструменты анализа для теории разорения
• Расширяет область применения методов сжимающих операторов
• Может стимулировать исследования более узких границ

Практическое влияние:

• Среднее: консервативность границ ограничивает прямое применение
• Итерационный метод имеет практическую ценность (быстрая сходимость)
• Требуется дополнительная работа для использования в управлении рисками

Воспроизводимость:

• Высокая: параметры численных примеров полные
• Теоретические доказательства подробны
• Отсутствует код, но формул достаточно для реализации

Вклад в область:

• Заполняет пробел в анализе непрерывности для моделей с диффузией
• Предоставляет систематическую методологическую рамку
• Создаёт основу для последующих исследований

Применимые сценарии

Подходящие сценарии:

1. Анализ чувствительности параметров: Оценка влияния ошибок оценки параметров
2. Проверка приближённых методов: Предоставление границ ошибок для численных методов
3. Теоретические исследования: Использование как эталона для более сложных моделей
4. Обучение: Демонстрация применения методов сжимающих операторов

Неподходящие сценарии:

1. Точный контроль ошибок: Границы могут быть чрезмерно консервативны
2. Распределения с тяжёлыми хвостами: Условия на моменты не выполняются
3. Высокомерные сложные модели: Расширение методов нетривиально
4. Вычисления в реальном времени: Вычисление интегралов может быть медленным

Рекомендации по улучшению:

• Комбинирование с методом Монте-Карло для проверки более общих случаев
• Разработка адаптивных границ (в зависимости от конкретных параметров)
• Создание программного пакета для реализации
• Интеграция с методами машинного обучения для оценки параметров

Избранные источники

Основополагающие работы:

• Asmussen & Albrecher (2010): Стандартный учебник по теории вероятности разорения
• Dufresne & Gerber (1991): Основополагающая работа по моделям с диффузионными возмущениями

Методология:

• Gordienko & Vázquez-Ortega (2016): Работа, непосредственно расширяемая в данной статье
• Gajek & Rudź (2018, 2025): Применение принципа Банаха о сжимающем отображении

Модели с диффузией:

• Tsai (2003, 2006): Структура и свойства K(u)
• Veraverbeke (1993): Анализ асимптотического поведения

Общая оценка: Это математически строгая и методологически инновационная теоретическая работа, предоставляющая новые инструменты анализа для теории разорения. Основной вклад заключается в систематическом применении методов вероятностных метрик к моделям с диффузионными возмущениями и предоставлении вычислимых границ непрерывности. Хотя численные границы относительно консервативны, что ограничивает прямую практическую применимость, вклад в теорию и методологию значителен. Статья подходит для публикации в специализированных журналах по актуарной математике или прикладной теории вероятностей и имеет важное справочное значение для исследователей в этой области. Будущие работы должны сосредоточиться на улучшении тесноты границ и расширении на более общие модели.