A two-sex matrix population model to represent harem structure

Population dynamic models often include males in the calculation of population change, but even in those cases males have rarely been introduced to represent polygyny (harem social structure), where it is particularly important to include males in the reproductive performance of the population. In this article we develop an adaptable matrix population modeling framework for species that have a harem-like social structure under an assumption that the transitions from newborn to juvenile and juvenile to adult both take one time step. We are able to calculate not only the growth rates and stable stage distributions, but also the mathematical expressions for harem size for this model. We then provide applications of this model to two mammal species with slightly different harem behavior.