[論文レビュー] A closed-form approximation for the median of the beta distribution

Jouni Kerman|arXiv (Cornell University)|Nov 2, 2011

Probabilistic and Robust Engineering Design被引用数 32

ひとこと要約

本稿では、$a, b > 1$ の場合に有効な、ベータ分布の中央値の簡単な閉形式近似を提案する：$(a - 1/3)/(a + b - 2/3)$。この近似の相対誤差は4%未満であり、形状母数が増加するに従い真の中央値に急速に収束する。他の $d$ 値を用いた $m(a,b;d) = (a-d)/(a+b-2d)$ の代替値と比較して、$d = 1/3$ がゼロ誤差への収束が最も速い。

ABSTRACT

A simple closed-form approximation for the median of the beta distribution Beta(a, b) is introduced: (a-1/3)/(a+b-2/3) for (a,b) both larger than 1 has a relative error of less than 4%, rapidly decreasing to zero as both shape parameters increase.

研究の動機と目的

• 一般解が存在しないベータ分布の中央値に対して、簡単で実用的かつ高精度な閉形式近似を構築すること。
• 不完全ベータ関数およびその逆関数に関する広範な文献にもかかわらず、中央値を計算するのに広く適用可能で計算効率の良い公式が欠如している問題に取り組むこと。
• 関数形 $m(a,b;d) = (a-d)/(a+b-2d)$ の $d$ の値を特定し、誤差を最小化し、対称性と真の中央値への収束を保証すること。
• 形状母数や分布の平均が変化する状況において、特に $a$ と $b$ が大きい極限において、近似の性能を評価すること。

提案手法

• 大規模な $a$ に対してガンマ分布の中央値が $a - 1/3$ であるという漸近的性質に基づき、近似 $m(a,b;1/3) = (a - 1/3)/(a + b - 2/3)$ を提案する。
• モード-中央値-平均の不等式を用いて、$m(a,b;d) = (a-d)/(a+b-2d)$ の関数形を正当化し、境界と対称性を保証する。
• 数値的に計算された中央値と比較して、$a$, $b$, $p = a/(a+b)$ のさまざまな値における相対誤差と絶対誤差を、図的および解析的手法を用いて評価する。
• ペイザー＝プリット近似を用いて、誤差の収束速度を分析し、$d = 1/3$ が $O(a^{-3/2})$ の誤差減衰を示し、他の $d$ 値よりも速いことを示す。
• 近似が中央値として適切であるかを確認するため、尾部確率 $\Pr(\theta \leq m(a,b;1/3))$ を評価し、確率が0.5に近いことを確認する。

実験結果

リサーチクエスチョン

• RQ1 ベータ分布の中央値に対して、簡潔で正確かつ計算効率の良い閉形式式を導出可能か？
• RQ2 関数形 $m(a,b;d) = (a-d)/(a+b-2d)$ の $d$ の選択が、近似誤差と収束速度にどのように影響するか？
• RQ3 近似 $m(a,b;1/3) = (a - 1/3)/(a + b - 2/3)$ は、相対誤差が4%未満であり、$a$ と $b$ が増加するに従い、真の中央値に一様に収束するか？
• RQ4 ペイザー＝プリットの漸近的解析が示唆するように、$d = 1/3$ は誤差減衰速度の観点で最適か？

主な発見

• すべての $a, b > 1$ に対して、近似 $(a - 1/3)/(a + b - 2/3)$ の相対誤差は4%未満であり、形状母数が増加するに従い急速に減少する。
• $a \geq 2$ の場合、相対誤差は1%未満に低下し、中程度の母数値に対しても高い精度を示す。
• $p < 0.5$ の場合に一貫して中央値を小さく見積もるが、$p > 0.5$ の場合は大きく見積もるが、$a$ と $b$ が増加するに従い、すべての $p$ に対して誤差は一様にゼロに収束する。
• 尾部確率 $\Pr(\theta \leq m(a,b;1/3))$ は、小さい方の形状母数が1以上であれば、$[0.4865, 0.5135]$ の区間内に保たれ、母数が増加するに従い0.5に急速に近づく。
• $d = 1/3$ の場合、誤差は $O(a^{-3/2})$ の速度で減少するが、$d = 0$（平均）などの他の $d$ 値では $O(a^{-1/2})$ の速度であるため、$d = 1/3$ が最適であることが確認される。
• 特に誤差減少の対数スケールにおいて、$d = 0.3$ よりも長期間にわたって誤差の低減と安定性に優れている。