菜单图标 主题
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 的标准伽玛变量。


另请参阅

β分布, χ² 分布, 厄尔朗分布

使用 Wolfram|Alpha 探索

参考文献

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.

在 Wolfram|Alpha 中被引用

伽玛分布

引用为

Weisstein, Eric W. “伽玛分布。” 来自 MathWorld--Wolfram 网络资源。 https://mathworld.net.cn/GammaDistribution.html

主题分类