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


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 iX j= - n p i p j,其中 i不等于 jnpi(1pi)i=1,2...,kCov(Xi,Xj)=npipjij


1
查看第二参考。
统计

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

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

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

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

1

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

有写着给定一个随机变量定理X=[X1,,Xm]TMultinom(n,p),对于m维矢量pipi=1iXi=n,那个;

Xdndiag(u)Q[Z1Zm10]+[np1npm],

n

  • uui=pi
  • ZiN(0,1)i=1,,m1
  • Qu

m1m1XXm

QI2vvTvi=(δimui)/2(1um)

m1Qm1m1X^Q^

X^dndiag(u^)Q^[Z1Zm1]+[np1npm1]N(μ,nΣ),

n

  • u^m1u
  • μ=[np1,,npm1]T
  • nΣ=nAATA=diag(u^)Q^

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

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

(nΣ)ij=npipj(δijpipj)

i,j=1,,m1m1m1

博客条目是我的起点。


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.