Poisson分布に従う少ない観測量に対する誤差統計量の近似式
Poisson 分布は観測値が大きければ中心極限定理により Guassian 分布に近づくわけですが,得てして天文業界では観測値が大きくない場合が多いです.そんなわけで Poisson 分布関数をそれなりにまじめに計算しないと本当はいけない場合が多いです.ざっと の場合に Poisson 分布的な振る舞いになります.その場合,以下のような近似式を使えば まで 1% 以下の精度で計算できます.ただし, は でゼロ割をしてしまうので三項間演算子を使うなどして逃げる必要があります. のように.
Confidence limits for small numbers of events in astrophysical data by Gehrels
http://ads.nao.ac.jp/abs/1986ApJ...303..336G
upper limit:
lower limit:
Parameters | CL=0.8413 | 0.9772 | 0.9987 |
---|---|---|---|
S | 1.000 | 2.000 | 3.000 |
0.0 | 0.062 | 0.222 | |
... | -2.19 | -1.88 |
S は gaussian の らしい.
精度は1%以下
参考
http://d.hatena.ne.jp/youtsumi/20090119/1232290800