如何有效地将Dpoibin分解为R中的求和项?

时间:2018-09-24 14:52:38

标签: r probability probability-theory probability-distribution bernoulli-probability

泊松二项式分布涉及一系列具有不同成功概率的独立伯努利试验中成功次数的概率。这是二项分布的一般化。

使用dpoibin包中的命令poibin,可以获得质量概率函数。例如,使用以下命令:

library(poibin)
n <- 100
Probs_Success <- runif(n)
dpoibin(kk = 30, pp = Probs_Success)
可以在向量Probs_Success中包含100个独立的伯努利试验序列中获得30个成功的概率。要计算此概率,必须将长度为100的所有可能序列的概率相加,其中有30个成功案例和70个失败案例。

问题:如何有效地获得在R中产生上述概率的所有求和项?非常感谢您的帮助。

对于那些引起该问题的问题感兴趣的人,请单击以下链接:

https://math.stackexchange.com/questions/2924831/bivariate-poisson-binomial-distribution

1 个答案:

答案 0 :(得分:0)

我解决了引发上述问题的问题。该解决方案基于论文

https://www.tandfonline.com/doi/abs/10.1080/03610918708812585

要查看此解决方案,只需单击以下链接

https://math.stackexchange.com/questions/2924831/bivariate-poisson-binomial-distribution

感谢所有向我提供建议的人。