r - 如何在 R 中有效地将 dpoibin 分解为其被加数？

Question

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

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

library(poibin)
n <- 100
Probs_Success <- runif(n)
dpoibin(kk = 30, pp = Probs_Success)

人们可以获得在 100 次独立伯努利试验序列中获得 30 次成功的概率，成功概率包含在向量 Probs_Success 中。为了计算这个概率，我们必须将长度为 100 的所有可能序列的概率相加，其中有 30 次成功，70 次失败。

问题:如何在 R 中有效地获取生成上述概率的所有被加数？非常感谢您的帮助。

对于那些对引发此问题的问题感兴趣的人，请点击以下链接:

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

Answer 1

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

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

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

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

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