Iid随机变量之和的平方根的中心极限定理


11

出于对math.stackexchange的一个问题的兴趣,并进行了实证研究,我想知道以下有关iid随机变量之和的平方根的陈述。

假设是具有有限非零均值和方差 iid随机变量,并且。中心极限定理说随着增加。X1,X2,,Xnμσ2Y=i=1nXiÑYnμnσ2 d N(0,1)n

如果,我也可以说类似随着增加?ž - Z=|Y|ÑZn|μ|σ24|μ|σ24|μ| d N(0,1)n

例如,假设是具有均值和方差伯努利(Bernoulli,那么是二项式的,我可以在R中模拟它,例如: p p 1 p Y p = 1Xipp(1p)Yp=13

set.seed(1)
cases <- 100000
n <- 1000
p <- 1/3
Y <- rbinom(cases, size=n, prob=p)
Z <- sqrt(abs(Y))

给出了的期望均值和方差Z

> c(mean(Z), sqrt(n*p - (1-p)/4))
[1] 18.25229 18.25285
> c(var(Z), (1-p)/4)
[1] 0.1680012 0.1666667

和一个看起来像高斯的QQ图

qqnorm(Z)

在此处输入图片说明


1
@MichaelM:感谢您的评论。我从非负开始,但我认为您所描述的直观渐近行为允许将其推广到更多分布。令我惊讶的是:(a)和的平方根的方差明显趋于一个常数,而不取决于;(b)看起来非常接近高斯分布的外观。一个反例是可以接受的,但是当我尝试其他最初似乎不是高斯的情况时,进一步增加 似乎会使分布回到CLT类型的结果。 ñ ñXinn
亨利

这方面的一个推论是根均方(或二次平均)独立同分布的随机变量适当地缩放(乘以的与算术平均值)也收敛于规定,高斯分布的阶矩基础分布是有限的。 4n4
亨利

3
简而言之:该声明是Delta方法的特例,请参见Casella&Berger在“统计推断”一书中的定理5.5.24。
Michael M

@Michael:也许您看到的是我目前不知道的东西,但是我认为这个特殊问题不适合经典Delta方法的假设(例如,如您所参考的定理中所述)。请注意,不会在分布上收敛(非平凡地在),因此“应用的Delta方法”不能满足必要的要求。但是,正如S. Catterall的答案所证明的那样,它提供了一种有用的启发式方法,可以得出正确的答案。R g y = YRg(y)=|y|
红衣主教

(我相信您可以将Delta方法的证明适用于与上述相似的情况,以便充分严格地执行上述启发式方法。)
红衣主教

Answers:


14

收敛到高斯确实是普遍现象。

假设是均值和方差 IID随机变量,并定义和。固定数字。通常的中心极限定理告诉我们为,其中是标准普通CDF。但是,限制cdf的连续性意味着我们也有μ > 0 σ 2 Ŷ Ñ = Σ Ñ = 1 X α P ý ñ - ñ μX1,X2,X3,...μ>0σ2Yn=i=1nXiαñ→交通ΦPýñ-ñμP(Ynnμσnα)Φ(α)nΦPÝÑασ

P(Ynnμσnα+α2σ24μσn)Φ(α)
因为不等式右侧的附加项趋于零。重新排列此表达式会导致
Pÿñασ2μ+ñμ2Φα

取平方根,并注意到表示,我们得到换句话说,。该结果证明了在极限为收敛到高斯。P Y n < 0 0 P μ>0Pÿñ<00

P|ÿñ|ασ2μ+ñμΦα
|ÿñ|-ñμσ/2μdñ01个ñ

这是否意味着是大的良好近似值?好吧,我们可以做得更好。正如@Henry指出的那样,假设一切都是正数,我们可以使用以及和近似值,以获得改进的近似值如上述问题所述,。还要注意,我们仍然有因为ñμË[|ÿñ|]ñË[ÿñ]=Ë[ÿñ]-VarÿñË[ÿñ]=ñμVarÿñσ24μË[|ÿñ|]ñμ-σ24μ

|ÿñ|-ñμ-σ24μσ/2μdñ01个
ñμ-σ24μ-ñμ0为。ñ

您可能需要将作为以获得我的结果ñμ-ñμ-σ24μ0ñ
亨利

@Henry您可以将用替换为任何常数,这不会更改限制分布,但可能会更改对于特定的大很好地近似于。您是如何提出?ñμñμ-ķķ|ÿñ|-ñμ-ķσ/2μñ01个ññμ-σ24μ
S. Catterall恢复莫妮卡(Monica)

我们有所以。假设一切都是正的,则而建议,然后将这些引线组合为。Varž=Ë[ž2]-Ë[ž]2Ë[ž]=Ë[ž2]-VaržË[ž2]=Ë[ÿ]=ñμ|ÿñ|-ñμσ/2μ ë[Ž]听,说:Varžσ24μË[ž]ñμ-σ24μ
亨利

好的,谢谢,我已经尝试在我的回答中对此进行介绍。
S. Catterall恢复莫妮卡(Monica)
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.