make_Q_Imperial {MGDrivE2}

Rate Matrix (Q) for Adult Mosquito SEI Dynamics

Description

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