pcount function

Fit the N-mixture model of Royle (2004)