如何有效地将 dpoibin 分解为 R 中的被加数？

Question

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

dpoibin命令，在poibin包中，可以得到质量概率函数。例如，使用此命令：

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

一个人可以获得在 100 次独立伯努利试验序列中获得 30 次成功的概率，成功概率包含在向量 Probs_Success 中。要计算此概率，必须对长度为 100 的所有可能序列的概率求和，其中有 30 次成功和 70 次失败。

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

对于引发此问题的问题感兴趣的人，请单击以下 link：

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

Answer 1

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

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

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

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

感谢所有给我建议的人。