#fit a finite Poisson mixture model of order 2 to this data set,
#plot the histogram of the observations and the fitted probability mass function.
data(earthquake)
out <- pmle.pois(earthquake,2,1)
plotmix.pois(earthquake,out)
Run the code above in your browser using DataLab