We study the local asymptotic normality (LAN) property for the likelihood function associated with discretely observed $d$-dimensional McKean-Vlasov stochastic differential equations over a fixed time interval. The model involves a joint parameter in both the drift and diffusion coefficients, introducing challenges due to its dependence on the process distribution. We derive a stochastic expansion of the log-likelihood ratio using Malliavin calculus techniques and establish the LAN property under appropriate conditions. The main technical challenge arises from the implicit nature of the transition densities, which we address through integration by parts and Gaussian-type bounds. This work extends existing LAN results for interacting particle systems to the mean-field regime, contributing to statistical inference in non-linear stochastic models
论文ID : 2511.13366标题 : Local asymptotic normality for discretely observed McKean-Vlasov diffusions作者 : Akram Heidari, Mark Podolskij (University of Luxembourg)分类 : math.ST, stat.TH (统计理论)提交时间 : 2025年11月17日论文链接 : https://arxiv.org/abs/2511.13366 本文研究了在固定时间区间内离散观测的d维McKean-Vlasov随机微分方程的似然函数的局部渐近正态性(LAN)性质。模型在漂移系数和扩散系数中都包含联合参数，由于对过程分布的依赖性引入了挑战。作者使用Malliavin微积分技术导出对数似然比的随机展开，并在适当条件下建立LAN性质。主要技术挑战源于转移密度的隐式性质，通过分部积分和高斯型界来解决。这项工作将现有的交互粒子系统的LAN结果扩展到平均场机制，为非线性随机模型的统计推断做出贡献。
本文研究McKean-Vlasov随机微分方程(SDEs)的参数估计问题，具体针对离散时间观测数据建立似然函数的局部渐近正态性(LAN)。McKean-Vlasov方程形式为：
d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i dX^{i,\theta}_t = b_{\theta_1}(X^{i,\theta}_t, \mu^\theta_t)dt + a_{\theta_2}(X^{i,\theta}_t)dW^i_t d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i
其中μ t θ \mu^\theta_t μ t θ X t i , θ X^{i,\theta}_t X t i , θ
应用广泛性 : McKean-Vlasov方程在统计物理、金融、平均场博弈等领域有广泛应用理论基础 : LAN性质是渐近统计推断的基础工具，由Le Cam引入，可用于推导估计量渐近方差的下界平均场理论 : 连接微观粒子系统与宏观平均场极限的统计推断连续观测 vs 离散观测 : 现有LAN结果主要针对连续观测情形13 ，此时可通过Girsanov定理获得似然函数的闭式表达式转移密度的隐式性 : 离散观测情形下，转移密度没有显式表达，需要新的技术处理交互粒子系统的挑战 : 对于交互粒子系统(3.18)，需要处理dN维联合转移密度，目前文献中缺乏相应的界填补离散观测McKean-Vlasov方程LAN理论的空白 发展处理隐式转移密度的Malliavin微积分技术 为平均场模型的统计推断提供理论基础，与最近的估计方法1 建立理论联系 建立LAN性质 : 首次为离散观测的McKean-Vlasov方程建立LAN性质，在渐近机制Δ n → 0 , N → ∞ \Delta_n \to 0, N \to \infty Δ n → 0 , N → ∞ T T T Malliavin微积分技术 : 使用Malliavin微积分中的分部积分公式导出转移密度对数导数的显式表示(命题3.1)随机展开 : 建立对数似然比的精确随机展开(命题3.2)，识别主项和余项渐近协方差矩阵 : 给出显式的渐近协方差矩阵Σ θ 0 \Sigma^{\theta_0} Σ θ 0 ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 不同估计速率 : 证明漂移参数以速率N \sqrt{N} N N / Δ n \sqrt{N/\Delta_n} N / Δ n 1 中的对比估计方法一致技术创新 : 通过高斯型上下界(命题4.2)和矩估计(引理4.1)处理隐式转移密度，克服了主要技术障碍观测数据 :
{ X t k i , θ } i = 1 , … , N k = 1 , … , n \{X^{i,\theta}_{t_k}\}_{i=1,\ldots,N}^{k=1,\ldots,n} { X t k i , θ } i = 1 , … , N k = 1 , … , n t k = T k / n t_k = Tk/n t k = T k / n Δ n = T / n \Delta_n = T/n Δ n = T / n
参数扰动 :
θ + = ( θ 1 + , θ 2 + ) = ( θ 1 0 + u N , θ 2 0 + v N / Δ n ) \theta^+ = (\theta_1^+, \theta_2^+) = \left(\theta_1^0 + \frac{u}{\sqrt{N}}, \theta_2^0 + \frac{v}{\sqrt{N/\Delta_n}}\right) θ + = ( θ 1 + , θ 2 + ) = ( θ 1 0 + N u , θ 2 0 + N / Δ n v )
目标 : 证明对数似然比
z ( θ 0 , θ + ) : = log d P θ + d P θ 0 z(\theta_0, \theta^+) := \log \frac{dP_{\theta^+}}{dP_{\theta_0}} z ( θ 0 , θ + ) := log d P θ 0 d P θ + z ( θ 0 , θ + ) → l a w ( u v ) ⊤ N θ 0 − 1 2 ( u v ) ⊤ Σ θ 0 ( u v ) z(\theta_0, \theta^+) \xrightarrow{law} \begin{pmatrix} u \\ v \end{pmatrix}^\top N_{\theta_0} - \frac{1}{2}\begin{pmatrix} u \\ v \end{pmatrix}^\top \Sigma_{\theta_0} \begin{pmatrix} u \\ v \end{pmatrix} z ( θ 0 , θ + ) l a w ( u v ) ⊤ N θ 0 − 2 1 ( u v ) ⊤ Σ θ 0 ( u v )
其中N θ 0 ∼ N ( 0 , Σ θ 0 ) N_{\theta_0} \sim N(0, \Sigma_{\theta_0}) N θ 0 ∼ N ( 0 , Σ θ 0 )
模型假设N N N d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i dX^{i,\theta}_t = b_{\theta_1}(X^{i,\theta}_t, \mu^\theta_t)dt + a_{\theta_2}(X^{i,\theta}_t)dW^i_t d X t i , θ = b θ 1 ( X t i , θ , μ t θ ) d t + a θ 2 ( X t i , θ ) d W t i
关键特征：
分布依赖 : 漂移项依赖于边际分布μ t θ = Law ( X t i , θ ) \mu^\theta_t = \text{Law}(X^{i,\theta}_t) μ t θ = Law ( X t i , θ ) 参数分离 : 漂移参数θ 1 \theta_1 θ 1 θ 2 \theta_2 θ 2 独立性 : 不同粒子的布朗运动( W i ) 1 ≤ i ≤ N (W^i)_{1\leq i \leq N} ( W i ) 1 ≤ i ≤ N 利用Markov性质：
z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N log p θ + p θ 0 ( t k , t k + 1 , X t k i , X t k + 1 i ) z(\theta_0, \theta^+) = \sum_{k=1}^n \sum_{i=1}^N \log \frac{p^{\theta^+}}{p^{\theta_0}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}}) z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N log p θ 0 p θ + ( t k , t k + 1 , X t k i , X t k + 1 i )
进一步分解为漂移和扩散两部分：
z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N ( ζ k i , θ 1 + ζ k i , θ 2 ) z(\theta_0, \theta^+) = \sum_{k=1}^n \sum_{i=1}^N (\zeta^{i,\theta_1}_k + \zeta^{i,\theta_2}_k) z ( θ 0 , θ + ) = ∑ k = 1 n ∑ i = 1 N ( ζ k i , θ 1 + ζ k i , θ 2 )
其中
ζ k i , θ 1 = u N ∫ 0 1 ∂ θ 1 p θ 1 ( l ) , θ 2 + p θ 1 ( l ) , θ 2 + ( t k , t k + 1 , X t k i , X t k + 1 i ) d l \zeta^{i,\theta_1}_k = \frac{u}{\sqrt{N}} \int_0^1 \frac{\partial_{\theta_1} p^{\theta_1(l), \theta_2^+}}{p^{\theta_1(l), \theta_2^+}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}})dl ζ k i , θ 1 = N u ∫ 0 1 p θ 1 ( l ) , θ 2 + ∂ θ 1 p θ 1 ( l ) , θ 2 + ( t k , t k + 1 , X t k i , X t k + 1 i ) d l
ζ k i , θ 2 = v N / Δ n ∫ 0 1 ∂ θ 2 p θ 1 0 , θ 2 ( l ) p θ 1 0 , θ 2 ( l ) ( t k , t k + 1 , X t k i , X t k + 1 i ) d l \zeta^{i,\theta_2}_k = \frac{v}{\sqrt{N/\Delta_n}} \int_0^1 \frac{\partial_{\theta_2} p^{\theta_1^0, \theta_2(l)}}{p^{\theta_1^0, \theta_2(l)}}(t_k, t_{k+1}, X^i_{t_k}, X^i_{t_{k+1}})dl ζ k i , θ 2 = N / Δ n v ∫ 0 1 p θ 1 0 , θ 2 ( l ) ∂ θ 2 p θ 1 0 , θ 2 ( l ) ( t k , t k + 1 , X t k i , X t k + 1 i ) d l
关键技术 : 对于转移密度的参数导数，使用Malliavin微积分得到：
∂ θ 1 p θ p θ ( t k , t k + 1 , x , y ) = 1 Δ n E t k , x θ [ ∑ r = 1 d δ ( ∂ θ 1 X r , Δ n i , θ U r i ) ∣ X t k + 1 i , θ = y ] \frac{\partial_{\theta_1} p^\theta}{p^\theta}(t_k, t_{k+1}, x, y) = \frac{1}{\Delta_n} E^\theta_{t_k,x}\left[\sum_{r=1}^d \delta(\partial_{\theta_1} X^{i,\theta}_{r,\Delta_n} U^i_r) \Big| X^{i,\theta}_{t_{k+1}} = y\right] p θ ∂ θ 1 p θ ( t k , t k + 1 , x , y ) = Δ n 1 E t k , x θ [ ∑ r = 1 d δ ( ∂ θ 1 X r , Δ n i , θ U r i ) X t k + 1 i , θ = y ]
其中：
δ \delta δ U s i = a θ 2 − 1 ( X t k + s i , θ ) Y s i , θ ( Y Δ n i , θ ) − 1 U^i_s = a^{-1}_{\theta_2}(X^{i,\theta}_{t_k+s}) Y^{i,\theta}_s (Y^{i,\theta}_{\Delta_n})^{-1} U s i = a θ 2 − 1 ( X t k + s i , θ ) Y s i , θ ( Y Δ n i , θ ) − 1 Y t i , θ Y^{i,\theta}_t Y t i , θ 参数导数过程 : ∂ θ 1 X t i , θ \partial_{\theta_1} X^{i,\theta}_t ∂ θ 1 X t i , θ ∂ θ 1 X t i , θ = ∫ 0 t ( ∂ θ 1 b θ 1 + ∇ x b θ 1 ∂ θ 1 X s i , θ + ∫ R d ∂ μ b θ 1 ( X t k + s i , θ , y , μ t k + s θ ) ∂ θ 1 μ t k + s θ ( d y ) ) d s + … \partial_{\theta_1} X^{i,\theta}_t = \int_0^t \left(\partial_{\theta_1} b_{\theta_1} + \nabla_x b_{\theta_1} \partial_{\theta_1} X^{i,\theta}_s + \int_{\mathbb{R}^d} \partial_\mu b_{\theta_1}(X^{i,\theta}_{t_k+s}, y, \mu^\theta_{t_k+s}) \partial_{\theta_1}\mu^\theta_{t_k+s}(dy)\right)ds + \ldots ∂ θ 1 X t i , θ = ∫ 0 t ( ∂ θ 1 b θ 1 + ∇ x b θ 1 ∂ θ 1 X s i , θ + ∫ R d ∂ μ b θ 1 ( X t k + s i , θ , y , μ t k + s θ ) ∂ θ 1 μ t k + s θ ( d y ) ) d s + …
注意第三项包含函数导数 ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1
漂移部分 : 证明
δ ( ∂ θ 1 X r , Δ n i , θ U r i ) = Δ n z r , θ 1 θ ( X t k i , θ ) [ a θ 2 − 2 ( X t k i , θ ) ( X t k + 1 i , θ − m t k , t k + 1 θ ( X t k i , θ ) ) ] r + H t k + 1 i \delta(\partial_{\theta_1} X^{i,\theta}_{r,\Delta_n} U^i_r) = \Delta_n z^{\theta}_{r,\theta_1}(X^{i,\theta}_{t_k}) [a^{-2}_{\theta_2}(X^{i,\theta}_{t_k})(X^{i,\theta}_{t_{k+1}} - m^\theta_{t_k,t_{k+1}}(X^{i,\theta}_{t_k}))]_r + H^i_{t_{k+1}} δ ( ∂ θ 1 X r , Δ n i , θ U r i ) = Δ n z r , θ 1 θ ( X t k i , θ ) [ a θ 2 − 2 ( X t k i , θ ) ( X t k + 1 i , θ − m t k , t k + 1 θ ( X t k i , θ )) ] r + H t k + 1 i
其中H t k + 1 i H^i_{t_{k+1}} H t k + 1 i ( E t k , x θ ∣ H t k + 1 i ∣ τ ) 1 / τ = R t k i ( Δ n 2 ) (E^\theta_{t_k,x}|H^i_{t_{k+1}}|^\tau)^{1/\tau} = R^i_{t_k}(\Delta_n^2) ( E t k , x θ ∣ H t k + 1 i ∣ τ ) 1/ τ = R t k i ( Δ n 2 )
关键量 :
z t θ ( x ) : = ∂ θ 1 b θ 1 ( x , μ t θ ) + ∫ R d ∂ μ b θ 1 ( x , y , μ t θ ) ∂ θ 1 μ t θ ( d y ) z^\theta_t(x) := \partial_{\theta_1} b_{\theta_1}(x, \mu^\theta_t) + \int_{\mathbb{R}^d} \partial_\mu b_{\theta_1}(x, y, \mu^\theta_t) \partial_{\theta_1}\mu^\theta_t(dy) z t θ ( x ) := ∂ θ 1 b θ 1 ( x , μ t θ ) + ∫ R d ∂ μ b θ 1 ( x , y , μ t θ ) ∂ θ 1 μ t θ ( d y )
这个量在渐近协方差矩阵中起核心作用。
技术路径 :
使用分部积分公式(2.7): δ ( F u ) = F δ ( u ) − ⟨ D F , u ⟩ H \delta(Fu) = F\delta(u) - \langle DF, u\rangle_H δ ( F u ) = F δ ( u ) − ⟨ D F , u ⟩ H 将U r i U^i_r U r i U ^ r i = a θ 2 − 1 ( X t k + r i , θ ) \hat{U}^i_r = a^{-1}_{\theta_2}(X^{i,\theta}_{t_k+r}) U ^ r i = a θ 2 − 1 ( X t k + r i , θ ) 证明各余项H n i , j , j = 1 , 2 , 3 H^{i,j}_n, j=1,2,3 H n i , j , j = 1 , 2 , 3 Δ n 2 \Delta_n^2 Δ n 2 扩散部分 : 类似地证明
δ ( ∂ θ 2 X r , Δ n i , θ U r i ) = [ ∂ θ 2 a θ 2 ( X t k i ) a θ 2 − 1 ( X t k i ) ( X t k + 1 i − m t k , t k + 1 θ ( X t k i ) ) ] r × [ … ] + 余项 \delta(\partial_{\theta_2} X^{i,\theta}_{r,\Delta_n} U^i_r) = [\partial_{\theta_2} a_{\theta_2}(X^i_{t_k}) a^{-1}_{\theta_2}(X^i_{t_k})(X^i_{t_{k+1}} - m^\theta_{t_k,t_{k+1}}(X^i_{t_k}))]_r \times [\ldots] + \text{余项} δ ( ∂ θ 2 X r , Δ n i , θ U r i ) = [ ∂ θ 2 a θ 2 ( X t k i ) a θ 2 − 1 ( X t k i ) ( X t k + 1 i − m t k , t k + 1 θ ( X t k i )) ] r × [ … ] + 余项
余项满足Δ n 3 / 2 \Delta_n^{3/2} Δ n 3/2
Gobet方法19,20 : 原始方法针对遍历扩散过程，依赖长时间观测的遍历性本文方法 : 不需要遍历性假设，渐近性由粒子数N → ∞ N \to \infty N → ∞ 函数导数 : McKean-Vlasov模型中出现∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 Aronson型上下界 :
1 L Δ n d / 2 exp ( − c ∥ x − y ∥ 2 Δ n ) exp ( − c Δ n ∥ x ∥ 2 ) ≤ p θ ( t k , t k + 1 , x , y ) \frac{1}{L\Delta_n^{d/2}} \exp\left(-c\frac{\|x-y\|^2}{\Delta_n}\right) \exp(-c\Delta_n\|x\|^2) \leq p^\theta(t_k, t_{k+1}, x, y) L Δ n d /2 1 exp ( − c Δ n ∥ x − y ∥ 2 ) exp ( − c Δ n ∥ x ∥ 2 ) ≤ p θ ( t k , t k + 1 , x , y ) ≤ L Δ n d / 2 exp ( − ∥ x − y ∥ 2 c Δ n ) exp ( c Δ n ∥ x ∥ 2 ) \leq \frac{L}{\Delta_n^{d/2}} \exp\left(-\frac{\|x-y\|^2}{c\Delta_n}\right) \exp(c\Delta_n\|x\|^2) ≤ Δ n d /2 L exp ( − c Δ n ∥ x − y ∥ 2 ) exp ( c Δ n ∥ x ∥ 2 )
参数导数界 :
E t k , x θ ˉ [ ∣ ∂ θ 1 p θ p θ ( t k , t k + 1 , x , X t k + 1 i ) ∣ m ] ≤ L Δ n m / 2 exp ( c Δ n ∥ x ∥ 2 ) ( 1 + ∥ x ∥ ) q E^{\bar{\theta}}_{t_k,x}\left[\left|\frac{\partial_{\theta_1} p^\theta}{p^\theta}(t_k, t_{k+1}, x, X^i_{t_{k+1}})\right|^m\right] \leq \frac{L}{\Delta_n^{m/2}} \exp(c\Delta_n\|x\|^2)(1+\|x\|)^q E t k , x θ ˉ [ p θ ∂ θ 1 p θ ( t k , t k + 1 , x , X t k + 1 i ) m ] ≤ Δ n m /2 L exp ( c Δ n ∥ x ∥ 2 ) ( 1 + ∥ x ∥ ) q
这些界对于证明余项的可忽略性（命题4.4）至关重要。
注 : 本文是纯理论论文，不包含数值实验。主要结果是理论定理的证明。
虽然没有数值实验，但论文通过以下方式验证理论的合理性：
与现有结果的一致性 : 渐近协方差矩阵Σ θ 0 \Sigma^{\theta_0} Σ θ 0 N Δ n → 0 N\Delta_n \to 0 N Δ n → 0 1 中的对比估计方法的渐近方差一致估计速率 :漂移参数: N \sqrt{N} N 扩散参数: N / Δ n \sqrt{N/\Delta_n} N / Δ n 1 一致特殊情况 : 当∂ μ b θ 1 = 0 \partial_\mu b_{\theta_1} = 0 ∂ μ b θ 1 = 0 LAN性质 : 在假设A1-A5下，
z ( θ 0 , θ + ) → P θ 0 − l a w ( u v ) ⊤ N θ 0 − 1 2 ( u v ) ⊤ Σ θ 0 ( u v ) z(\theta_0, \theta^+) \xrightarrow{P^{\theta_0}-law} \begin{pmatrix} u \\ v \end{pmatrix}^\top N_{\theta_0} - \frac{1}{2}\begin{pmatrix} u \\ v \end{pmatrix}^\top \Sigma_{\theta_0} \begin{pmatrix} u \\ v \end{pmatrix} z ( θ 0 , θ + ) P θ 0 − l a w ( u v ) ⊤ N θ 0 − 2 1 ( u v ) ⊤ Σ θ 0 ( u v )
渐近协方差矩阵 :
Σ θ 0 = ( Σ b θ 0 0 0 Σ a θ 0 ) \Sigma_{\theta_0} = \begin{pmatrix} \Sigma^{\theta_0}_b & 0 \\ 0 & \Sigma^{\theta_0}_a \end{pmatrix} Σ θ 0 = ( Σ b θ 0 0 0 Σ a θ 0 )
其中
Σ b θ 0 = ∫ 0 T ∫ R d z s θ 0 ( x ) ⊤ a θ 2 0 − 2 ( x ) z s θ 0 ( x ) μ s θ 0 ( d x ) d s \Sigma^{\theta_0}_b = \int_0^T \int_{\mathbb{R}^d} z^{\theta_0}_s(x)^\top a^{-2}_{\theta_2^0}(x) z^{\theta_0}_s(x) \mu^{\theta_0}_s(dx)ds Σ b θ 0 = ∫ 0 T ∫ R d z s θ 0 ( x ) ⊤ a θ 2 0 − 2 ( x ) z s θ 0 ( x ) μ s θ 0 ( d x ) d s
Σ a θ 0 = 2 ∫ 0 T ∫ R d tr ( ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ) μ s θ 0 ( d x ) d s \Sigma^{\theta_0}_a = 2\int_0^T \int_{\mathbb{R}^d} \text{tr}(\partial_{\theta_2} a_{\theta_2^0}(x) a^{-1}_{\theta_2^0}(x) \partial_{\theta_2} a_{\theta_2^0}(x) a^{-1}_{\theta_2^0}(x)) \mu^{\theta_0}_s(dx)ds Σ a θ 0 = 2 ∫ 0 T ∫ R d tr ( ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x ) ∂ θ 2 a θ 2 0 ( x ) a θ 2 0 − 1 ( x )) μ s θ 0 ( d x ) d s
对角结构 : Σ θ 0 \Sigma_{\theta_0} Σ θ 0 函数导数的作用 : z s θ 0 ( x ) z^{\theta_0}_s(x) z s θ 0 ( x ) ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 与交互粒子系统的差异 :McKean-Vlasov模型(1.1): 协方差包含∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 交互粒子系统(3.18): 协方差中z t θ ( x ) z^\theta_t(x) z t θ ( x ) ∂ θ 1 b θ 1 ( x , μ t θ ) \partial_{\theta_1} b_{\theta_1}(x, \mu^\theta_t) ∂ θ 1 b θ 1 ( x , μ t θ ) 定理3.4的证明通过验证6个收敛条件(4.30)-(4.36)：
条件(4.30) : 漂移部分的一阶矩
∑ k = 1 n ∑ i = 1 N E t k θ 0 [ ζ ^ k i , θ 1 ] → P θ 0 − 1 2 u 2 Σ b θ 0 \sum_{k=1}^n \sum_{i=1}^N E^{\theta_0}_{t_k}[\hat{\zeta}^{i,\theta_1}_k] \xrightarrow{P^{\theta_0}} -\frac{1}{2}u^2 \Sigma^{\theta_0}_b ∑ k = 1 n ∑ i = 1 N E t k θ 0 [ ζ ^ k i , θ 1 ] P θ 0 − 2 1 u 2 Σ b θ 0
关键步骤：利用Taylor展开
m t k , t k + 1 θ 0 ( X t k i ) − m t k , t k + 1 θ 1 ( l ) , θ 2 + ( X t k i ) = − l u Δ n N z t k θ 0 ( X t k i ) + R t k i ( ε n , N Δ n / N ) m^{\theta_0}_{t_k,t_{k+1}}(X^i_{t_k}) - m^{\theta_1(l),\theta_2^+}_{t_k,t_{k+1}}(X^i_{t_k}) = -\frac{lu\Delta_n}{\sqrt{N}} z^{\theta_0}_{t_k}(X^i_{t_k}) + R^i_{t_k}(\varepsilon_{n,N}\Delta_n/\sqrt{N}) m t k , t k + 1 θ 0 ( X t k i ) − m t k , t k + 1 θ 1 ( l ) , θ 2 + ( X t k i ) = − N l u Δ n z t k θ 0 ( X t k i ) + R t k i ( ε n , N Δ n / N )
条件(4.31) : 漂移部分的二阶矩收敛到u 2 Σ b θ 0 u^2\Sigma^{\theta_0}_b u 2 Σ b θ 0
技巧：
交叉项∑ i 1 ≠ i 2 ζ ^ k i 1 , θ 1 ζ ^ k i 2 , θ 1 \sum_{i_1 \neq i_2} \hat{\zeta}^{i_1,\theta_1}_k \hat{\zeta}^{i_2,\theta_1}_k ∑ i 1 = i 2 ζ ^ k i 1 , θ 1 ζ ^ k i 2 , θ 1 主项来自∑ i = 1 N ( ζ ^ k i , θ 1 ) 2 \sum_{i=1}^N (\hat{\zeta}^{i,\theta_1}_k)^2 ∑ i = 1 N ( ζ ^ k i , θ 1 ) 2 条件(4.32) : 四阶矩条件
∑ k = 1 n E t k θ 0 [ ∣ ∑ i = 1 N ζ ^ k i , θ 1 ∣ 4 ] → P θ 0 0 \sum_{k=1}^n E^{\theta_0}_{t_k}\left[\left|\sum_{i=1}^N \hat{\zeta}^{i,\theta_1}_k\right|^4\right] \xrightarrow{P^{\theta_0}} 0 ∑ k = 1 n E t k θ 0 [ ∑ i = 1 N ζ ^ k i , θ 1 4 ] P θ 0 0
通过将四次项分解为不同指标组合（全不同、两对、全相同）并证明每部分都趋于0。
条件(4.33)-(4.35) : 扩散部分的类似条件，使用条件方差的展开
V t k , t k + 1 θ 0 ( x ) − V t k , t k + 1 θ 1 0 , θ 2 ( l ) ( x ) = − 2 l v Δ n 3 / 2 N ∂ θ 2 a θ 2 0 ( X t k i ) a θ 2 0 ( X t k i ) + … V^{\theta_0}_{t_k,t_{k+1}}(x) - V^{\theta_1^0,\theta_2(l)}_{t_k,t_{k+1}}(x) = -\frac{2lv\Delta_n^{3/2}}{\sqrt{N}} \partial_{\theta_2} a_{\theta_2^0}(X^i_{t_k}) a_{\theta_2^0}(X^i_{t_k}) + \ldots V t k , t k + 1 θ 0 ( x ) − V t k , t k + 1 θ 1 0 , θ 2 ( l ) ( x ) = − N 2 l v Δ n 3/2 ∂ θ 2 a θ 2 0 ( X t k i ) a θ 2 0 ( X t k i ) + …
条件(4.36) : 漂移和扩散的渐近独立性，通过证明交叉项渐近为0
离散观测 :
1 Amorino等(2023): 交互粒子系统的对比估计，建立一致性和渐近正态性6 Bishwal(2011): 交互扩散的估计9 Chen(2021): 单轨迹数据的势能最大似然估计16,17 Genon-Catalot & Larédo(2021): 小方差和长时间McKean-Vlasov模型27 Liu & Qiao(2022): 路径依赖McKean-Vlasov SDE31 Sharrock等(2021): 在线参数估计连续观测 :
13 Della Maestra & Hoffmann(2023): 平均场机制下McKean-Vlasov模型的LAN性质（本文直接相关）
区别：连续观测可用Girsanov定理，似然函数有闭式表达 2 Amorino等(2024): 通过反卷积的多项式速率4 Belomestny等(2022): McKean-Vlasov SDE的半参数估计11 Comte等(2024): 非参数矩方法12 Della Maestra & Hoffmann(2022): 交互粒子系统的非参数估计29 Nickl等(2025): McKean-Vlasov模型的贝叶斯非参数推断19,20 Gobet(2001,2002):
椭圆扩散的局部渐近混合正态性 遍历扩散离散观测的LAN性质 本文基础方法 : 使用Malliavin微积分导出转移密度导数的表示vs 13 （连续观测） :vs 1 （对比估计） :提供似然方法的理论基础 建立LAN性质，可导出估计量的渐近最优性 vs 20 （经典SDE） :扩展到McKean-Vlasov设定 不需要遍历性假设 处理函数导数∂ μ b \partial_\mu b ∂ μ b vs 交互粒子系统 :避免高维联合转移密度的界（注释3.3指出这是交互粒子系统LAN的主要障碍） 利用i.i.d.结构简化分析 LAN性质的建立 : 首次为离散观测McKean-Vlasov方程建立LAN性质，填补了该领域的理论空白渐近协方差矩阵的显式形式 :
Σ θ 0 = diag ( Σ b θ 0 , Σ a θ 0 ) \Sigma_{\theta_0} = \text{diag}(\Sigma^{\theta_0}_b, \Sigma^{\theta_0}_a) Σ θ 0 = diag ( Σ b θ 0 , Σ a θ 0 ) ∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 估计速率的确认 :漂移: N \sqrt{N} N 扩散: N / Δ n \sqrt{N/\Delta_n} N / Δ n 1 一致技术贡献 : 发展了处理隐式转移密度的Malliavin微积分技术，结合高斯型界和分部积分假设条件较强 :A2 : 系数有界且Lipschitz连续A3 : 系数的高阶光滑性（C 2 C^2 C 2 A5 : 扩散矩阵的一致椭圆性一维参数限制 : 虽然论文指出可扩展到多参数情形，但只详细处理了θ 1 , θ 2 ∈ R \theta_1, \theta_2 \in \mathbb{R} θ 1 , θ 2 ∈ R 交互粒子系统的空白 :注释3.3指出，对于交互粒子系统(3.18)，由于缺乏dN维转移密度的界，无法建立LAN性质 这是一个重要的开放问题 渐近机制 : 需要Δ n → 0 , N → ∞ \Delta_n \to 0, N \to \infty Δ n → 0 , N → ∞ N Δ n N\Delta_n N Δ n 初始分布 : 假设A1要求初始分布μ 0 \mu_0 μ 0 交互粒子系统的LAN : 发展高维转移密度的界，建立模型(3.18)的LAN性质放松假设条件 :研究非椭圆扩散的情形 允许系数无界或仅局部Lipschitz 多参数扩展 : 完整处理θ 1 ∈ R p , θ 2 ∈ R q \theta_1 \in \mathbb{R}^p, \theta_2 \in \mathbb{R}^q θ 1 ∈ R p , θ 2 ∈ R q 最优估计量构造 : 利用LAN性质构造渐近有效的估计量假设检验 : 基于LAN性质发展McKean-Vlasov模型的假设检验理论非遍历情形 : 扩展到非遍历McKean-Vlasov过程高频数据 : 研究Δ n → 0 \Delta_n \to 0 Δ n → 0 理论严密性 :证明完整详细（第4节占论文一半篇幅） 每个技术步骤都有清晰的论证 适当使用了现代随机分析工具（Malliavin微积分） 方法创新性 :Malliavin微积分的巧妙应用 : 将Skorohod积分展开为主项加余项的技术（命题3.2）是核心创新函数导数的处理 : 正确识别和处理∂ μ b θ 1 \partial_\mu b_{\theta_1} ∂ μ b θ 1 余项控制 : 通过命题4.4统一处理各种余项的可忽略性理论贡献 :填补了离散观测McKean-Vlasov方程LAN理论的空白 连接了似然方法和对比估计方法（与1 的联系） 为McKean-Vlasov模型的渐近统计推断提供理论基础 写作清晰度 :结构清晰：假设→主要结果→证明 记号系统完善（第2.1节） 关键困难和解决思路说明充分（引言和注释3.3） 文献综述全面 : 准确定位本文在McKean-Vlasov统计推断文献中的位置实用性有限 :假设条件较强，实际数据可能不满足 没有数值模拟验证理论结果 没有讨论如何在实践中验证假设 技术细节的可读性 :第4节的证明非常技术化，对非专家不友好 某些关键不等式（如命题4.2的证明）引用了20 但没有详细说明如何适配到McKean-Vlasov设定 结果的局限性 :交互粒子系统(3.18)的LAN性质仍然开放（注释3.3） 只处理了固定时间区间T T T T → ∞ T \to \infty T → ∞ 多参数情形的处理 :虽然声称可扩展到多参数，但只给出了框架性说明 多参数情形的技术细节（特别是Σ θ 0 \Sigma_{\theta_0} Σ θ 0 与应用的脱节 :没有具体的应用例子 没有讨论如何在金融、神经科学等领域使用结果 对领域的贡献 :理论基础 : 为McKean-Vlasov模型的统计推断提供坚实的理论基础方法论 : Malliavin微积分在McKean-Vlasov统计中的系统应用开放问题 : 明确指出交互粒子系统LAN的技术障碍（高维转移密度界），为未来研究指明方向实用价值 :估计量评价 : 可用于评价现有估计量（如1 中的对比估计）的渐近效率下界 : LAN性质提供估计量渐近方差的下界（Cramér-Rao型）最优估计 : 可指导渐近有效估计量的构造可复现性 :✅ 理论结果完全可验证（证明完整） ❌ 没有代码或数值实验 ✅ 假设清晰明确 ⚠️ 某些技术细节需要查阅19,20,30 等参考文献 预期引用情况 :短期 : 平均场统计推断领域的专家会引用中期 : 可能成为McKean-Vlasov统计推断的标准参考长期 : 如果交互粒子系统的问题得到解决，影响力会进一步扩大理论研究 :McKean-Vlasov模型的统计理论 平均场博弈的参数估计 非线性SDE的渐近统计 潜在应用领域 :金融 : 系统性风险模型18 、期权定价21 神经科学 : 神经元网络模型3 统计物理 : 粒子系统的平均场极限社会动力学 : 意见动力学模型8 方法适用性 :✅ 大样本（N N N ✅ 高频观测（Δ n \Delta_n Δ n ✅ 固定时间区间 ✅ 光滑系数 ❌ 小样本或低频观测 ❌ 非椭圆扩散 与其他方法的比较 :vs 对比估计1 : LAN提供理论最优性，对比估计更易计算vs 贝叶斯方法29 : LAN是频率学派，贝叶斯方法更灵活但计算量大vs 非参数方法12 : LAN针对参数模型，非参数方法适用于模型不确定的情况1 Amorino等(2023)13 Della Maestra & Hoffmann(2023)19,20 Gobet(2001,2002)30 Nualart(1995)22,25,26 Le Cam系列本文是McKean-Vlasov随机微分方程统计推断领域的重要理论贡献。通过巧妙应用Malliavin微积分，作者成功建立了离散观测情形下的局部渐近正态性，填补了该领域的理论空白。论文的技术水平高，证明严密，为平均场模型的渐近统计推断提供了坚实的理论基础。
主要价值在于：(1) 理论完整性：系统建立了McKean-Vlasov模型的LAN理论；(2) 方法创新：发展了处理隐式转移密度的技术；(3) 理论指导：为估计量的渐近最优性提供基准。
主要局限：(1) 假设条件较强；(2) 缺乏数值验证；(3) 交互粒子系统的LAN仍然开放。
对于从事McKean-Vlasov模型统计推断的研究者，这是一篇必读文献。对于应用研究者，需要结合具体问题判断假设条件是否满足。
