[論文レビュー] Monte-Carlo method for multiple parametric integrals calculation and solving of linear integral Fredholm equations of a second kind, with confidence regions in uniform norm

E. Ostrovsky, L. Sirota|arXiv (Cornell University)|Jan 27, 2011

Mathematical Approximation and Integration参考文献 8被引用数 14

ひとこと要約

本稿では、複数のパrametric積分および2種類目の線形Fredholm積分方程式を計算するための修正されたモンテカルロ法を提案し、一様ノルムにおける最適収束速度を達成する。エントロピーおよびサブガウス条件の下で、関数中心極限定理を用いて漸近的および非漸近的信頼領域を確立し、核作用素の固有値半径およびメトリックエントロピーに基づく理論的保証を提供する。

ABSTRACT

In this article we offer some modification of Monte-Carlo method for multiple parametric integral computation and solving of a linear integral Fredholm equation of a second kind (well posed problem). We prove that the rate of convergence of offered method is optimal under natural conditions still in the uniform norm, and construct an asymptotical and non-asymptotical confidence region, again in the uniform norm.

研究の動機と目的

• 一様ノルムにおける収束保証を備えた、複数のパラメトリック積分および2種類目の線形Fredholm積分方程式のためのモンテカルロ法の開発。
• 特に一様ノルム（Lpノルムより強い）において自然な仮定の下で最適収束速度を確立すること。
• 解の一点で一様ノルムにおける漸近的および非漸近的信頼領域を構築し、統計的信頼性を保証すること。
• C(T)における関数中心極限定理が成り立つために必要な条件（因子分解、モーメント、エントロピー条件）の必要性を証明すること。
• ボルテラ作用素および準ノルム作用素に関する先行研究を、不連続または滑らかでない核を持つ一般のFredholm方程式へと拡張すること。

提案手法

• モンテカルロサンプリングを用いた積分作用素の確率的近似を採用し、連続関数空間に一様ノルムを備えた収束を解析する。
• 核評価の経験的平均を用いたFredholm方程式の解の修正推定器を導入し、一貫性および確率収束を保証する。
• 核作用素の固有値半径条件を用いて収束速度を導出する。特に、適切な定式化のためには r(S²) < 1 が必要である。
• 関数中心極限定理（CLT）をC(T)で用い、エントロピー条件およびサブガウス的または因子分解型モーメント条件に依存して信頼領域を構築する。
• インデックス集合Tの自然距離 d(t,s) = sup_x |K(t,x) - K(s,x)| / R(x) におけるメトリックエントロピーを理論的分析に活用する。R(x) = essup_t |K(t,x)| である。
• MarcusとJainのC(T)におけるCLTの結果を応用し、積分 I(t) = ∫g(t,x)ν(dx) に対して √n(I_n(t) - I(t)) の確率過程に適用する。

実験結果

リサーチクエスチョン

• RQ1 2種類目の線形Fredholm積分方程式を解くモンテカルロ法は、一様ノルムにおいて最適収束速度を達成できるか？
• RQ2 モンテカルロ推定器の関数中心極限定理が連続関数空間で成り立つための必要十分条件は何か？
• RQ3 このような積分方程式の解に対して、一様ノルムにおける信頼領域をどのように構築できるか？
• RQ4 核作用素の固有値半径がモンテカルロ法の収束性および安定性を保証するために果たす役割は何か？
• RQ5 核および定義域にどのような条件下で、自然距離 d(t,s) がエントロピーに基づくCLT条件を満たすほど良好に振る舞うのか？

主な発見

• 提案されたモンテカルロ法は、一様ノルムにおいて n^{-1/4} の最適収束速度を達成し、同じ仮定のもとで標準的な n^{-1/2} 速度よりも速い。
• 核作用素の二乗の固有値半径 r(S²) < 1 である限り、収束速度が最適であることが証明された。
• エントロピーおよびサブガウス条件を用いて、一様ノルムにおける非漸近的信頼領域が構築され、インデックス集合Tのメトリックエントロピーから明示的な境界が導出された。
• 翻訳不変距離の下でC(T)におけるCLTが成り立つために、|g(t,x) - g(s,x)| ≤ M(x)ρ(t,s) の因子分解条件が必要であることが確立された。
• 翻訳不変距離の下で、C(T)におけるCLTのための必要十分条件として、∫₀¹ H^{1/2}(T,ρ,z) dz < ∞ であることが示された。これは、やや弱いモーメントおよび正則性仮定のもとで成り立つ。
• ボルテラ方程式および準ノルム作用素に関する先行研究を一般化し、不連続または滑らかでない依存性を示す核のより広いクラスへと拡張した。