期望值1 / x1/x什么时候Xx遵循 Beta 分布

机器算法验证 可能性 期望值 贝塔分布 力矩生成函数
2022-03-18 04:27:48

x有概率密度:

f(x)=xα1(1x)1βB(α,β)

在哪里α,β是两个正参数和0x1是域x. 的期望值是多少1/x? 也就是说,什么是价值:

1/x=01f(x)xdx

2个回答

首先请注意,Beta 的 pdf(α,β)分布仅定义为α,β>0. 这意味着当α0或者β0

01xα1(1x)1βB(α,β)
不是有限的。现在,
E(1X)=011xxα1(1x)1βB(α,β)dx=01x(α1)1(1x)1βB(α1,β)Beta(α1,β)B(α1,β)B(α,β)dx=B(α1,β)B(α,β) if α1 > 0.

因此当α>1, 期望E(1/X)是有限的并且如上所述(这可以根据 Nate Pope 的评论进一步简化)。否则,它是未定义的。

我想指出另一种有趣的解决方法,它也将结果推广到期望Xm对于整数m=1,2,3,. 我将使用矩生成函数 (mgf) 和 N Cressie 等人的论文中的结果http://amstat.tandfonline.com/doi/abs/10.1080/00031305.1981.10479334?journalCode=utas20 “矩生成函数和负整数矩”。

他们给出的结果是,当X是一个随机变量,mgfMX(t)它定义在原点的开放邻域中,那么我们有

EXm=Γ(m)10tm1MX(t)dt
对于正整数m.

众所周知,对于 beta 分布,mgf 由一个汇合的超几何函数给出

MX(t)=1F1(α;α+β;t)
所以使用上面的结果给出了
EXm=Γ(m)10tm11F1(α;α+β;t)dt
我在 maple 的帮助下评估了这个积分:

assume( a>0, b>0 );assume(m-1,posint)
GAMMA(m)^(-1) * int( t^(m-1)*hypergeom([a],[a+b],-t), t=0..infinity )
                   GAMMA(a + b) GAMMA(a - m)
                   -------------------------
                   GAMMA(a) GAMMA(a + b - m)

所以最后我们可以把结果写成

EXm=Γ(α+β)Γ(αm)Γ(α)Γ(α+βm)
与其他答案一致m=1. 然后需要一些人类数学来得出结论,我们需要假设α>m为有效。