set.seed(27)
X <- PoissonBinomial(0.5, 0.3, 0.8)
X
mean(X)
variance(X)
skewness(X)
kurtosis(X)
random(X, 10)
pdf(X, 2)
log_pdf(X, 2)
cdf(X, 2)
quantile(X, 0.8)
cdf(X, quantile(X, 0.8))
quantile(X, cdf(X, 2))
## equivalent definitions of four Poisson binomial distributions
## each summing up three Bernoulli probabilities
p <- cbind(
p1 = c(0.1, 0.2, 0.1, 0.2),
p2 = c(0.5, 0.5, 0.5, 0.5),
p3 = c(0.8, 0.7, 0.9, 0.8))
PoissonBinomial(p)
PoissonBinomial(p[, 1], p[, 2], p[, 3])
PoissonBinomial(p[, 1:2], p[, 3])
Run the code above in your browser using DataLab