[论文解读] Distribution of the Height of Local Maxima of Gaussian Random Fields

Dan Cheng, Armin Schwartzman|arXiv (Cornell University)|Jul 22, 2013

Morphological variations and asymmetry参考文献 16被引用 4

一句话总结

本文推导了在欧几里得空间和黎曼流形上的非平稳各向同性高斯随机场中局部极大值高度分布的精确公式与渐近公式。通过利用局部几何与随机矩阵理论（特别是高斯正交系综，GOE），该文给出了峰值高度条件概率的闭式表达式，从而在无需平稳性或全局场性质假设的前提下，实现了峰值检测中的精确统计推断。

ABSTRACT

Let $\{f(t): t\in T\}$ be a smooth Gaussian random field over a parameter space $T$, where $T$ may be a subset of Euclidean space or, more generally, a Riemannian manifold. For any local maximum of $f(t)$ located at $t_0$ in the interior of $T$, we provide general formulae and asymptotic approximations for both the tail distribution of the height of a local maximum $\mathbb{P}\{f(t_0)>u | t_0 ext{is a local maximum of} f(t) \}$ and the overshoot distribution of a local maximum $\mathbb{P}\{f(t_0)>u+v | t_0 ext{is a local maximum of} f(t) ext{and} f(t_0)>v\}$. Assuming further that $f$ is isotropic, we apply techniques from random matrix theory related to the Gaussian orthogonal ensemble to compute such conditional probabilities explicitly when $T$ is Euclidean or a sphere of arbitrary dimension. Such calculations are motivated by the statistical problem of detecting peaks in the presence of smooth Gaussian noise.

研究动机与目标

• 推导定义在黎曼流形或欧几里得空间上的非平稳高斯随机场中局部极大值高度的条件分布的一般公式。
• 在仅依赖于局部场性质的前提下，建立当预阈值较高时局部极大值超额分布的渐近近似。
• 将现有平稳场结果推广至非平稳且各向同性的场，包括任意维度的球面上的情形。
• 通过应用随机矩阵理论（特别是高斯正交系综，GOE）的工具，为各向同性场下高度分布提供精确的、非渐近的闭式表达式。

提出的方法

• 通过点附近邻域的极限定义条件概率，来定义局部极大值的高度分布，从而避免零概率的条件事件。
• 基于场在某一点的联合分布（场值、梯度与海森矩阵）进行局部分析，通过指标与行列式条件来刻画局部极大值。
• 利用高斯正交系综（GOE）对各向同性下海森矩阵的分布进行建模，从而实现对海森矩阵行列式期望的精确计算。
• 通过高斯向量的条件矩推导关键表达式，特别是给定场值与零梯度时海森矩阵的分布。
• 利用随机矩阵理论计算海森矩阵特征值函数的期望，尤其关注在指标约束下的绝对行列式。
• 通过正交基与基于坐标图的论证，将结果从欧几里得空间推广至黎曼流形，适应局部几何结构。

实验结果

研究问题

• RQ1 在黎曼流形上的非平稳高斯随机场中，局部极大值高度的精确分布是什么？
• RQ2 当场值超过一个高预阈值时，如何对局部极大值的超额分布进行渐近近似？
• RQ3 能否利用随机矩阵理论，对球面或欧几里得空间上的各向同性高斯场，以闭式表达计算局部极大值的高度分布？
• RQ4 极限超额分布是否依赖于场的相关结构，还是在各向同性下具有普适性？
• RQ5 在保持局部依赖性的同时，结果能多大程度地从欧几里得空间推广至一般黎曼流形？

主要发现

• 推导出非平稳高斯场在黎曼流形上局部极大值高度分布的一般公式，仅依赖于局部场性质（如邻域内局部极大值的期望数量）。
• 对于球面或欧几里得空间上的各向同性高斯场，局部极大值高度分布以高斯正交系综（GOE）形式给出闭式表达，明确依赖于场的方差与相关参数。
• 当预阈值 $ v \to \bigvee $ 时，超额分布渐近趋近于一个与相关函数无关的分布，且在各向同性下近似误差为超指数级小。
• 在各向同性下，给定场值与零梯度时，海森矩阵的条件分布被证明等价于一个平移后的GOE矩阵，从而可精确计算行列式期望。
• 通过水平 $ v $ 上的上穿集的欧拉示性数近似邻域内局部极大值的期望数量，从而支持超额分布的渐近分析。
• 通过正交基调整局部几何，将结果推广至一般黎曼流形，证明依赖于保持原始欧几里得推导结构的基于坐标图的论证。