[論文レビュー] A Truncation Approach for Fast Computation of Distribution Functions

Xinjia Chen|ArXiv.org|Feb 23, 2008

Algorithms and Data Compression参考文献 2被引用数 5

ひとこと要約

本稿では、統合または和集合の定義域を計算的に効率的な小さな部分集合に制限することにより、多次元確率分布関数の計算を高速化するための切断手法を提案する。チェルノフ型の不等式と尾確率の制御を活用することで、切断された確率がユーザーが定めた誤差範囲内に保たれることを保証し、計算コストを著しく削減しながらも精度を維持する。

ABSTRACT

In this paper, we propose a general approach for improving the efficiency of computing distribution functions. The idea is to truncate the domain of summation or integration.

研究の動機と目的

• 大きな定義域上で多次元確率分布関数を評価する際の高い計算コストに対処すること。
• 精度を損なわせることなく統合／和集合の定義域を削減する一般的な手法を開発すること。
• 切断された確率が真の確率からの許容可能な誤差範囲内に保たれることを保証すること。
• モーメント母関数を用いて尾確率の境界を系統立てて計算する方法を提供すること。
• 高次元空間における和または積分に関する応用において、効率的かつ正確に分布関数を計算できるようにすること。

提案手法

• 新しい境界 $a_i' = \max(a_i, u_i)$ および $b_i' = \min(b_i, v_i)$ を定義することで、統合または和集合の定義域を切断する。ここで $u_i$ と $v_i$ は $\Pr\{X_i \leq u_i\} \leq \alpha_i$ および $\Pr\{X_i \geq v_i\} \leq \beta_i$ を満たすように選ばれる。
• 尾確率を境界付けるために、モーメント母関数 $\mathscr{C}(z) = \inf_{t \in \mathbb{R}} \mathbb{E}[e^{t(X-z)}]$ を用い、$\alpha_i, \beta_i$ が与えられたときの $u_i$ と $v_i$ を決定する。
• ボンフェローニの不等式を用いて誤差境界を導出し、$P \leq P' + \sum_{i=1}^m (\alpha_i + \beta_i)$ が成り立つことを示す。ここで $P'$ は切断された確率を表す。
• 所定の誤差許容誤差 $\eta$ に対して、$\alpha_i = \beta_i = \eta / (2m)$ と設定することで、$P' \leq P \leq P' + \eta$ を保証する。
• 集中不等式を用いて二項和への応用を行い、$\mathscr{M}(z, \mu)$ を用いて標本平均の尾確率を境界付ける。
• 誤差が $\eta$ 以内に保たれるようにするための切断閾値 $\mathrm{T}^-$ と $\mathrm{T}^+$ を定義する $z_1$ と $z_2$ の明示的表現を導出する。

実験結果

リサーチクエスチョン

• RQ1 多次元確率分布関数の計算複雑性を、精度を損なわせることなくどのように低減できるか？
• RQ2 多次元確率分布関数の定義域を切断する際の許容可能な最大誤差は何か？
• RQ3 モーメント母関数を用いて、どのように効率的に尾確率を境界付けることができるか？
• RQ4 ユーザーが定めた許容誤差内に真の確率に近い切断確率を系統立てて計算するための方法を開発できるか？
• RQ5 誤差が所定の $\eta$ 以内に保たれるようにするための切断閾値の明示的境界は何か？

主な発見

• 本手法は、切断確率 $P'$ が $P' \leq P \leq P' + \sum_{i=1}^m (\alpha_i + \beta_i)$ を満たすことを保証し、厳密な誤差境界を提供する。
• 所定の誤差許容誤差 $\eta > 0$ に対して、$\alpha_i = \beta_i = \eta / (2m)$ と設定することで、$P \leq P' + \eta$ が保証される。
• モーメント母関数 $\mathscr{C}(z)$ は尾確率に対してきめ細やかな境界を提供し、切断点 $u_i$ と $v_i$ の効率的決定を可能にする。
• $z_1$ と $z_2$ の明示的表現が $n$, $p$, および $\eta$ の関数として導出され、$\mathrm{T}^-$ と $\mathrm{T}^+$ の直接計算が可能になる。
• 計算の定義域が有限区間に削減され、和集合や積分における項の数が著しく減少する。
• 二項和への集中不等式の適用により、$\mathscr{M}(z, \mu)$ を用いて尾確率が指数関数的に減少することを実証した。