Parameter distributions for N and p were calculated using MCMC simulation using a standard Metropolis algorithm with flat prior information. It was necessary to compute a correlation matrix to define a joint proposal since N and p are closely correlated. We monitored convergence using two chains. Once they had converged, we had a burn in period of 10 5 samples.