A survival-analysis-based simulation model for Russian wheat aphid population dynamics

A simulation model for Russian wheat aphid (RWA), Diuraphis noxia (Mordvilko), populations is built by integrating survival-analysis-based development and survivor functions and the same-shape reproduction distribution model in the framework of Leslie [Leslie, P.H., 1945. On the use of matrices in certain population mathematics. Biometrika 33, 183–212] matrix structure. Survival analysis is utilized to model both the development and survival of RWA populations, and the Cox (1972) proportional hazards model is fitted with the data sets from our laboratory observation of 1800 RWA individuals under 25 factorial combinations of five temperature regimes and five barley plant-growth stages. Rather than using simple age-specific survivor rates as in the traditional Leslie matrix, the survivor functions based on survival analysis describe age-specific, temperature and plant stage-dependent RWA survival probabilities. Similarly, a probability model from survival analysis to estimate the probability that an individual will reach mature adult stage is utilized to describe the development process; this makes the transition from nymphal stage to mature adult stage dependent on RWA age as well as temperature and plant-growth stage.Inspired by the same-shape distribution and rate-summation approach for modeling insect development, a similar approach for modeling insect reproduction under variable temperature is developed. This new same-shape reproduction distribution model incorporates individual variation in reproduction capability, as well as the effects of RWA age, temperature and plant-growth stage. Consequently, the same-shape reproduction distribution model replaces the simple age-specific fecundities in Leslie matrix model. To the best of our knowledge, this work is the first to introduce survival analysis to simulation modeling in entomology and ecology and also the first to integrate our newly developed same-shape reproduction distribution model into application.