UA MATH564 概率论 高阶矩的计算:有限差算子方法2
UA MATH564 概率論 高階矩的計算:有限差算子方法2
- 二項分布的高階矩
- Poisson分布的高階矩
- 負二項分布的高階矩
這一講介紹有限差算子方法怎么用來計算常見的離散型隨機變量的高階矩。這里提到的高階矩特指高階原點矩,原點矩與中心矩之間可以互相轉換:
E[X?μ]r=E(∑j=0rCrjXjμr?j)=∑j=0rCrjμr?jEXjE[X-\mu]^r = E \left( \sum_{j=0}^r C_r^j X^j \mu^{r-j} \right) = \sum_{j=0}^r C_r^j \mu^{r-j} EX^jE[X?μ]r=E(j=0∑r?Crj?Xjμr?j)=j=0∑r?Crj?μr?jEXj
因為高階原點矩都是用第二類Stirling數表示的,這里給一個S2(n,k)S_2(n,k)S2?(n,k)的表
后續的數值可以用遞推關系計算:
S2(n+1,k)=kS2(n,k)+S2(n,k?1)S_2(n+1,k) = kS_2(n,k) + S_2(n,k-1)S2?(n+1,k)=kS2?(n,k)+S2?(n,k?1)
二項分布的高階矩
二項分布X~Binom(n,p)X \sim Binom(n,p)X~Binom(n,p)的mass function為
fX(x)=Cnxpx(1?p)n?x,x=0,1,?,nf_X(x) = C_n^x p^x (1-p)^{n-x},x=0,1,\cdots,nfX?(x)=Cnx?px(1?p)n?x,x=0,1,?,n
下面計算它的rrr階原點矩
EXr=∑x=0nxrCnxpx(1?p)n?x=∑x=0nCnxpx(1?p)n?xMx0r=(pM+q)n0r=(1+pΔ)n0r=∑x=0rCnxpxΔx0r=∑x=0rCnxpxS2(r,x)x!=∑x=0rAnxpxS2(r,x)EX^r = \sum_{x=0}^{n} x^rC_n^x p^x (1-p)^{n-x} = \sum_{x=0}^{n} C_n^x p^x (1-p)^{n-x}M^x0^r \\ = (pM+q)^n 0^r = (1+p\Delta)^n 0^r = \sum_{x=0}^r C_n^xp^x \Delta^x 0^r \\= \sum_{x=0}^r C_n^x p^x S_2(r,x)x! = \sum_{x=0}^r A_n^x p^x S_2(r,x)EXr=x=0∑n?xrCnx?px(1?p)n?x=x=0∑n?Cnx?px(1?p)n?xMx0r=(pM+q)n0r=(1+pΔ)n0r=x=0∑r?Cnx?pxΔx0r=x=0∑r?Cnx?pxS2?(r,x)x!=x=0∑r?Anx?pxS2?(r,x)
Poisson分布的高階矩
Poisson分布X~Pois(λ)X \sim Pois(\lambda)X~Pois(λ)的mass function為
fX(x)=λxe?λx!,x=0,1,?,∞f_X(x) = \frac{\lambda^xe^{-\lambda}}{x!},x=0,1,\cdots,\inftyfX?(x)=x!λxe?λ?,x=0,1,?,∞
下面計算它的rrr階原點矩
EXr=∑x=0∞xrλxe?λx!=∑x=0∞λxe?λx!Mx0r=e?λeλM0r=eλΔ0r=∑x=0rλxx!Δx0r=∑x=0rλxx!S2(r,x)x!=∑x=0rλxS2(r,x)EX^r = \sum_{x=0}^{\infty} x^r\frac{\lambda^xe^{-\lambda}}{x!} = \sum_{x=0}^{\infty} \frac{\lambda^xe^{-\lambda}}{x!}M^x0^r = e^{-\lambda}e^{\lambda M} 0^r \\= e^{\lambda\Delta} 0^r = \sum_{x=0}^r \frac{\lambda^x}{x!} \Delta^x 0^r = \sum_{x=0}^r \frac{\lambda^x}{x!} S_2(r,x)x! = \sum_{x=0}^r \lambda^xS_2(r,x)EXr=x=0∑∞?xrx!λxe?λ?=x=0∑∞?x!λxe?λ?Mx0r=e?λeλM0r=eλΔ0r=x=0∑r?x!λx?Δx0r=x=0∑r?x!λx?S2?(r,x)x!=x=0∑r?λxS2?(r,x)
負二項分布的高階矩
負二項分布X~NegBinom(r,p)X \sim NegBinom(r,p)X~NegBinom(r,p)的mass function為
fX(x)=Cr+x?1r?1pr(1?p)x,x=0,1,?,∞f_X(x) = C_{r+x-1}^{r-1}p^{r}(1-p)^x,x=0,1,\cdots,\inftyfX?(x)=Cr+x?1r?1?pr(1?p)x,x=0,1,?,∞
下面計算它的rrr階原點矩
EXr=∑x=0∞xrCr+x?1r?1pr(1?p)x=∑x=0∞Cr+x?1r?1pr(1?p)xMx0r=pr[1?(1?p)M]?r0r=pr[p?(1?p)Δ]?r0r=[1?1?ppΔ]?r0r=∑x=0rCr+x?1r?1(1?pp)xΔx0r=∑x=0rAr+x?1x(1?pp)xS2(r,x)EX^r = \sum_{x=0}^{\infty} x^rC_{r+x-1}^{r-1}p^{r}(1-p)^x = \sum_{x=0}^{\infty} C_{r+x-1}^{r-1}p^{r}(1-p)^xM^x0^r \\ = p^r[1-(1-p)M]^{-r}0^r = p^r[p-(1-p)\Delta]^{-r}0^r = [1-\frac{1-p}{p}\Delta]^{-r}0^r \\ = \sum_{x = 0}^r C_{r+x-1}^{r-1}(\frac{1-p}{p})^x \Delta ^x 0^r = \sum_{x = 0}^r A_{r+x-1}^{x}(\frac{1-p}{p})^x S_2(r,x) EXr=x=0∑∞?xrCr+x?1r?1?pr(1?p)x=x=0∑∞?Cr+x?1r?1?pr(1?p)xMx0r=pr[1?(1?p)M]?r0r=pr[p?(1?p)Δ]?r0r=[1?p1?p?Δ]?r0r=x=0∑r?Cr+x?1r?1?(p1?p?)xΔx0r=x=0∑r?Ar+x?1x?(p1?p?)xS2?(r,x)
總結
以上是生活随笔為你收集整理的UA MATH564 概率论 高阶矩的计算:有限差算子方法2的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 脑与认知科学1 脑的基本结构下
- 下一篇: UA MATH564 概率论 多项分布