[論文レビュー] On Log-concavity of the Generalized Marcum Q Function

Yaming Yu|arXiv (Cornell University)|May 29, 2011

Advanced Combinatorial Mathematics参考文献 16被引用数 8

ひとこと要約

本稿は、Sun, Baricz, and Zhou (2010) が提起した予想を確認し、一般化されたMarcum Q関数 $ Q_{\nu}(a,b) $ が $ b \in [0,\infty) $ においてかつただかつそのときに限り $ \nu \geq 1/2 $ のとき対数凹型であることを証明する。証明は、修正ベッセル関数を含む被積分関数の対数凹型性を分析し、比 $ r_{\nu}(t) = I_{\nu}(t)/I_{\nu-1}(t) $ に基づく新規な十分条件を用いる。臨界パrameter値の数値的検証も行われている。

ABSTRACT

It is shown that, if nu >= 1/2 then the generalized Marcum Q function Q_nu(a, b) is log-concave in 0<=b

研究の動機と目的

• Sun, Baricz, and Zhou (2010) が提起した、一般化されたMarcum Q関数 $ Q_{\nu}(a,b) $ が $ b $ に関して対数凹型であるかどうかの予想を解決すること。
• 積分の被積分関数が $ t $ に関して対数凹型となる十分条件を確立することにより、積分の対数凹型性定理の適用を可能にすること。
• すべての $ a \geq 0 $ に対して被積分関数が $ t $ に関して対数凹型となる臨界閾値 $ \nu_0 \approx 0.7845 $ を特定すること。
• 結果を非心 $ \chi^2 $ 分布および信頼性理論と統計理論における生存関数の既知の性質と結びつけること。

提案手法

• 修正ベッセル関数およびその比の性質を用いて、$ Q_{\nu}(a,b) $ の被積分関数 $ f(t) $ の $ \log f(t) $ の2階微分を導出する。
• Lemma 1 を適用し、モードにおける連続性、単峰性、および減少部分における対数凹型性に基づく、生存関数の対数凹型性の十分条件を適用する。
• $ h_{\nu}(t) = 1 - \frac{2\nu-1}{t^2} - \frac{2\nu-1}{t}r_{\nu}(t) - r_{\nu}^2(t) $ を導入し、$ f(t) $ が $ t $ に関して対数凹型であるための必要十分条件が $ h_{\nu}(t) \leq 0 $ がすべての $ t > 0 $ で成り立つことであることを示す。
• 数値的検証により、方程式 $ \frac{I_{\nu}(\sqrt{5-2\nu})}{I_{\nu-1}(\sqrt{5-2\nu})} = \frac{3-2\nu}{\sqrt{5-2\nu}} $ の $ (1/2, 3/2) $ 内の唯一の解である $ \nu_0 \approx 0.78449776 $ の存在と値を特定する。
• $ h_{\nu}(t) $ および補助関数 $ l_{\nu}(t) $ の符号の変化と単調性を分析し、$ \nu \geq \nu_0 $ のとき $ h_{\nu}(t) $ が正の値を取らないことを排除することで、グローバルな対数凹型性を保証する。

実験結果

リサーチクエスチョン

• RQ1 一般化されたMarcum Q関数 $ Q_{\nu}(a,b) $ は、すべての $ \nu \geq 1/2 $ に対して $ b \in [0,\infty) $ で対数凹型であるか？
• RQ2 $ \nu $ と $ a $ がどのような条件下で、$ Q_{\nu}(a,b) $ の被積分関数が $ t \in (0,\infty) $ に関して対数凹型になるか？
• RQ3 すべての $ a \geq 0 $ に対して被積分関数が $ t $ に関して対数凹型となる最小の $ \nu $ は何か？
• RQ4 被積分関数の対数凹型性を用いて $ Q_{\nu}(a,b) $ の対数凹型性を確立できるか。このアプローチの限界は何か？
• RQ5 結果は、非心 $ \chi^2 $ 分布およびその生存関数の既知の対数凹型性の性質とどのように関連するか？

主な発見

• 一般化されたMarcum Q関数 $ Q_{\nu}(a,b) $ は、$ b \in [0,\infty) $ においてかつただかつそのときに限り $ \nu \geq 1/2 $ のとき対数凹型である。これはSun, Baricz, and Zhou (2010) の予想1を確認する。
• $ \nu \geq 1/2 $ のとき、$ Q_{\nu}(a,b) $ の被積分関数は $ t \in (0,\infty) $ に関して対数凹型であるための必要十分条件は $ 0 \leq a \leq 1 $ であり、これは命題1で示されている。
• すべての $ a \geq 0 $ に対して被積分関数が $ t $ に関して対数凹型であるための必要十分条件は $ \nu \geq \nu_0 $ であり、ここで $ \nu_0 \approx 0.78449776 $ は、$ (1/2, 3/2) $ 内の修正ベッセル関数の比を含む方程式の唯一の解である。
• 補題1により、$ 1 - Q_{\nu}(a,b) $ が $ b $ に関して対数凹型であることが示され、2つの条件が成り立つ：(i) $ \nu \geq 1/2 $ かつ $ 0 \leq a \leq 1 $、または (ii) $ \nu \geq \nu_0 $。
• 補題2により、$ \nu = 1 $ のとき、Rice分布の確率密度関数、累積分布関数、生存関数がすべて対数凹型であることが示され、$ Q_1(a,b) $ がRice分布に従う確率変数の生存関数に対応する。
• Finner and Roters (2001) の定理3は、本稿の結果と整合的であり、特に非心 $ \chi^2 $ 累積分布関数の対数凹型性に関して、部分的に包含されている。