N <- 100 # sample size
x <- 10 # sample deaths within 30 days
pi <- 0.08 # national proportion of death
lambda <- (N * pi)
1 - ppois( q = 10,
lambda = 8,
lower.tail = TRUE
)
## [1] 0.1841142
?dbinom
dbinom(x = 10,
size = 100,
prob = 0.08
)
## [1] 0.102355
There is higher probability for deaths to occur within the 30 days in the hospital than the national average. There is discrepancy between the answers given by the poisson and binomial methods.