## Maximum likelihood estimation for randomly generated data
## from the Poisson-Lindley distribution.
require(stats4)
set.seed(100)
pl.data <- rpoislind(n = 500, theta = 0.5)
out.pl <- poislind.ll(pl.data)
coef(out.pl)
vcov(out.pl)Run the code above in your browser using DataLab