[論文レビュー] A note on the distribution of the maximum of a set of Poisson random variables

Keith Briggs, Li Song+1 more|ArXiv.org|Mar 25, 2009

Mathematical functions and polynomials参考文献 4被引用数 12

ひとこと要約

この論文は、n 個の i.i.d. ポisson(λ) 確率変数の最大値のモードに対する漸近的近似を精緻化し、ラムゼット・W関数とニュートン・ラプソンの修正を用いたより正確な公式を導入している。主な貢献は、テストされたすべての n および λ に対して誤差が1未塔である非常に高精度の最大値の最も確率の高い値の推定であり、キムバー（1983）の先行研究を著しく上回っている。

ABSTRACT

Given a set of independent Poisson random variables with common mean, we study the distribution of their maximum and obtain an accurate asymptotic formula to locate the most probable value of the maximum. We verify our analytic results with very precise numerical computations.

研究の動機と目的

• n 個の i.i.d. ポアソン(λ) 確率変数の最大値のモードに対するキムバー（1983）の漸近的近似を改善すること。
• 特に大きな n に対して、最大値の最も確率の高い値のより正確で解析的に取り扱いやすい公式を導出すること。
• 大きな n に対して、最大値が高々2つの隣接整数に集中する「集約効果」を特定すること。
• gλ(x) = 1/n の根に対する高精度な漸近的展開を提供すること。ここで gλ(x) はポアソン最大値の累積分布関数である。
• 広範な n および λ の範囲で高精度な数値計算を用いて、改良された近似の妥当性を検証すること。

提案手法

• n 個のポアソン(λ) 変数の最大値の CDF を近似する連続的補間関数 gλ(x) = 1 − Γ(x+1, λ)/Γ(x+1) を定義する。
• log(gλ(x)) の x, λ, n に関する漸近的展開を、O(x⁻²) 項まで保持して行い、x の対数および逆数の項を含む。
• 主分枝を用いた一次近似 x₀ = log n / W(log n / (eλ)) を導出する。
• log(gλ(x)) 展開の完全な微分を用いて、x₀ に対して1回のニュートン・ラプソン補正ステップを適用し、精錬された推定値 x₁ を得る。
• 構成上、精錬された推定値 x₁ は誤差が1未塔であることを保証し、モードの正しく特定可能である。
• n = 10⁰ から 10⁴⁰ および λ = 0.5, 1, 2, 5 の範囲で、x₀ および x₁ の精度を数値的に検証し、x₁ が正確なモードの2つの整数ステップの間に常に存在することを確認した。

実験結果

リサーチクエスチョン

• RQ1 n 個の i.i.d. ポアソン(λ) 確率変数の最大値のモードを、既存の漸近的公式を上回る精度でどのように近似できるか？
• RQ2 n → ∞ のとき、最大値の最も確率の高い値の漸近的挙動は何か？ そしてそれは λ にどのように依存するか？
• RQ3 gλ(x) = 1/n の根に対する精錬された漸近的展開は、すべての n および λ に対して誤差が1未塔の推定値をもたらすか？
• RQ4 キムバーが以前に指摘した最大値の確率質量の振動的挙動は、精錬された近似でどのように現れるか？
• RQ5 新しい公式は、収束性および精度の観点から、キムバーの近似をどの程度上回るか？

主な発見

• ラムゼット・W関数に基づく近似 x₀ にニュートン・ラプソン補正を施した精錬された漸近的推定値 x₁ は、テストされたすべての n および λ に対して誤差が1未塔である。
• 公式 x₁ = x₀ + [log λ − λ − log(2π)/2 − (3/2)log x₀] / [log x₀ − log λ] は、最大ポアソン変数のモードを非常に正確に近似する。
• 数値結果から、x₁ が正確なモードの2つの整数ステップの間に常に存在することが示され、その精度が確認され、モードの正しく特定可能であることが裏付けられた。
• ラムゼット・W関数を用いた近似 x₀ は、キムバーの log n / log log n 公式に比べて著しく改善されており、特に λ 依存性を的確に捉えている。
• 最大値がその2つの最も確率の高い値のいずれかをとる確率 Pn = Pr[Mn ∈ {In, In+1}] は、n → ∞ のときも持続的に振動する。この挙動は論文の図で確認され、可視化された。
• 最大値が高々2つの隣接整数に集中する「集約効果」は、数値的に確認され、高精度な漸近的近似の可能性を裏付けるものである。