選單圖示 主題
Search

伽瑪分佈


GammaDistribution

伽瑪分佈是一種通用的統計分佈型別,它與β分佈相關,並且自然地出現在泊松分佈事件之間等待時間相關的過程中。伽瑪分佈有兩個自由引數,標記為 alphatheta,上面展示了其中一些。

考慮分佈函式 D(x),表示在給定變化率為 lambda泊松分佈下,直到第 h 個泊松事件發生的等待時間:

D(x)=P(X<=x)
(1)
=1-P(X>x)
(2)
=1-sum_(k=0)^(h-1)((lambdax)^ke^(-lambdax))/(k!)
(3)
=1-e^(-lambdax)sum_(k=0)^(h-1)((lambdax)^k)/(k!)
(4)
=1-(Gamma(h,xlambda))/(Gamma(h))
(5)

對於 x in [0,infty),其中 Gamma(x) 是完全伽瑪函式,而 Gamma(a,x) 是不完全伽瑪函式。當 h 為整數時,此分佈是稱為厄爾朗分佈的特殊情況。

相應的機率函式 P(x),表示直到第 h 個泊松事件發生的等待時間的機率函式,透過對 D(x) 求導獲得:

P(x)=D^'(x)
(6)
=lambdae^(-lambdax)sum_(k=0)^(h-1)((lambdax)^k)/(k!)-e^(-lambdax)sum_(k=0)^(h-1)(k(lambdax)^(k-1)lambda)/(k!)
(7)
=lambdae^(-lambdax)+lambdae^(-lambdax)sum_(k=1)^(h-1)((lambdax)^k)/(k!)-e^(-lambdax)sum_(k=1)^(h-1)(k(lambdax)^(k-1)lambda)/(k!)
(8)
=lambdae^(-lambdax)-lambdae^(-lambdax)sum_(k=1)^(h-1)[(k(lambdax)^(k-1))/(k!)-((lambdax)^k)/(k!)]
(9)
=lambdae^(-lambdax){1-sum_(k=1)^(h-1)[((lambdax)^(k-1))/((k-1)!)-((lambdax)^k)/(k!)]}
(10)
=lambdae^(-lambdax){1-[1-((lambdax)^(h-1))/((h-1)!)]}
(11)
=(lambda(lambdax)^(h-1))/((h-1)!)e^(-lambdax).
(12)

現在令 alpha=h(不一定是整數),並定義 theta=1/lambda 為變化間隔時間。那麼上述方程可以寫成:

 P(x)=(x^(alpha-1)e^(-x/theta))/(Gamma(alpha)theta^alpha)
(13)

對於 x in [0,infty)。這是伽瑪分佈的機率函式,相應的分佈函式是:

 D(x)=P(alpha,x/theta),
(14)

其中 P(a,z) 是正則化伽瑪函式

它在Wolfram 語言中作為函式實現:GammaDistribution[alpha, theta].

描述此分佈的特徵函式是:

phi(t)=F_x{(x^(-x/theta)x^(alpha-1))/(Gamma(alpha)theta^alpha)[1/2(1+sgnx)]}(t)
(15)
=(1-ittheta)^(-alpha),
(16)

其中 F_x[f](t) 是引數為 a=b=1傅立葉變換,而矩量生成函式是:

M(t)=int_0^infty(e^(tx)x^(alpha-1)e^(-x/theta)dx)/(Gamma(alpha)theta^alpha)
(17)
=int_0^infty(x^(alpha-1)e^(-(1-thetat)x/theta)dx)/(Gamma(alpha)theta^alpha).
(18)

給出關於 0 的矩為:

 mu_r^'=(theta^rGamma(alpha+r))/(Gamma(alpha))
(19)

(Papoulis 1984,第 147 頁)。

為了明確找到使用矩量生成函式的分佈,令:

y=((1-thetat)x)/theta
(20)
dy=(1-thetat)/thetadx,
(21)

因此:

M(t)=int_0^infty((thetay)/(1-thetat))^(alpha-1)(e^(-y))/(Gamma(alpha)theta^alpha)(thetady)/(1-thetat)
(22)
=1/((1-thetat)^alphaGamma(alpha))int_0^inftyy^(alpha-1)e^(-y)dy
(23)
=1/((1-thetat)^alpha),
(24)

給出對數矩量生成函式為:

R(t)=-alphaln(1-thetat)
(25)
R^'(t)=(alphatheta)/(1-thetat)
(26)
R^('')(t)=(alphatheta^2)/((1-thetat)^2).
(27)

均值方差偏度超額峰度分別是:

mu=alphatheta
(28)
sigma^2=alphatheta^2
(29)
gamma_1=2/(sqrt(alpha))
(30)
gamma_2=6/alpha.
(31)

伽瑪分佈與其他統計分佈密切相關。如果 X_1X_2、...、X_n 是獨立的隨機變數,它們具有引數為 (alpha_1,theta)(alpha_2,theta)、...、(alpha_n,theta) 的伽瑪分佈,則 sum_(i=1)^(n)X_i 也服從伽瑪分佈,引數為:

alpha=sum_(i=1)^(n)alpha_i
(32)
theta=theta.
(33)

此外,如果 X_1X_2 是獨立的隨機變數,它們具有引數為 (alpha_1,theta)(alpha_2,theta) 的伽瑪分佈,則 X_1/(X_1+X_2) 是引數為 (alpha_1,alpha_2)β分佈變數。兩者都可以如下推匯出來。

 P(x_1,x_2)=1/(Gamma(alpha_1)Gamma(alpha_2))e^(x_1+x_2)x_1^(alpha_1-1)x_2^(alpha_2-1).
(34)

設:

 u=x_1+x_2    x_1=uv
(35)
 v=(x_1)/(x_1+x_2)    x_2=u(1-v),
(36)

雅可比行列式為:

 J((x_1,x_2)/(u,v))=|v u; 1-v -u|=-u,
(37)

因此:

 g(u,v)dudv=f(x,y)dxdy=f(x,y)ududv.
(38)
g(u,v)=u/(Gamma(alpha_1)Gamma(alpha_2))e^(-u)(uv)^(alpha_1-1)u^(alpha_2-1)(1-v)^(alpha_2-1)
(39)
=1/(Gamma(alpha_1)Gamma(alpha_2))e^(-u)u^(alpha_1+alpha_2-1)v^(alpha_1-1)(1-v)^(alpha_2-1).
(40)

因此,和 X_1+X_2 的分佈為:

 f(u)=f(x_1+x_2)=int_0^1g(u,v)dv=(e^(-u)u^(alpha_1+alpha_2-1))/(Gamma(alpha_1+alpha_2)),
(41)

這是一個伽瑪分佈,比率 X_1/(X_1+X_2) 的分佈為:

h(v)=h((x_1)/(x_1+x_2))
(42)
=int_0^inftyg(u,v)du
(43)
=(v^(alpha_1-1)(1-v)^(alpha_2-1))/(B(alpha_1,alpha_2)),
(44)

其中 Bβ函式,這是一個β分佈

如果 XY 是引數分別為 alpha_1alpha_2 的伽瑪變數,則 X/Y 是引數為 alpha_1alpha_2β' 分佈變數。設:

 u=x+y    v=x/y,
(45)

雅可比行列式為:

 J((u,v)/(x,y))=|1 1; 1/y -x/(y^2)|=-(x+y)/(y^2)=-((1+v)^2)/u,
(46)

因此:

 dxdy=u/((1+v)^2)dudv
(47)
g(u,v)=1/(Gamma(alpha_1)Gamma(alpha_2))e^(-u)((uv)/(1+v))^(alpha_1-1)(u/(1+v))^(alpha_2-1)u/((1+v)^2)
(48)
=1/(Gamma(alpha_1)Gamma(alpha_2))e^(-u)u^(alpha_1+alpha_2-1)v^(alpha_1-1)(1+v)^(-alpha_1-alpha_2).
(49)

因此,比率 X/Y 的分佈為:

 h(v)=int_0^inftyg(u,v)du=(v^(alpha_1-1)(1+v)^(-alpha_1-alpha_2))/(B(alpha_1,alpha_2)),
(50)

這是一個引數為 (alpha_1,alpha_2)β' 分佈

伽瑪分佈的“標準形式”透過令 y=x/theta 給出,因此 dy=dx/theta 且:

P(y)dy=(x^(alpha-1)e^(-x/theta))/(Gamma(alpha)theta^alpha)dx
(51)
=((thetay)^(alpha-1)e^(-y))/(Gamma(alpha)theta^alpha)(thetady)
(52)
=(y^(alpha-1)e^(-y))/(Gamma(alpha))dy,
(53)

因此,關於 0 的為:

nu_r=1/(Gamma(alpha))int_0^inftye^(-x)x^(alpha-1+r)dx
(54)
=(Gamma(alpha+r))/(Gamma(alpha))
(55)
=(alpha)_r,
(56)

其中 (alpha)_r波赫哈默爾符號。關於 mu=mu_1為:

mu_1=alpha
(57)
mu_2=alpha
(58)
mu_3=2alpha
(59)
mu_4=3alpha^2+6alpha.
(60)

矩量生成函式為:

 M(t)=1/((1-t)^alpha),
(61)

累積量生成函式為:

 K(t)=alphaln(1-t)=alpha(t+1/2t^2+1/3t^3+...),
(62)

因此累積量為:

 kappa_r=alphaGamma(r).
(63)

如果 X正態變數,均值mu標準差sigma,則:

 Y=((X-mu)^2)/(2sigma^2)
(64)

是引數為 alpha=1/2 的標準伽瑪變數。


另請參閱

β分佈, χ² 分佈, 厄爾朗分佈

使用 探索

參考文獻

Beyer, W. H. CRC Standard Mathematical Tables, 28th ed. Boca Raton, FL: CRC Press, p. 534, 1987.Jambunathan, M. V. "Some Properties of Beta and Gamma Distributions." Ann. Math. Stat. 25, 401-405, 1954.Papoulis, A. Probability, Random Variables, and Stochastic Processes, 2nd ed. New York: McGraw-Hill, pp. 103-104, 1984.

在 中被引用

伽瑪分佈

引用為

Weisstein, Eric W. “伽瑪分佈。” 來自 網路資源。 https://mathworld.tw/GammaDistribution.html

主題分類