欢迎访问 生活随笔!

生活随笔

当前位置: 首页 > 编程资源 > 编程问答 >内容正文

编程问答

UA MATH564 概率论VI 数理统计基础3 卡方分布中

发布时间:2025/4/14 编程问答 32 豆豆
生活随笔 收集整理的这篇文章主要介绍了 UA MATH564 概率论VI 数理统计基础3 卡方分布中 小编觉得挺不错的,现在分享给大家,帮大家做个参考.

UA MATH564 概率论VI 数理统计基础3 卡方分布中

  • 卡方分布的基本性质

上一讲介绍了卡方分布的定义:假设X1,⋯,XnX_1,\cdots,X_nX1,,Xn互相独立,并且Xi∼N(ai,1)X_i \sim N(a_i,1)XiN(ai,1),则称
∑i=1nXi2∼χ2(n,δ)\sum_{i=1}^n X_i^2 \sim \chi^2(n,\delta)i=1nXi2χ2(n,δ)
其中nnn代表样本数,δ\deltaδ是非中心化参数
δ=∑i=1nai2\delta = \sqrt{\sum_{i=1}^n a_i^2}δ=i=1nai2
如果样本来自标准正态总体,则δ=0\delta=0δ=0,称之为中心化的卡方分布。并推导了它的概率密度为:
k(x∣n,δ)=e−δ2+x2∑i=0∞δ2i2ii!xi+n/2−12i+n/2Γ(i+n/2)k(x|n,\delta) = e^{-\frac{\delta^2+x}{2}}\sum_{i=0}^{\infty}\frac{\delta^{2i}}{2^ii!} \frac{x^{i+n/2-1}}{2^{i+n/2}\Gamma(i+n/2)}k(xn,δ)=e2δ2+xi=02ii!δ2i2i+n/2Γ(i+n/2)xi+n/21
δ=0\delta=0δ=0时,
k(x∣n,0)=K′(x∣n,0)=(1/2)n/2Γ(n/2)xn2−1e−x/2=dΓ(n2,12)k(x|n,0) = K'(x|n,0) = \frac{(1/2)^{n/2}}{\Gamma(n/2)}x^{\frac{n}{2}-1}e^{-x/2}=_d \Gamma(\frac{n}{2},\frac{1}{2})k(xn,0)=K(xn,0)=Γ(n/2)(1/2)n/2x2n1ex/2=dΓ(2n,21)
下面介绍卡方分布的一些常用性质:

卡方分布的基本性质

  • X∼χ2(n)X\sim \chi^2(n)Xχ2(n),则EX=n,Var(X)=2nEX=n,\ Var(X)=2nEX=n, Var(X)=2n
  • Xn∼χ2(n),n=1,2,⋯X_n\sim \chi^2(n),n=1,2,\cdotsXnχ2(n),n=1,2,,则n→∞n \to \inftynXn−n2n→dN(0,1),2Xn−2n∼dN(0,1)\frac{X_n-n}{\sqrt{2n}} \to_d N(0,1),\ \ \sqrt{2X_n}-\sqrt{2n} \sim_d N(0,1)2nXnndN(0,1),  2Xn2ndN(0,1)
  • Y1,⋯,YmY_1,\cdots,Y_mY1,,Ym互相独立,Yi∼χni,δi2Y_{i}\sim \chi^2_{n_i,\delta_i}Yiχni,δi2,则∑i=1mYi∼χn,δ2,n=∑i=1mni,δ2=∑i=1nδi2\sum_{i=1}^m Y_i \sim \chi^2_{n,\delta},\ n = \sum_{i=1}^m n_i, \delta^2 = \sum_{i=1}^n \delta_i^2i=1mYiχn,δ2, n=i=1mni,δ2=i=1nδi2
  • 用一般的正态分布构造卡方分布的方法在数理统计基础1中已经介绍过了,这里不再重复,性质1可以用gamma分布的规律直接得到,性质2的前半部分就是Classical CLT,性质3描述的是卡方分布的可加性。下面给出性质2的后半部分和性质3的证明。
    证明
    性质2第二个式子,考虑
    P(2Xn−2n≤x)=P(Xn≤x2+22nx+2n2)=P(Xn−n2n≤x+x222n)P(\sqrt{2X_n}-\sqrt{2n} \le x) = P(X_n \le \frac{x^2 + 2\sqrt{2n}x+2n}{2}) = P(\frac{X_n-n}{\sqrt{2n}}\le x + \frac{x^2}{2\sqrt{2n}})P(2Xn2nx)=P(Xn2x2+22nx+2n)=P(2nXnnx+22nx2)
    n→∞n\to \inftyn时,Xn−n2n→dN(0,1)\frac{X_n-n}{\sqrt{2n}} \to_d N(0,1)2nXnndN(0,1)x222n→0\frac{x^2}{2\sqrt{2n}} \to 022nx20,因此2Xn−2n→dN(0,1)\sqrt{2X_n}-\sqrt{2n} \to_d N(0,1)2Xn2ndN(0,1)

    性质3,中心化的卡方分布的可加性根据gamma分布的可加性可以直接得到,下面考虑非中心化的卡方分布的可加性,两两可加推广到有限可加是非常平凡的,这里证明两两可加的情况:因为Y1,Y2Y_1,Y_2Y1,Y2是互相独立的,因此
    fY1+Y2(y)=fY1∗fY2(y)=∫z=0yfY1(y−z)fY2(z)dzf_{Y_1+Y_2}(y) = f_{Y_1}*f_{Y_2}(y) = \int_{z=0}^{y}f_{Y_1}(y-z)f_{Y_2}(z)dzfY1+Y2(y)=fY1fY2(y)=z=0yfY1(yz)fY2(z)dz

    经过上一讲对非中心化卡方分布的介绍,相信大家对它的密度函数已经驾轻就熟了,我们把被积函数写出来
    fY1(y−z)fY2(z)=(e−δ12+y−z2∑i=0∞δ12i2ii!(y−z)i+n1/2−12i+n1/2Γ(i+n1/2))(e−δ22+z2∑j=0∞δ22j2jj!zj+n2/2−12j+n2/2Γ(j+n2/2))f_{Y_1}(y-z)f_{Y_2}(z) = \left( e^{-\frac{\delta_1^2+y-z}{2}}\sum_{i=0}^{\infty}\frac{\delta_1^{2i}}{2^ii!} \frac{(y-z)^{i+n_1/2-1}}{2^{i+n_1/2}\Gamma(i+n_1/2)} \right) \left( e^{-\frac{\delta_2^2+z}{2}}\sum_{j=0}^{\infty}\frac{\delta_2^{2j}}{2^jj!} \frac{z^{j+n_2/2-1}}{2^{j+n_2/2}\Gamma(j+n_2/2)} \right)fY1(yz)fY2(z)=(e2δ12+yzi=02ii!δ12i2i+n1/2Γ(i+n1/2)(yz)i+n1/21)(e2δ22+zj=02jj!δ22j2j+n2/2Γ(j+n2/2)zj+n2/21)
    虽然看上去很吓人但不要慌,需要的技巧我们在上一讲已经学会了——构造beta函数求积。先把上面的式子展开,
    RHS=e−(δ12+δ22)+y2∑i,j=0∞δ12iδ22j2i+ji!j!(y−z)i+n1/2−1zj+n2/2−12i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)fY1+Y2(y)=e−(δ12+δ22)+y2∑i,j=0∞δ12iδ22j2i+ji!j!∫0y(y−z)i+n1/2−1zj+n2/2−1dz2i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)RHS = e^{-\frac{(\delta_1^2+ \delta^2_2)+y}{2}}\sum_{i,j=0}^{\infty} \frac{\delta_1^{2i}\delta_2^{2j}}{2^{i+j}i!j!} \frac{(y-z)^{i+n_1/2-1}z^{j+n_2/2-1}}{2^{i+j+(n_1+n_2)/2}\Gamma(i+n_1/2)\Gamma(j+n_2/2)} \\ f_{Y_1+Y_2}(y) = e^{-\frac{(\delta_1^2+ \delta^2_2)+y}{2}}\sum_{i,j=0}^{\infty} \frac{\delta_1^{2i}\delta_2^{2j}}{2^{i+j}i!j!} \frac{\int_{0}^y(y-z)^{i+n_1/2-1}z^{j+n_2/2-1}dz}{2^{i+j+(n_1+n_2)/2}\Gamma(i+n_1/2)\Gamma(j+n_2/2)}RHS=e2(δ12+δ22)+yi,j=02i+ji!j!δ12iδ22j2i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)(yz)i+n1/21zj+n2/21fY1+Y2(y)=e2(δ12+δ22)+yi,j=02i+ji!j!δ12iδ22j2i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)0y(yz)i+n1/21zj+n2/21dz

    t=z/yt=z/yt=z/y,则
    ∫0y(y−z)i+n1/2−1zj+n2/2−1dz=yi+j+(n1+n2)/2−1∫01(1−t)i+n1/2−1tj+n2/2−1dz=yi+j+(n1+n2)/2−1B(i+n1/2,j+n2/2)=yi+j+(n1+n2)/2−1Γ(i+n1/2)Γ(j+n2/3)Γ(i+j+(n1+n2)/2)\int_{0}^y(y-z)^{i+n_1/2-1}z^{j+n_2/2-1}dz = y^{i+j+(n_1+n_2)/2-1}\int_{0}^1 (1-t)^{i+n_1/2-1}t^{j+n_2/2-1}dz \\ = y^{i+j+(n_1+n_2)/2-1}B(i+n_1/2,j+n_2/2) =y^{i+j+(n_1+n_2)/2-1}\frac{\Gamma(i+n_1/2)\Gamma(j+n_2/3)}{\Gamma(i+j+(n_1+n_2)/2)}0y(yz)i+n1/21zj+n2/21dz=yi+j+(n1+n2)/2101(1t)i+n1/21tj+n2/21dz=yi+j+(n1+n2)/21B(i+n1/2,j+n2/2)=yi+j+(n1+n2)/21Γ(i+j+(n1+n2)/2)Γ(i+n1/2)Γ(j+n2/3)

    带回到fY1+Y2(y)f_{Y_1+Y_2}(y)fY1+Y2(y)的表达式中,
    fY1+Y2(y)=e−(δ12+δ22)+y2∑i,j=0∞δ12iδ22j2i+ji!j!∫0y(y−z)i+n1/2−1zj+n2/2−1dz2i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)=e−(δ12+δ22)+y2∑i,j=0∞δ12iδ22j2i+ji!j!yi+j+(n1+n2)/2−12i+j+(n1+n2)/2Γ(i+j+(n1+n2)/2)f_{Y_1+Y_2}(y) = e^{-\frac{(\delta_1^2+ \delta^2_2)+y}{2}}\sum_{i,j=0}^{\infty} \frac{\delta_1^{2i}\delta_2^{2j}}{2^{i+j}i!j!} \frac{\int_{0}^y(y-z)^{i+n_1/2-1}z^{j+n_2/2-1}dz}{2^{i+j+(n_1+n_2)/2}\Gamma(i+n_1/2)\Gamma(j+n_2/2)} \\ = e^{-\frac{(\delta_1^2+ \delta^2_2)+y}{2}}\sum_{i,j=0}^{\infty} \frac{\delta_1^{2i}\delta_2^{2j}}{2^{i+j}i!j!} \frac{y^{i+j+(n_1+n_2)/2-1}}{2^{i+j+(n_1+n_2)/2}\Gamma(i+j+(n_1+n_2)/2)}fY1+Y2(y)=e2(δ12+δ22)+yi,j=02i+ji!j!δ12iδ22j2i+j+(n1+n2)/2Γ(i+n1/2)Γ(j+n2/2)0y(yz)i+n1/21zj+n2/21dz=e2(δ12+δ22)+yi,j=02i+ji!j!δ12iδ22j2i+j+(n1+n2)/2Γ(i+j+(n1+n2)/2)yi+j+(n1+n2)/21

    接下来重新安排一下求和的指标,因为{(i,j):i,j=0,1,⋯}\{(i,j):i,j=0,1,\cdots\}{(i,j):i,j=0,1,}的势和{k=0,1,⋯}\{k=0,1,\cdots\}{k=0,1,}是一样的,不妨记i+j=ki+j=ki+j=k,而δ12iδ22ji!j!\frac{\delta_1^{2i}\delta_2^{2j}}{i! j!}i!j!δ12iδ22j正好是多项式(δ12+δ22)k(\delta_1^2 + \delta_2^2)^k(δ12+δ22)k展开的第i,ji,ji,j项,定义
    n=∑i=1mni,δ2=∑i=1nδi2n = \sum_{i=1}^m n_i, \delta^2 = \sum_{i=1}^n \delta_i^2n=i=1mni,δ2=i=1nδi2

    上式可以化简为
    fY1+Y2(y)=e−(δ12+δ22)+y2∑k=0∞(δ12+δ22)k2kk!yk+(n1+n2)/2−12k+(n1+n2)/2Γ(k+(n1+n2)/2)=e−δ2+y2∑k=0∞δ2k2kk!yk+n/2−12k+n/2Γ(k+n/2)f_{Y_1+Y_2}(y) = e^{-\frac{(\delta_1^2+ \delta^2_2)+y}{2}}\sum_{k=0}^{\infty} \frac{(\delta_1^2 + \delta_2^2)^k}{2^{k}k!} \frac{y^{k+(n_1+n_2)/2-1}}{2^{k+(n_1+n_2)/2}\Gamma(k+(n_1+n_2)/2)} \\ = e^{-\frac{\delta^2+y}{2}}\sum_{k=0}^{\infty} \frac{\delta^{2k}}{2^{k}k!} \frac{y^{k+n/2-1}}{2^{k+n/2}\Gamma(k+n/2)}fY1+Y2(y)=e2(δ12+δ22)+yk=02kk!(δ12+δ22)k2k+(n1+n2)/2Γ(k+(n1+n2)/2)yk+(n1+n2)/21=e2δ2+yk=02kk!δ2k2k+n/2Γ(k+n/2)yk+n/21

    显然Y1+Y2∼χn,δ2Y_1 + Y_2 \sim \chi^2_{n,\delta}Y1+Y2χn,δ2

    总结

    以上是生活随笔为你收集整理的UA MATH564 概率论VI 数理统计基础3 卡方分布中的全部内容,希望文章能够帮你解决所遇到的问题。

    如果觉得生活随笔网站内容还不错,欢迎将生活随笔推荐给好友。