A comparison of gamma and lognormal maximum likelihood estimators in a sequential population analysis
We analyze the model used to assess most major commercial marine fish populations, namely, sequential population analysis (SPA). This model estimates population abundance by combining catch-at-age data with research surveys or commercial catch per unit effort indices of abundance. We examine two maximum likelihood estimators of SPA parameters. These estimators are based on assuming that the stock-size indices are from lognormal or gamma distributions. Using simulations, we find that both types of estimators can have significant biases; however, our results indicate that it is preferable to use the gamma model, because it tends to have lower bias and variability, even when the true distribution of the stock-size indices is lognormal.