Rate Matrix (Q) for Adult Mosquito SEI Dynamics — make_Q

Construct the infinitesimal generator matrix for (individual) adult female infection dynamics. Adult females follow SEI (Susceptible-Exposed-Infectious) style dynamics with a Gamma distributed EIP, with a mean duration 1/q and variance 1/nq^2 (following shape-scale parameterization, EIP ~ Gamma(n,1/nq)). This function only constructs the rate matrix for either a single mosquito or cohort that all emerged at the same time (the rate matrix for a population with emergence is infinite in dimension).

Usage

make_Q_Imperial(q, n, mu, FOIv)

Arguments

q: related to scale parameter of Gamma distributed EIP (1/q is mean length of EIP)
n: shape parameter of Gamma distributed EIP
mu: mosquito mortality rate
FOIv: equilibrium force of infection on mosquitos

Value

rate matrix for a single (emergence) cohort of SEI mosquito