2025-11-25T15:43:18.160640

On goodness-of-fit testing for volatility in McKean-Vlasov models

Heidari, Podolskij
This paper develops a statistical framework for goodness-of-fit testing of volatility functions in McKean-Vlasov stochastic differential equations, which describe large systems of interacting particles with distribution-dependent dynamics. While integrated volatility estimation in classical SDEs is now well established, formal model validation and goodness-of-fit testing for McKean-Vlasov systems remain largely unexplored, particularly in regimes with both large particle limits and high-frequency sampling. We propose a test statistic based on discrete observations of particle systems, analysed in a joint regime where both the number of particles and the sampling frequency increase. The estimators involved are proven to be consistent, and the test statistic is shown to satisfy a central limit theorem, converging in distribution to a centred Gaussian law.
academic

On goodness-of-fit testing for volatility in McKean-Vlasov models

基本信息

  • 论文ID: 2510.12607
  • 标题: On goodness-of-fit testing for volatility in McKean-Vlasov models
  • 作者: Akram Heidari, Mark Podolskij (University of Luxembourg)
  • 分类: stat.ME (Statistics - Methodology)
  • 发表时间: 2025年10月14日 (arXiv预印本)
  • 论文链接: https://arxiv.org/abs/2510.12607

摘要

本文为McKean-Vlasov随机微分方程中的波动率函数开发了一个统计拟合优度检验框架。McKean-Vlasov方程描述了具有分布依赖动力学的大型相互作用粒子系统。虽然经典SDE中的积分波动率估计已经得到很好的建立，但McKean-Vlasov系统的正式模型验证和拟合优度检验在很大程度上仍未被探索，特别是在大粒子极限和高频采样的联合情形下。作者提出了基于粒子系统离散观测的检验统计量，在粒子数量和采样频率同时增加的联合情形下进行分析。证明了相关估计量的一致性，并表明检验统计量满足中心极限定理，分布收敛到中心化高斯律。

研究背景与动机

问题描述

McKean-Vlasov随机微分方程是描述大型相互作用粒子系统的重要数学工具，其中每个粒子的动力学不仅依赖于其个体状态，还依赖于整个系统的统计分布。这种分布依赖性使得McKean-Vlasov模型特别适合捕捉系统性相互作用和涌现行为。

研究重要性

  1. 应用广泛性: McKean-Vlasov模型在金融、物理、工程等领域有广泛应用，包括系统性风险建模、均值场博弈、大规模相互作用系统分析等
  2. 理论缺口: 虽然经典SDE的波动率估计和检验理论已经成熟，但McKean-Vlasov系统的模型验证方法仍然缺乏
  3. 实际需求: 在实际应用中，需要严格验证波动率函数的结构假设，因为错误的规格设定会显著影响下游预测和风险度量

现有方法局限性

  1. 经典方法不适用: 现有的SDE波动率检验方法(如Dette & Podolskij, 2008)仅适用于非相互作用系统
  2. 研究不充分: 现有McKean-Vlasov文献主要关注漂移函数的参数估计，波动率检验基本未被涉及
  3. 方法论空白: 缺乏处理分布依赖性和非线性效应的统计检验框架

核心贡献

  1. 首创性框架: 提出了McKean-Vlasov模型中波动率函数拟合优度检验的首个严格统计框架
  2. 双重渐近理论: 在粒子数N→∞和采样频率增加(Δₙ→0)的联合渐近框架下建立理论
  3. 一致性证明: 证明了所涉及估计量的一致性和检验统计量的中心极限定理
  4. 实用检验程序: 构造了具有正确渐近水平且对任何固定备择假设一致的检验程序
  5. 技术创新: 克服了分布依赖性带来的非线性和路径依赖效应的技术挑战

方法详解

任务定义

考虑N个相互作用粒子系统：

dX^i_t = b(X^i_t, μₜ)dt + a(X^i_t, μₜ)dW^i_t, i = 1,...,N, t ∈ [0,T]

其中μₜ是X^i_t的分布，目标是检验波动率函数a(x,μ)是否属于给定的参数族。

模型架构

假设检验框架

原假设:

H₀: L := min_{(λ₁,...,λₐ)∈ℝᵈ} ∫₀ᵀ ∫_ℝ (a²(x,μₜ) - Σᵈₖ₌₁ λₖa²ₖ(x,μₜ))² μₜ(dx)dt = 0

备择假设: H₁: L > 0

检验统计量构造

距离测度L的闭式表达：

L = B - (Γ₁,...,Γₐ)Λ⁻¹(Γ₁,...,Γₐ)ᵀ

其中：

  • B = ∫₀ᵀ ∫_ℝ a⁴(x,μₜ)μₜ(dx)dt
  • Γₖ = ∫₀ᵀ ∫_ℝ a²ₖ(x,μₜ)a²(x,μₜ)μₜ(dx)dt
  • Λₖ,ₗ = ∫₀ᵀ ∫_ℝ a²ₖ(x,μₜ)a²ₗ(x,μₜ)μₜ(dx)dt

经验估计量

基于离散观测(X^i_{tⱼ})构造估计量：

B̂ := 1/(3NΔₙ) Σᵢ₌₁ᴺ Σⱼ₌₁ⁿ |X^i_{tⱼ₊₁} - X^i_{tⱼ}|⁴

Γ̂ₖ := 1/N Σᵢ₌₁ᴺ Σⱼ₌₁ⁿ a²ₖ(X^i_{tⱼ}, μᴺ_{tⱼ})|X^i_{tⱼ₊₁} - X^i_{tⱼ}|²

Λ̂ₖ,ₗ := Δₙ/N Σᵢ₌₁ᴺ Σⱼ₌₁ⁿ a²ₖ(X^i_{tⱼ}, μᴺ_{tⱼ})a²ₗ(X^i_{tⱼ}, μᴺ_{tⱼ})

最终检验统计量：

ŜN = B̂ - Γ̂ᵀΛ̂⁻¹Γ̂

技术创新点

  1. 线性泛函导数: 利用线性泛函导数处理分布依赖性，这是处理McKean-Vlasov方程的关键技术工具
  2. 双重渐近分析: 同时处理N→∞和Δₙ→0的渐近行为，要求NΔ²ₙ→0的平衡条件
  3. U-统计量分解: 使用Hoeffding分解技术处理经验分布与真实分布的差异
  4. 半鞅理论应用: 利用Itô公式和半鞅性质处理高频统计的误差估计

实验设置

理论验证框架

本文主要是理论性工作，通过数学证明验证方法的有效性，而非传统的数值实验。

关键假设条件

  1. Assumption 1: 初始分布的矩条件
  2. Assumption 2: 系数的Lipschitz连续性和线性增长条件
  3. Assumption 3: 波动率函数的线性泛函导数存在性和光滑性

渐近条件

  • 粒子数N→∞
  • 采样间隔Δₙ→0
  • 平衡条件：NΔ²ₙ→0

实验结果

主要理论结果

定理4.1 (一致性)

在假设条件下，有：

√N(Λ̂ - Λ) = √NMΛ + oP(1)

定理4.2 (随机展开)

√N(Γ̂ₖ - Γₖ) = √NMₖ + oP(1)
√N(B̂ - B) = √NMB + oP(1)

推论4.3 (渐近正态性)

√N(ŜN - L) →^L N(0, τ²)

检验程序

在显著性水平α下，当以下条件成立时拒绝原假设：

√NŜN/τ̂ > z₁₋α

其中τ̂²是τ²的一致估计量。

理论保证

  1. 正确渐近水平: 检验程序在原假设下达到正确的α水平
  2. 一致性: 对任何固定备择假设H₁: L > 0，有√NŜN →^P +∞
  3. 相对测度: 引入标准化统计量G = L/B ∈ 0,1，便于解释

相关工作

McKean-Vlasov估计理论

  • Amorino等(2024): 非参数估计的多项式收敛率
  • Belomestny等(2022): 半参数估计
  • Comte & Genon-Catalot(2024): 参数推断
  • Della Maestra & Hoffmann(2022): 非参数估计

经典SDE检验方法

  • Dette & Podolskij(2008): 经典扩散模型的波动率检验
  • Ait-Sahalia(1996): 连续时间模型检验
  • Corradi & White(1999): 扩散方差的规格检验

本文相对优势

  1. 首次处理McKean-Vlasov: 现有方法仅适用于经典SDE
  2. 分布依赖性: 能够处理波动率对整个分布的依赖
  3. 双重渐近: 联合考虑高频和大样本渐近

结论与讨论

主要结论

  1. 成功建立了McKean-Vlasov模型波动率函数拟合优度检验的完整统计理论
  2. 证明了在双重渐近框架下估计量的一致性和检验统计量的渐近正态性
  3. 构造了具有正确渐近性质的实用检验程序

局限性

  1. 理论性工作: 缺乏数值实验验证理论结果的有效性
  2. 假设条件: 需要较强的光滑性和矩条件假设
  3. 计算复杂性: 实际实施中需要计算泛函导数，可能较为复杂
  4. 有限样本性质: 未提供有限样本下的性能分析

未来方向

  1. 数值验证: 通过蒙特卡洛模拟验证理论结果
  2. 实际应用: 在金融数据上测试方法的实用性
  3. 扩展: 推广到多维情形和更一般的相互作用结构
  4. 计算优化: 开发高效的数值算法

深度评价

优点

  1. 理论严谨性: 数学证明完整，技术处理精细，特别是对分布依赖性的处理创新性强
  2. 问题重要性: 填补了McKean-Vlasov模型统计检验的重要空白
  3. 方法创新: 巧妙结合了泛函分析、随机过程理论和高频统计
  4. 实用价值: 提供了可实施的检验程序，具有良好的渐近性质

不足

  1. 缺乏数值验证: 纯理论工作，未提供仿真实验验证
  2. 假设较强: 线性泛函导数等假设在实际中可能难以验证
  3. 计算挑战: 实际实施可能面临计算复杂性问题
  4. 应用指导不足: 缺乏实际应用的具体指导

影响力

  1. 学术贡献: 在McKean-Vlasov统计推断领域具有开创性意义
  2. 理论价值: 为后续研究提供了重要的理论基础
  3. 应用潜力: 在金融风险管理、系统性风险建模等领域有应用前景

适用场景

  1. 金融建模: 系统性风险、均值场博弈模型的验证
  2. 物理系统: 大规模相互作用粒子系统的建模验证
  3. 社会科学: 群体行为模型的统计检验
  4. 工程应用: 复杂网络系统的动力学建模验证

参考文献

论文引用了30篇相关文献，主要包括：

  • McKean-Vlasov理论基础文献（Sznitman, 1991; Carmona & Delarue, 2018）
  • 统计估计方法（Amorino等, 2024; Belomestny等, 2022）
  • 经典SDE检验方法（Dette & Podolskij, 2008; Corradi & White, 1999）
  • 高频统计理论（Barndorff-Nielsen等, 2006）

本论文在McKean-Vlasov随机微分方程的统计检验领域做出了重要的理论贡献，为这一新兴交叉领域提供了坚实的数学基础。虽然缺乏数值验证，但其理论框架的建立为后续的应用研究奠定了基础。