Last chance! 50% off unlimited learning
Sale ends in
dpospois(x, lambda)
ppospois(q, lambda)
qpospois(p, lambda)
rpospois(n, lambda)
dpospois
gives the density,
ppospois
gives the distribution function,
qpospois
gives the quantile function, and
rpospois
generates random deviates.lambda
is not permitted here.pospoisson
,
zapoisson
,
rpois
.lambda = 2
y = rpospois(n=1000, lambda)
table(y)
mean(y) # Sample mean
lambda / (1-exp(-lambda)) # Population mean
(i = dpospois(0:7, lambda))
cumsum(i) - ppospois(0:7, lambda) # Should be 0s
table(rpospois(100, lambda))
table(qpospois(runif(1000), lambda))
round(dpospois(1:10, lambda) * 1000) # Should be similar
x = 0:7
plot(x, dpospois(x, lambda), type="h", ylab="Probability",
main=paste("Positive Poisson(", lambda, ") (blue) vs",
" Poisson(", lambda, ") (red and shifted slightly)", sep=""),
lwd=2, col="blue", las=1)
lines(x+0.05, dpois(x, lambda), type="h", lwd=2, col="red")
Run the code above in your browser using DataLab