ここでは、[離散分布のあてはめ]メニューで用意されている統計分布について説明します。
Poisson分布では、尺度パラメータλ > 0が推定されます。
確率関数: ただし 0 ≤ λ < ∞; x = 0,1,2,...
E(x) = λ
Var(x) = λ
Poisson分布は離散分布なので、分布関数は整数ごとに変化する階段関数になります。
負の二項分布は、特定の失敗回数に達するまでの成功回数をモデル化するのに便利です。以下のパラメータ表現では、平均パラメータλおよび過分散パラメータσが含まれています。
確率関数:
E(x) = λ
Var(x) = λ + σλ2
ここで、Γ(·)は、ガンマ関数です。
負の二項分布はガンマPoisson分布と同等です。ガンマPoisson分布は、平均μが異なる複数のPoisson分布から、データが生成されているときに役立ちます。
ガンマPoisson分布は、x|μが平均μのPoisson分布に従い、その平均がパラメータ(α,τ)のガンマ分布に従うという仮定から導出できます。ガンマPoisson分布には、λ = ατとσ = τ+1の2つのパラメータがあります。パラメータσは過分散パラメータです。σ > 1の場合は、過大分散、つまり通常のPoisson分布よりも分散が大きくなります。σ = 1の場合、xの分布は平均λのPoisson分布になります。
確率関数: ただし 0 < λ; 1 ≤ σ; x = 0,1,2,...
E(x) = λ
Var(x) = λσ
ここで、Γ(·)は、ガンマ関数です。
ガンマPoisson分布は、σnegbin = (σgp - 1) / λgpの負の二項分布と同等です。
「Samples/Scripts」フォルダの中の「demoGammaPoisson.jsl」を実行すると、パラメータλとσのガンマPoisson分布とパラメータλのPoisson分布を比較できます。
ゼロ強調 Poisson分布には、尺度パラメータ(λ > 0)と、ゼロ強調パラメータ(π)があります。
確率関数:
E(x) = (1 - π)λ
Var(x) = λ(1 - π)(1 + λπ)
ゼロ強調 負の二項分布には、尺度パラメータ(λ > 0)、過分散パラメータ(σ > 0)、ゼロ強調パラメータ(π)があります。
確率関数:
E(x) = (1 - π)λ
Var(x) = λ(1 - π)[1 + λ(σ + π)]
[二項のあてはめ]オプションでは、一定の標本サイズか、標本サイズを含む列のいずれかを指定することができます。
確率関数: ただし 0 ≤ p ≤ 1; x = 0,1,2,...,n
E(x) = np
Var(x) = np(1-p)
上の式で、nは、独立した試行の回数を示します。
メモ: 二項分布のパラメータに対する信頼区間の計算には、スコア信頼区間が使用されています。Agresti and Coull(1998)を参照してください。
ベータ二項分布は、確率πが異なる複数の二項分布から、データが生成されているときに役立ちます。たとえば、複数の製造ラインにおける不適合品の個数などが挙げられます。この場合、不適合となる確率(π)は、製造ラインによって異なると考えられます。
ベータ二項分布は、x|πが二項分布(n,π)に従い、πがベータ分布(α,β)に従うという仮定から導出できます。ベータ二項分布には、p = α/(α+β)とδ = 1/(α+β+1)の2つのパラメータがあります。パラメータδは、過分散パラメータです。δ > 0の場合は、過大分散、つまり通常の二項分布よりも分散が大きくなります。δ < 0の場合は、過小分散です。δ = 0の場合、xは二項分布(n,p)に従います。ベータ二項分布は、n ≥ 2の場合にのみ、あてはめることができます。
確率関数:
ただし ; ; x = 0,1,2,...,n
E(x) = np
Var(x) = np(1-p)[1+(n-1)δ]
ここで、Γ(·)は、ガンマ関数です。
前述のとおり、x|πは二項分布(n,π)に従い、πはベータ分布(α,β)に従うと仮定されています。このプラットフォームでは、p = α/(α+β)とδ = 1/(α+β+1)のパラメータ推定値が計算されます。αとβの推定値は、次式で計算できます。
ただし、δの推定値が0の場合は、上記の式が使用できません。その場合、ベータ二項分布は、確率pの二項分布となり、がそのpの推定値となります。
なお、ベータ二項分布におけるパラメータの信頼区間は、プロファイル尤度法によって計算されます。
「Samples/Scripts」フォルダにある「demoBetaBinomial.jsl」を実行すると、nを20に固定しパラメータpを可変として、二項分布と過分散パラメータδのベータ二項分布を比較できます。