Towards a replicator dynamics model of age structured populations

We present a new modelling framework combining replicator dynamics, the standard model of frequency dependent selection, with an age-structured population model. The new framework allows for the modelling of populations consisting of competing strategies carried by individuals who change across their life cycle. Firstly the discretization of the McKendrick von Foerster model is derived. We show that the Euler–Lotka equation is satisfied when the new model reaches a steady state (i.e. stable frequencies between the age classes). This discretization consists of unit age classes where the timescale is chosen so that only a fraction of individuals play a single game round. This implies a linear dynamics and individuals not killed during the round are moved to the next age class; linearity means that the system is equivalent to a large Bernadelli–Lewis–Leslie matrix. Then we use the methodology of multipopulation games to derive two, mutually equivalent systems of equations. The first contains equations describing the evolution of the strategy frequencies in the whole population, completed by subsystems of equations describing the evolution of the age structure for each strategy. The second contains equations describing the changes of the general population’s age structure, completed with subsystems of equations describing the selection of the strategies within each age class. We then present the obtained system of replicator dynamics in the form of the mixed ODE-PDE system which is independent of the chosen timescale, and much simpler. The obtained results are illustrated by the example of the sex ratio model which shows that when different mortalities of the sexes are assumed, the sex ratio of 0.5 is obtained but that Fisher’s mechanism, driven by the reproductive value of the different sexes, is not in equilibrium.