Analyzing distortion riskmetrics and weighted entropy for unimodal and symmetric distributions under partial information constraints

基本信息

• 论文ID: 2504.19725
• 标题: Analyzing distortion riskmetrics and weighted entropy for unimodal and symmetric distributions under partial information constraints
• 作者: Baishuai Zuo (Qufu Normal University & Southern University of Science and Technology), Chuancun Yin (Qufu Normal University)
• 分类: q-fin.RM (Quantitative Finance - Risk Management), stat.AP (Statistics - Applications)
• 发表时间: 2025年11月21日 (arXiv v2)
• 论文链接: https://arxiv.org/abs/2504.19725v2

摘要

本文在已知均值和方差信息的条件下，针对单峰分布、对称分布以及对称单峰分布，推导了最坏情况下扭曲风险度量(distortion riskmetrics)和加权熵(weighted entropy)的上下界。特别地，对于对称分布，本文给出了精确的上界。这些理论结果被应用于多种(加权)熵、短缺量(shortfalls)及其他风险度量，包括累积Tsallis过去熵、α阶累积残差Tsallis熵、扩展Gini系数、分数广义累积残差熵、扩展Gini短缺、n阶期望短缺等。

研究背景与动机

研究问题

本文要解决的核心问题是：在部分信息约束下（仅知道均值和方差），如何确定具有特定分布特征（单峰性、对称性）的随机变量的扭曲风险度量和加权熵的最坏情况界。

问题重要性

1. 实际需求：在现实中，许多数据集是无模型的或存在分布不确定性，但我们通常能获得其部分信息（如均值、方差）和分布的结构特征（如单峰性、对称性）
2. 风险管理应用：扭曲风险度量在行为经济学和风险管理中有广泛应用，特别是在保费原则构建和风险厌恶建模方面
3. 理论价值：熵在生命表研究、年金成本计算和保费原则构建中扮演重要角色

现有方法的局限性

1. 一般分布研究：Zuo and Yin (2025)研究了一般分布下的最坏情况扭曲风险度量和加权熵，但未考虑分布的特殊结构
2. 缺少对称约束：Psarrakos et al. (2024)和Zuo and Yin (2025)提供了部分信息下的最坏情况结果，但都没有考虑对称约束
3. 单峰分布研究不足：虽然Bernard et al. (2020, 2025)研究了单峰分布的RVaR界，但针对一般扭曲风险度量和加权熵的单峰、对称单峰分布研究尚未完成

研究动机

1. 扭曲风险度量和加权熵具有广泛应用价值
2. 实践中"已知部分信息和分布特殊结构来快速估计风险度量"是常见需求
3. 填补特殊分布结构下最坏情况风险度量研究的空白

核心贡献

1. 单峰分布的界：推导了单峰分布下扭曲风险度量和加权熵的下界和上界（定理1、推论1、定理4）
2. 对称分布的精确界：给出了对称分布下扭曲风险度量和加权熵的精确上界及最坏情况分布（定理2、推论2、命题1、定理5）
3. 对称单峰分布的界：建立了对称单峰分布下扭曲风险度量和加权熵的下界和上界（定理3、推论3、定理6）
4. 广泛应用：将理论结果应用于多种具体的熵度量、短缺量和风险度量，包括：
   • 熵：累积Tsallis过去熵、累积残差Tsallis熵、Gini均值半差、分数广义累积残差熵、累积残差熵等
   • 短缺量：扩展Gini短缺、Gini短缺、累积残差熵短缺等
   • 其他风险度量：n阶期望短缺、对偶幂原则、比例风险原则等
5. 数值验证：通过真实股票收益数据验证了理论结果的实用性

方法详解

任务定义

问题1：对于扭曲风险度量 $\rho_g(X) = \int_0^\infty g(P(X > x))dx + \int_{-\infty}^0 [g(P(X > x)) - g(1)]dx$

求解优化问题： $\sup_{X \in \mathcal{V}(\mu, \sigma)} \rho_g(X)$

其中$\mathcal{V}(\mu, \sigma) \in \{\mathcal{V}_U(\mu, \sigma), \mathcal{V}_S(\mu, \sigma), \mathcal{V}_{SU}(\mu, \sigma)\}$分别表示具有均值$\mu$和方差$\sigma^2$的单峰、对称、对称单峰随机变量集合。

问题2：对于加权熵 $\int_{-\infty}^{+\infty} \psi(x)g(\bar{F}_X(x))dx$

在加权约束$E[\Psi(F_X^{-1}(U))] = \mu_\Psi$, $\text{Var}[\Psi(F_X^{-1}(U))] = \sigma_\Psi^2$下，求解类似的优化问题。

核心方法框架

1. 分位数表示（Quantile Representation）

利用引理1和引理2，将扭曲风险度量和加权熵转换为分位数函数的积分形式： $\rho_g(X) = \int_0^1 F_X^{-1}(u)d\hat{g}(u)$$\int_{-\infty}^{+\infty} \psi(x)g(\bar{F}_X(x))dx = \int_0^1 \Psi(F_X^{-1}(u))d\hat{g}(u)$

其中$\hat{g}(u) = g(1) - g(1-u)$。

2. 极端分布构造

单峰分布：定义两类极端分布集合$U_R$和$U_L$（或加权版本$U_{\Psi R}$和$U_{\Psi L}$）： $U_R = \left\{X: F_X^{-1}(u) = \begin{cases} a, & u \in [0, b) \\ c(u-b) + a, & u \in [b, 1] \end{cases}\right\}$

这些分布在给定约束下提供了下界。

对称单峰分布：定义集合$\mathcal{S}(\alpha)$： $\mathcal{S}(\alpha) = \left\{X: F_X^{-1}(u) = \begin{cases} a + c(u-1+b), & u \in (0, 1-b) \\ a, & u \in [1-b, b] \\ a + c(u-b), & u \in (b, 1) \end{cases}\right\}$

3. 凸包络技术（Convex Envelope）

对于上界，利用凸包络$\hat{g}_*$（$\hat{g}$的最大凸下界）和引理3（单调性引理）： $\rho_g(X) = \int_0^1 F_X^{-1}(u)d\hat{g}(u) \leq \int_0^1 F_X^{-1}(u)d\hat{g}_*(u)$

4. TVaR连接

通过将问题转换为尾部条件期望（TVaR）的形式： $\sup_{X \in \mathcal{V}_U(\mu, \sigma)} \rho_g(X) \leq \mu(\hat{g}_*)^\prime(0) + \int_0^1 (1-u)\sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{TVaR}_u(X) d(\hat{g}_*)^\prime(u)$

利用Bernard et al. (2020)关于单峰分布TVaR的已知结果： $\sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{TVaR}_\alpha(X) = \begin{cases} \mu + \frac{\sigma}{3}\sqrt{\frac{\alpha(8-9\alpha)}{1-\alpha}}, & \alpha \in (0, \frac{1}{2}) \\ \mu + \sigma\sqrt{\frac{8}{9(1-\alpha)} - 1}, & \alpha \in [\frac{1}{2}, 1) \end{cases}$

技术创新点

1. 统一框架：将扭曲风险度量和加权熵纳入统一的数学框架，通过分位数表示和凸包络技术系统地处理
2. 精确对称界：对于对称分布，通过利用对称性约束$F_X^{-1}(u) + F_X^{-1}(1-u) = 2\mu$，得到精确的闭式解
3. 分段线性逼近：使用分段线性函数序列逼近凸包络，结合单调收敛定理建立上界
4. 残差寿命和过去寿命扩展：将结果扩展到残差寿命$X_t = [X-t|X>t]$和过去寿命$X_{(t)} = [X|X \leq t]$的情况（推论4-9）

实验设置

数据集

作者使用了三家公司的日股票收益数据：

1. Coca-Cola Company (KO)：$\mu_1 = 0.052128514$, $\sigma_1^2 = 0.416460169$
2. Cincinnati Financial (CINF)：$\mu_2 = -0.009317269$, $\sigma_2^2 = 3.604281185$
3. eBay Inc. (EBAY)：$\mu_3 = 0.09911245$, $\sigma_3^2 = 0.900906979$

数据期间：2024-05-20至2024-05-18的工作日数据（来自nasdaq.com）

评价指标

计算Gini短缺（$\text{GS}_p^\tau$）的精确上界，参数设置：

• $\tau = 0.25$（加载参数）
• $r = 2$
• $p \in \{0.80, 0.90, 0.95, 0.98, 0.99\}$（置信水平）

对比方法

比较四种分布约束下的上界：

1. 一般分布（$\mathcal{V}(\mu, \sigma)$）
2. 单峰分布（$\mathcal{V}_U(\mu, \sigma)$）
3. 对称分布（$\mathcal{V}_S(\mu, \sigma)$）
4. 对称单峰分布（$\mathcal{V}_{SU}(\mu, \sigma)$）

实现细节

• 通过直方图分析验证数据的单峰性和对称性
• 使用数值积分和优化方法计算理论界
• 对于涉及不完全Gamma函数的情况，使用标准数学函数库

实验结果

主要结果

表1-3：三只股票的Gini短缺上界

以KO股票为例（表1）：

关键发现

1. 单调性：所有上界随置信水平$p$单调递增，符合风险度量的直觉
2. 约束效应的排序：对于相同的$p$值， $\sup_{\mathcal{V}} \text{GS}_p^\tau > \sup_{\mathcal{V}_S} \text{GS}_p^\tau > \sup_{\mathcal{V}_U} \text{GS}_p^\tau > \sup_{\mathcal{V}_{SU}} \text{GS}_p^\tau$
   这表明单峰性约束比对称性约束对上界的影响更大
3. 方差的影响：CINF（$\sigma_2^2 = 3.604$）的上界最大，KO（$\sigma_1^2 = 0.416$）的上界最小，反映了方差对风险度量的直接影响
4. 分布形状不变性：对于不同的股票，约束的相对效应保持一致，说明理论结果的稳健性

具体应用案例

例1：累积Tsallis过去熵（$\text{CT}_\alpha(X)$）

对于单峰分布（$\alpha = 2$，即Gini均值半差）： $\sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{Gini}(X) \geq \sigma\sqrt{3}$

对于对称分布： $\sup_{X \in \mathcal{V}_S(\mu, \sigma)} \text{Gini}(X) = \sigma\sqrt{3}$ 最坏情况分布为：$F_X^{-1}(u) = \mu + \sigma\sqrt{3}(2u-1)$

例2：累积残差熵（$\alpha = 1$）

对于单峰分布： $\sup_{X \in \mathcal{V}_U(\mu, \sigma)} E(X) \geq 0.927341\sigma$

对于对称分布： $\sup_{X \in \mathcal{V}_S(\mu, \sigma)} E(X) = \frac{\pi\sigma}{2\sqrt{3}}$

对于对称单峰分布： $\sup_{X \in \mathcal{V}_{SU}(\mu, \sigma)} E(X) \geq 0.878282\sigma$

例3：期望短缺（$\text{ES}_p(X)$）

对于单峰分布： $\sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{ES}_p(X) = \begin{cases} \mu + \frac{\sqrt{p(8-9p)}}{3(1-p)}\sigma, & p \in (0, \frac{1}{2}] \\ \mu + \frac{1}{3}\sqrt{\frac{9p-1}{1-p}}\sigma, & p \in (\frac{1}{2}, 1) \end{cases}$

这一结果推广了Li et al. (2018)的引理2.6（仅要求$p \in (\frac{5}{6}, 1)$）。

理论验证

备注5：通过不等式链 $\sup_{X \in \mathcal{V}_{SU}(\mu, \sigma)} \text{Gini}(X) \leq \sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{Gini}(X) \leq \sup_{X \in \mathcal{V}(\mu, \sigma)} \text{Gini}(X)$

结合Zuo and Yin (2025)的结果$\sup_{X \in \mathcal{V}(\mu, \sigma)} \text{Gini}(X) = \sigma\sqrt{3}$和本文的下界$\sup_{X \in \mathcal{V}_{SU}(\mu, \sigma)} \text{Gini}(X) \geq \sigma\sqrt{3}$，得出： $\sup_{X \in \mathcal{V}_{SU}(\mu, \sigma)} \text{Gini}(X) = \sup_{X \in \mathcal{V}_U(\mu, \sigma)} \text{Gini}(X) = \sigma\sqrt{3}$

这说明对于Gini系数，单峰约束和对称单峰约束给出相同的精确界。

相关工作

部分信息下的风险度量

1. VaR和TVaR：El Ghaoui et al. (2003)和Chen et al. (2011)分别给出了已知前两阶矩时VaR和TVaR的闭式解
2. RVaR：Li et al. (2018)研究了RVaR的最坏情况；Bernard et al. (2020, 2025)进一步研究了单峰和对称单峰分布的RVaR界
3. 谱风险度量：Li (2018)研究了谱风险度量（凸扭曲风险度量的重要子类）
4. 对称信息：Zhu and Shao (2018)在对称信息下推广了Li (2018)的结果

扭曲风险度量

1. 理论基础：Wang et al. (2020a,b)建立了扭曲风险度量的一般理论框架
2. 高阶矩：Liu et al. (2020)和Cai et al. (2025a)研究了已知高阶矩时的最坏情况
3. Wasserstein球：Bernard et al. (2024)和Cai et al. (2025b)研究了分布位于Wasserstein球中的情况
4. 极端情况：Shao and Zhang (2023, 2024)和Zhao et al. (2024)给出了已知前两阶矩和对称性时的闭式解

本文的定位

本文是首次系统研究单峰和对称单峰约束下扭曲风险度量和加权熵的工作，填补了：

• Zuo and Yin (2025)未考虑分布特殊结构的空白
• Bernard et al. (2020, 2025)仅针对RVaR的局限
• Zhao et al. (2024)未给出单峰和对称单峰分布精确界的不足

结论与讨论

主要结论

1. 单峰分布：建立了扭曲风险度量和加权熵的可计算上下界，下界由极端分布$U_R$和$U_L$给出，上界通过凸包络和TVaR界得到
2. 对称分布：得到了精确的闭式上界 $\sup_{X \in \mathcal{V}_S(\mu, \sigma)} \rho_g(X) = \mu g(1) + \frac{\sigma}{2}\sqrt{\int_0^1 [(\hat{g}(u) + \hat{g}(1-u))^\prime_*]^2 du}$
3. 对称单峰分布：建立了上下界，其中下界由集合$\mathcal{S}(\alpha)$给出，上界通过凸包络技术得到
4. 广泛适用性：理论结果成功应用于15种以上的具体风险度量和熵度量

局限性

1. 计算复杂性：
   • 单峰和对称单峰分布的界不是闭式的，需要数值优化
   • 涉及复杂的积分和特殊函数（如不完全Gamma函数）
2. 界的紧性：
   • 对于单峰和对称单峰分布，上下界之间可能存在间隙
   • 某些情况下（如Gini系数）界是紧的，但一般情况下紧性未知
3. 分布假设：
   • 仅考虑单峰分布（单个众数）
   • 未涉及多峰或其他更复杂的分布结构
4. 矩约束：
   • 仅使用前两阶矩信息
   • 更高阶矩可能提供更紧的界

未来方向

作者明确提出了两个研究方向：

1. 优化问题：在部分信息下，研究单峰和对称单峰分布的扭曲风险度量最优再保险问题
2. 多峰分布：研究多峰分布下扭曲风险度量和加权熵的精确上界
3. 其他扭曲度量：将结果扩展到其他扭曲度量族，如参数化Wang保费（Fiori and Rosazza Gianin, 2023）

深度评价

优点

1. 理论严谨性：
   • 证明完整且严密，充分利用了凸分析、分位数函数理论和单调收敛定理
   • 对于对称分布给出了精确解，而非仅是界
2. 统一框架：
   • 将扭曲风险度量和加权熵纳入统一框架
   • 系统地处理单峰、对称、对称单峰三类分布
3. 广泛适用性：
   • 理论结果应用于15+种具体度量
   • 涵盖熵、短缺量、风险度量三大类应用
   • 扩展到残差寿命和过去寿命
4. 方法创新：
   • 巧妙结合凸包络技术和TVaR界
   • 构造的极端分布集合具有明确的几何意义
5. 实证验证：
   • 使用真实股票数据验证理论
   • 清晰展示了不同约束的相对效应

不足

1. 可计算性问题：
   • 单峰和对称单峰分布的界需要求解复杂的优化问题
   • 例1和例2中的方程$M_\alpha(b_0) = 0$等需要数值求解
   • 缺少算法实现细节和计算复杂度分析
2. 界的紧性分析不足：
   • 除了少数特例（如Gini系数），未系统分析上下界之间的间隙
   • 未提供界的渐近行为分析
3. 实验部分薄弱：
   • 仅使用3只股票的数据，样本量较小
   • 数据期间较短（约一年）
   • 缺少对数据单峰性和对称性的定量检验（仅依赖直方图）
   • 未与其他方法（如经验估计、参数方法）进行比较
4. 写作问题：
   • 论文较长（30页），结构可以更紧凑
   • 部分证明冗长，可移至附录
   • 图表质量一般（如图1-5的分辨率和标注）
5. 理论局限：
   • 单峰性定义基于分布函数的凸凹性，较为限制
   • 未考虑众数位置已知的情况
   • 对于加权熵，约束$E[\Psi(F_X^{-1}(U))] = \mu_\Psi$的实际意义不够清晰

影响力

1. 学术贡献：
   • 填补了单峰和对称单峰分布下扭曲风险度量研究的空白
   • 为后续研究提供了理论基础和技术工具
   • 预期被风险管理、精算学和决策理论领域引用
2. 实用价值：
   • 为保险和金融机构提供了保守的风险评估工具
   • 在分布不确定但知道结构信息时特别有用
   • 可用于压力测试和情景分析
3. 可复现性：
   • 理论证明详细，可以复现
   • 缺少代码和数据公开，实验部分可复现性一般
   • 部分数值结果（如$b_0 \approx 0.582812$）未说明计算方法

适用场景

1. 风险管理：
   • 保险公司的准备金计算
   • 银行的资本充足率评估
   • 投资组合的风险限额设定
2. 精算定价：
   • 保费原则的稳健设计
   • 年金产品定价
   • 再保险合同设计
3. 监管合规：
   • 当完整分布模型不可用但知道分布特征时
   • 需要保守估计的情况
   • Solvency II等监管框架下的风险评估
4. 决策分析：
   • 分布不确定性下的稳健决策
   • 行为经济学中的风险厌恶建模
5. 不适用场景：
   • 数据充分且可以可靠估计完整分布时
   • 分布明显多峰或不满足假设时
   • 需要精确值而非界的应用

参考文献（关键文献）

1. Bernard, C., Kazzi, R., Vanduffel, S. (2020, 2025)：单峰和对称单峰分布的RVaR界，本文的重要基础
2. Zuo, B., Yin, C. (2025)：一般分布下的扭曲风险度量，本文直接推广的工作
3. Wang, Q., Wang, R., Wei, Y. (2020a)：扭曲风险度量的一般理论框架
4. Shao, H., Zhang, Z.G. (2023, 2024)：参数模糊下的扭曲风险度量，提供了技术方法
5. Zhao, M., et al. (2024)：部分信息约束下的最优和最坏情况扭曲风险度量

总体评价

这是一篇理论严谨、应用广泛的风险管理领域论文。主要优势在于：(1)建立了单峰和对称分布下扭曲风险度量的系统理论框架；(2)对对称分布给出了精确解；(3)应用范围广泛。主要不足在于：(1)单峰和对称单峰分布的界需要数值计算，可操作性有待提高；(2)实验部分相对薄弱；(3)界的紧性分析不够充分。

论文适合发表在精算、风险管理或应用概率领域的期刊（如Insurance: Mathematics and Economics, ASTIN Bulletin等）。对于从业者，本文提供了实用的保守风险评估工具；对于研究者，本文提供了丰富的技术方法和未来研究方向。

推荐指数：4/5（理论贡献显著，但实验和可操作性有提升空间）