多项式分布的正态近似是多少?


Answers:


21

您可以使用多元正态分布对其进行近似,就像通过二元正态分布近似二项分布一样。检查分布理论和多项式分布的元素,第15-16-17页。

设是你的概率的向量。然后多元正态分布的均值向量是Ñ p = (Ñ p 1,Ñ p 2,。。。,Ñ p ķ)。协方差矩阵是k × k对称矩阵。对角元素实际上是X i的方差;即ñ p 我P=(p1,...,pk)np=(np1,np2,...,npk)k×kXi,我= 1 ,2 ... ,ķ。第i行第j列的对角线元素为 Cov (X i,X j)= - n p i p j,其中 i不等于 j。npi(1−pi)i=1,2...,kCov(Xi,Xj)=−npipjij


1
查看第二参考。
— 统计

3
统计信息,以便该答案可以独立存在(并且可以防止链接腐烂),您介意给出解决方案的摘要吗?
— ub

4
这需要连续性校正吗?您将如何应用它?
— 杰克·艾德利2014年

2
协方差矩阵不是正定的,而是正半定的,并且不是满秩的。这使得所得的多正态分布不确定。这是我面临的问题。知道如何处理吗?
— Mohammad Alaggan '16

2
@ M.Alaggan:这里定义的均值/协方差矩阵有一个小问题:对于具有变量的多项式分布,等效多变量正态具有k − 1个变量。这在简单的二项式示例中很明显,该示例由(普通)正态分布近似。有关更多讨论,请参见《分布理论的要素》示例12.7 。kk−1
— MS Dousti

1

此答案中给出的密度是退化的,因此我使用以下公式来计算由正态近似得出的密度:

有写着给定一个随机变量定理X=[X1,…,Xm]T∼Multinom(n,p),对于m维矢量p与∑ipi=1和∑iXi=n,那个;

X→dndiag(u)Q[Z1⋮Zm−10]+[np1⋮npm],

n

  • uui=pi
  • Zi∼N(0,1)i=1,…,m−1
  • Qu

m−1m−1XXm

QI−2vvTvi=(δim−ui)/2(1−um)

m−1Qm−1m−1X^Q^

X^→dndiag(u^)Q^[Z1⋮Zm−1]+[np1⋮npm−1]∼N(μ,nΣ),

n

  • u^m−1u
  • μ=[np1,…,npm−1]T
  • nΣ=nAATA=diag(u^)Q^

该最终方程式的右侧是计算中使用的非简并密度。

不出所料,当您插入所有内容时,您将获得以下协方差矩阵:

(nΣ)ij=npipj(δij−pipj)

i,j=1,…,m−1m−1m−1

此博客条目是我的起点。


1
另一个有用的资源是中提供的链接:stats.stackexchange.com/questions/2397/...
— stephematician

1
好的答案(+1)---请注意,您可以使用语法嵌入链接[textual description](hyperlink)。我已自由编辑此答案以嵌入您的链接。
— 本-恢复莫妮卡
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.