Simulations of a mortality plateau in the sexual Penna model for biological aging

The Penna model is a strategy to simulate the genetic dynamics of age-structured populations, in which the individual genomes are represented by bit strings. It provides a simple metaphor for the evolutionary process in terms of the mutation accumulation theory. In its original version, an individual dies due to inherited diseases when its current number of accumulated mutations, n, reaches a threshold value T. Since the mean number of diseases increases with age, the probability to die is zero for very young ages (n < T) and equals 1 for the old ones (n > or = T). Here, instead of using a step function to determine the genetic death age, we test several other functions that may or may not slightly increase the death probability at young ages (n < T), but that decrease this probability at old ones. Our purpose is to study the oldest old effect, that is, a plateau in the mortality curves at advanced ages. By imposing certain conditions, it has been possible to obtain a clear plateau using the Penna model. However, a more realistic one appears when a modified version, that keeps the population size fixed without fluctuations, is used. We also find a relation between the birth rate, the age structure of the population, and the death probability.