Finding Reproduction Numbers for Epidemic Models & Predator-Prey Models of Arbitrary Finite Dimension Using The Generalized Linear Chain Trick

TL;DR

This paper introduces a method to derive general reproduction number formulas for diverse epidemic and ecological models of varying dimensions using the generalized linear chain trick and the next generation operator approach.

Contribution

It presents a novel approach combining the generalized linear chain trick with the next generation operator to compute reproduction numbers across model families of different sizes.

Findings

Derived general reproduction number formulas for finite-dimensional models.

Connected mean-field ODE models with CTMCs and phase-type distributions.

Applied method to generalized SEIRS and predator-prey models.

Abstract

Reproduction numbers, like the basic reproduction number $\mathcal{R}_0$, play an important role in the analysis and application of dynamic models, including contagion models and ecological population models. One difficulty in deriving these quantities is that they must be computed on a model-by-model basis, since it is typically impractical to obtain general reproduction number expressions applicable to a family of related models, especially if these are of different dimensions. For example, this is typically the case for SIR-type infectious disease models derived using the linear chain trick (LCT). Here we show how to find general reproduction number expressions for such models families (which vary in their number of state variables) using the next generation operator approach in conjunction with the generalized linear chain trick (GLCT). We further show how the GLCT enables modelers to draw insights from these results by leveraging theory and intuition from continuous time Markov chains (CTMCs) and their absorption time distributions (i.e., phase-type probability distributions). To do this, we first review the GLCT and other connections between mean-field ODE model assumptions, CTMCs, and phase-type distributions. We then apply this technique to find reproduction numbers for two sets of models: a family of generalized SEIRS models of arbitrary finite dimension, and a generalized family of finite dimensional predator-prey (Rosenzweig-MacArthur type) models. These results highlight the utility of the GLCT for the derivation and analysis of mean field ODE models, especially when used in conjunction with theory from CTMCs and their associated phase-type distributions.