[論文レビュー] 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).

研究の動機と目的

• スケーリングパラメータ $ h \to 0 $ のとき、拡大する多様体上での非一様・局所定常ガウス場の上限の極限分布を導出すること。
• 区間または立方体上でのガウス過程に対する古典的極値結果を、滑らかでコンパクトな $ r $ 次元多様体の状況に一般化すること。
• 非パラメトリック統計におけるリッジ、等高線集合、密度境界などの幾何的対象に対する一様推論の理論的基盤を提供すること。
• スケーリングパrameter $ h \to 0 $ に従ってスケーリングされる多様体上での、領域を越えた上昇確率の均一漸近的結果を確立し、信頼区間の構築に役立つこと。
• 局所的等方性を持つ非定常・多様体インデックス付きガウス場への極値理論の適用範囲を拡張すること。

提案手法

• 中心化され、連続的で、局所定常なガウス場 $ \{X_h(t), t \in \mathcal{M}_h\} $ のクラスを形式化する。ここで $ \mathcal{M}_h = \frac{1}{h}\mathcal{M} $ はスケーリングされた滑らかな多様体である。
• 多様体上の各点の周囲での共分散構造の局所的展開を用い、行列関数 $ D_t $ を用いて局所定常性をモデル化する。
• 多様体 $ \mathcal{M}_h $ の離散化に三角形分割とメッシュ化技術を適用し、連続場の上限を離散グリッド上の最大値に近似する。
• 共分散の $ \alpha $-正則減衰 $ r(t) = 1 - \|t\|^\alpha + o(\|t\|^\alpha) $ を持つガウス場に、Pickands型極値指数とPickands–Berman法を適用する。
• Piterbarg（1996）およびMikhaleva–Piterbarg（1997）の固定多様体上での上昇確率に関する極限定理を、スケーリングによって拡大された状況に拡張する。
• 有界なハウスドルフ測度を持つジョルダン可測集合の族全体にわたる一様収束を確立し、多様体の形状やサイズにかかわらず頑健性を確保する。

実験結果

リサーチクエスチョン

• RQ1 スケーリングパラメータ $ h \to 0 $ のとき、スケーリングされた多様体上での局所定常ガウス場の上限の極限分布は何か？
• RQ2 このような場の極値挙動は、その背後にある多様体の局所的幾何学的構造と共分散構造にどのように依存するか？
• RQ3 区間上でのカーネル密度推定器の上限に関する古典的Bickel–Rosenblatt結果は、多様体にインデックスされた場へ一般化可能か？
• RQ4 上昇確率 $ \mathbb{P}(\sup_{t \in \mathcal{M}_h} X_h(t) > x) $ の漸近的挙動は $ h \to 0 $ のときどのように変化し、多様体の次元と曲率にどのように依存するか？
• RQ5 正規化された上限分布の極限が非退化なWeibull型分布に収束するための条件は何か？

主な発見

• 極限における正規化された上限分布は、$ x \to \infty $ のとき $ \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 $ と漸近的に一致する。ここで $ r $ は多様体の次元を表す。
• 漸近的スケーリング係数 $ x^{2r/\alpha} \Psi(x) $ は、共分散の局所的ホルダー正則性 $ \alpha $ と多様体の内因的次元 $ r $ の間の相互作用を捉えている。
• 極限強度は、多様体上での局所スケーリング行列 $ D_s $ の積分によって決定され、その重みは点 $ s $ における接空間の方向性を表す $ M_s $（接空間の正規直交基底行列）によって与えられる。
• 収束は、有界な $ r $ 次元ハウスドルフ測度を持つジョルダン可測 $ r $ 次元集合の族にわたって一様であるため、統計的推論において頑健性が保証される。
• 本結果は、区間から多様体へのBickelとRosenblatt（1973a）の結果の一般化であり、固定多様体から成長する多様体へのMikhalevaとPiterbarg（1997）の結果の拡張である。
• 極値極限は非退化的であり、Weibull型である。正規化定数 $ a_h $ と $ b_h $ は（1.1）式において、$ x^{2r/\alpha} \Psi(x) $ のスケーリングによって暗黙的に定義される。