[论文解读] Extrema of locally stationary Gaussian fields on growing manifolds

Wanli Qiao, Wolfgang Polonik|arXiv (Cornell University)|Oct 23, 2015

Hydrology and Drought Analysis参考文献 27被引用 3

一句话总结

本文建立了在缩放参数 $ h \to 0 $ 时，对在不断增长的流形上定义的局部平稳高斯场，其极值的渐近分布。通过推广经典结果 Bickel 和 Rosenblatt（1973a）以及 Mikhaleva 和 Piterbarg（1997），推导出此类场的上确界极限定理，表明其收敛至一种由流形维度和局部协方差结构决定几何强度的 Weibull 型极值分布。

ABSTRACT

We consider a class of non-homogeneous, continuous, centered Gaussian random fields $\{X_h(t), t \in {\cal M}_h;\,0 < h \le 1\}$ where ${\cal M}_h$ denotes a rescaled smooth manifold, i.e. ${\cal M}_h = \frac{1}{h} {\cal M},$ and study the limit behavior of the extreme values of these Gaussian random fields when $h$ tends to zero, which means that the manifold is growing. Our main result can be thought of as a generalization of a classical result of Bickel and Rosenblatt (1973a), and also of results by Mikhaleva and Piterbarg (1997).

研究动机与目标

• 推导当流形增长时，非平稳、局部平稳高斯场在缩放流形上的上确界极限分布。
• 将高斯过程在区间或立方体上的经典极值结果推广至光滑、紧致的 $ r $-维流形情形。
• 为非参数统计中几何对象（如脊线、等值集和密度边界）的统一推断提供理论基础。
• 建立在 $ h \to 0 $ 缩放下，流形上泛函上穿概率的统一渐近结果，适用于置信带构造。
• 将极值理论的应用范围扩展至具有局部各向同性性质的非平稳、流形索引高斯场。

提出的方法

• 形式化定义一类均值为零、连续且局部平稳的高斯场 $ \{X_h(t), t \in \mathcal{M}_h\} $，其中 $ \mathcal{M}_h = \frac{1}{h}\mathcal{M} $ 为缩放后的光滑流形。
• 通过在流形上每一点附近对协方差结构进行局部展开，利用矩阵函数 $ D_t $ 建模局部平稳性。
• 采用三角剖分与网格化技术对缩放流形 $ \mathcal{M}_h $ 进行离散化，将连续场的上确界近似为离散网格上最大值的上确界。
• 利用 Pickands 型极值指数及 Pickands–Berman 方法处理具有 $ \alpha $-正则协方差衰减 $ r(t) = 1 - \|t\|^\alpha + o(\|t\|^\alpha) $ 的高斯场。
• 利用 Piterbarg（1996）以及 Mikhaleva–Piterbarg（1997）关于固定流形上穿概率的极限结果，并通过缩放将其推广至增长情形。
• 在具有有界 $ r $-维 Hausdorff 测度的 Jordan 可测集族上建立一致收敛性，确保对不同形状与大小的流形具有鲁棒性。

实验结果

研究问题

• RQ1 当缩放参数 $ h \to 0 $ 时，局部平稳高斯场在缩放流形上的上确界极限分布为何？
• RQ2 此类场的极值行为如何依赖于底层流形的局部几何与协方差结构？
• RQ3 经典 Bickel–Rosenblatt 结果（关于核密度估计器上确界的结论）能否推广至流形索引场？
• RQ4 当 $ h \to 0 $ 时，上穿概率 $ \mathbb{P}(\sup_{t \in \mathcal{M}_h} X_h(t) > x) $ 的渐近行为如何？其与流形维度和曲率的关系为何？
• RQ5 在何种条件下，归一化上确界分布的极限收敛至非退化的 Weibull 型分布？

主要发现

• 归一化上确界分布的极限为 $ \mathbb{P}(\sup_{t \in \mathcal{M}} X(t) > x) \sim x^{2r/\alpha} \Psi(x) \cdot H_\alpha^{(r)} \int_{\mathcal{M}} \|D_s M_s\|_r \, ds $，当 $ x \to \infty $ 时成立，其中 $ r $ 为流形的维度。
• 渐近缩放因子 $ x^{2r/\alpha} \Psi(x) $ 反映了协方差的局部 Hölder 正则性 $ \alpha $ 与流形内在维度 $ r $ 之间的相互作用。
• 极限强度由流形上局部缩放矩阵 $ D_s $ 的积分决定，其权重由 $ M_s $（即点 $ s $ 处切空间的正交基矩阵）所体现的切空间方向决定。
• 收敛性在具有有界 $ r $-维 Hausdorff 测度的 Jordan 可测 $ r $-维集合族上一致成立，确保了统计推断的鲁棒性。
• 该结果将 Bickel 和 Rosenblatt（1973a）从区间推广至流形，同时将 Mikhaleva 和 Piterbarg（1997）从固定流形推广至增长流形。
• 极值极限为非退化的 Weibull 型分布，且（1.1）式中的归一化常数 $ a_h $ 与 $ b_h $ 隐含地由 $ x^{2r/\alpha} \Psi(x) $ 的缩放方式决定。