MCPで接続 →

計算を入力してください

公式

広告

結果

Probability density f at x = 2
0.172252
ν = 3, λ = 1
x 確率密度 f
0 0
0.2 0.10121143
0.4 0.13381672
0.6 0.15315904
0.8 0.165206
1 0.17247566
1.2 0.1763617
1.4 0.17774925
1.6 0.17724876
1.8 0.17530452
2 0.17225201
2.2 0.16835122
2.4 0.16380739
2.6 0.15878474
2.8 0.15341592
3 0.14780871
3.2 0.14205106
3.4 0.13621485
3.6 0.13035878
3.8 0.12453071
4 0.11876944
4.2 0.11310625
4.4 0.10756608
4.6 0.10216859
4.8 0.09692892
5 0.09185846
5.2 0.08696543
5.4 0.08225536
5.6 0.07773156
5.8 0.07339542
6 0.06924682
6.2 0.06528429
6.4 0.06150532
6.6 0.05790652
6.8 0.05448379
7 0.05123249
7.2 0.04814754
7.4 0.04522352
7.6 0.04245479
7.8 0.03983554
8 0.03735987
8.2 0.03502185
8.4 0.03281554
8.6 0.03073508
8.8 0.02877465
9 0.02692857
9.2 0.02519127
9.4 0.02355732
9.6 0.02202148
9.8 0.02057864
10 0.01922389
10.2 0.0179525
10.4 0.01675992
10.6 0.01564179
10.8 0.01459393
11 0.01361235
11.2 0.01269324
11.4 0.01183297
11.6 0.01102809
11.8 0.01027531
12 0.00957151
12.2 0.00891374
12.4 0.0082992
12.6 0.00772523
12.8 0.00718933
13 0.00668912
13.2 0.00622237
13.4 0.00578697
13.6 0.00538093
13.8 0.00500237
14 0.00464952
14.2 0.00432073
14.4 0.00401444
14.6 0.00372916
14.8 0.00346354
15 0.00321626
15.2 0.00298612
15.4 0.00277198
15.6 0.00257277
15.8 0.00238748
16 0.00221518
16.2 0.00205499
16.4 0.0019061
16.6 0.00176772
16.8 0.00163914
17 0.0015197
17.2 0.00140876
17.4 0.00130573
17.6 0.00121007
17.8 0.00112127
18 0.00103884
18.2 0.00096235
18.4 0.00089137
18.6 0.00082552
18.8 0.00076445
19 0.0007078
19.2 0.00065527
19.4 0.00060657
19.6 0.00056142
19.8 0.00051957
20 0.00048079

非心カイ2乗分布とは

非心カイ2乗分布は、通常の(中心)カイ2乗分布に非心度パラメータ \(\lambda\) を加えて一般化した分布です。平均が 0 でない独立な正規確率変数の二乗和の分布を表し、検定の検出力分析(パワー分析)、信号検出、仮説検定など幅広い分野で用いられます。本計算はあくまで純粋な数学に基づくもので、特定の国・地域に依存するルールはなく、どこでも同じように利用できます。

非心度の増加に伴い右へずれていく非心カイ二乗密度曲線の一群
非心度λが大きくなると、非心カイ二乗分布の密度曲線は右へずれて平坦になります。

使い方

まず出力したい量を選びます。確率密度 \(f\)、下側累積確率 \(P\)、上側累積確率 \(Q\) のいずれかです。続いて自由度 \(\nu\)(0 より大きい値)、非心度 \(\lambda\)(0 以上の値)、基準となる \(x\) の値を入力します。さらに \(x\) の初期値、刻み幅、繰り返し回数(行数)を設定すると、指定した範囲にわたって(\(x\), 計算値)の表を一括で生成できます。

計算式の解説

非心カイ2乗分布は、ポアソン分布 Poisson(\(\lambda/2\)) を重みとした中心カイ2乗分布の混合分布として表せます。第 \(j\) 項の重みは \(w_j = e^{-\lambda/2}(\lambda/2)^j / j!\) です。確率密度は、自由度 \(\nu+2j\) の中心カイ2乗分布の密度に重み \(w_j\) を掛けて総和をとった $$f(x;\nu,\lambda)=\sum_{j=0}^{\infty}\frac{e^{-\lambda/2}\left(\lambda/2\right)^{j}}{j!}\,\frac{x^{\,k/2-1}e^{-x/2}}{2^{k/2}\,\Gamma\!\left(k/2\right)}\quad\text{where}\;\; k=\nu+2j,\;\; \lambda=\lambda$$ となります。下側累積確率は、同じ混合を中心カイ2乗分布の累積分布関数(CDF)に適用したもので、正規化された下側不完全ガンマ関数を用いて計算されます。 $$F(x;\nu,\lambda)=\sum_{j=0}^{\infty}\frac{e^{-\lambda/2}\left(\lambda/2\right)^{j}}{j!}\,P\!\left(\frac{\nu+2j}{2},\;\frac{x}{2}\right)\quad\text{where}\;\; \lambda=\lambda$$ 上側累積確率は単純に \(Q = 1 - P\) です。

広告
垂直線で左側の下側累積P領域と右側の上側累積Q領域に分けられた密度曲線
点xにおいて、下側累積Pは左側の面積、上側累積Qは右側の面積です。

計算例

\(\nu = 3\)、\(\lambda = 1\)、\(x = 2\) の場合を見てみましょう。Poisson(0.5) の重みは \(0.6065\)、\(0.3033\)、\(0.0758\)、\(0.0126\)、\(0.0016\) となります。\(x=2\) における自由度 3、5、7、9、11 の中心カイ2乗分布の CDF はそれぞれ \(0.4276\)、\(0.1511\)、\(0.0387\)、\(0.0074\)、\(0.0011\) です。これらの重み付き総和より \(P \approx 0.3082\)、したがって \(Q \approx 0.6918\) が得られます。同じ点での確率密度 \(f\) は約 \(0.173\) です。

よくある質問

\(\lambda = 0\) のときはどうなりますか? 重み 1 をもつ \(j=0\) の項だけが残るため、分布は自由度 \(\nu\) の中心カイ2乗分布そのものに一致します。

\(\nu\) は整数でなくてもよいですか? はい。ガンマ関数は 0 より大きい任意の \(\nu\) を扱えるため、小数(非整数)の自由度も有効です。

なぜ \(x = 0\) で密度が 0 になるのですか? \(\nu\) が 2 以上のとき原点での密度は 0 になります。\(\nu\) が 2 未満の場合は発散しますが、本計算では実用上の取り扱いとして \(x = 0\) では 0 を返します。

最終更新: