The BQ module for the MY component

Source: R/adult-MY-BQ.R

Implements the MY component using a BQ (Blood-feeding / egg-laying Queue) model of adult mosquito ecology and infection dynamics. Mosquitoes alternate between a blood-feeding state (B) and a gravid / egg-laying state (Q), and infection is tracked across both states.

State Variables

Bu

density of uninfected, blood-feeding mosquitoes

Qu

density of uninfected, gravid (egg-laying) mosquitoes

By

density of infected (exposed), blood-feeding mosquitoes

Qy

density of infected (exposed), gravid mosquitoes

Bz

density of infectious, blood-feeding mosquitoes

Qz

density of infectious, gravid mosquitoes

Parameters

f

blood feeding rate

q

human blood fraction

g

mosquito mortality rate

sigma_b

emigration rate during blood feeding

sigma_q

emigration rate during egg laying

mu

fraction lost during emigration

nu

oviposition rate (per mosquito)

eip

extrinsic incubation period (\(\tau\)); \(\phi = 1/\tau\)

eggsPerBatch

eggs laid per oviposition bout

Dynamics

$$ \begin{array}{rl} dB_u/dt &= \Lambda + \nu Q_u - f B_u - \Omega_b \cdot B_u \\ dQ_u/dt &= f(1 - q\kappa) B_u - \nu Q_u - \Omega_q \cdot Q_u \\ dB_y/dt &= \nu Q_y - (f + \phi) B_y - \Omega_b \cdot B_y \\ dQ_y/dt &= fq\kappa B_u + f B_y - (\nu + \phi) Q_y - \Omega_q \cdot Q_y \\ dB_z/dt &= \phi B_y + \nu Q_z - f B_z - \Omega_b \cdot B_z \\ dQ_z/dt &= \phi Q_y + f B_z - \nu Q_z - \Omega_q \cdot Q_z \\ \end{array} $$