Implements the MY component using a BRQS (Blood-feeding / Resting / gravid Queue / Sugar-feeding) model of adult mosquito ecology and infection dynamics. Mosquitoes cycle through four behavioural states, and infection is tracked across all states with rate \(\tau\) governing progression from exposed (y) to infectious (z).
Su, Bu, Ru, Quuninfected mosquitoes in the sugar-feeding, blood-feeding, resting, and gravid states
Sy, By, Ry, Qyinfected (exposed) mosquitoes in each behavioural state
Sz, Bz, Rz, Qzinfectious mosquitoes in each behavioural state
fblood feeding rate
qhuman blood fraction
gbackground mortality rate
rhorate of leaving the sugar-feeding state to blood-feed
zetarate at which blood-feeding mosquitoes take a sugar meal
xirate of leaving the resting state
etarate of transition from resting to gravid
nu, thetaparameters governing the egg-laying / gravid state
deltafraction of gravid mosquitoes returning to blood feeding
omegafraction of emerging mosquitoes entering the sugar-feeding state
taurate of progression from exposed to infectious (\(1/\tau\) = EIP)
sigma_b, sigma_q, sigma_semigration rates in each state
$$ \begin{array}{rl} dB_u/dt &= (1-\omega)\Lambda + \rho S_u + (\nu+\theta)\delta Q_u + \xi R - (f+\zeta) B_u - \Omega_b \cdot B_u \\ dR_u/dt &= f(1-q\kappa) B_u - (\xi + \eta + g) R_u \\ dQ_u/dt &= \eta R_u - (\nu + \theta) Q_u - \Omega_q \cdot Q_u \\ dS_u/dt &= \omega\Lambda + \zeta B_u + (\nu+\theta)(1-\delta) Q_u - \rho S_u - \Omega_s \cdot S_u \\ \end{array} $$ with analogous equations for the infected (y) and infectious (z) classes, where the term \(\tau\) governs progression between infection stages.