Распределение Пуассона-Бинома касается вероятности числа успехов в последовательности независимых испытаний Бернулли с различными вероятностями успеха.Это обобщение биномиального распределения.
С помощью команды dpoibin в пакете poibin можно получить функцию массовой вероятности.Например, с помощью этой команды:
library(poibin)
n <- 100
Probs_Success <- runif(n)
dpoibin(kk = 30, pp = Probs_Success)
можно получить вероятность получения 30 успехов в последовательности из 100 независимых испытаний Бернулли с вероятностями успеха, содержащимися в векторе Probs_Success.Чтобы вычислить эту вероятность, нужно сложить вероятности всех возможных последовательностей длины 100, где есть 30 успехов и 70 неудач.
Вопрос: Как я могу эффективно получить все слагаемые, которые генерируют вышеупомянутую вероятность в R?Большое спасибо за вашу помощь.
Для тех, кто заинтересован в проблеме, которая мотивировала этот вопрос, пожалуйста, нажмите на следующую ссылку:
https://math.stackexchange.com/questions/2924831/bivariate-poisson-binomial-distribution