bogaotory / comoodecpp Goto Github PK

Hi CoMo team. I can't see contributing guidelines on the repo so not sure the best way to make a pull request - so I'll raise an issue instead.

In comoOde.cpp the derivative for the QVR compartment does not include any quarantine release. Does it need an additional term releasing those from quarantine?

For instance, should the derivative be the following?

      arma::vec dQVRdt = quarantine_rate*VR
                        - (1.0-vaccine_eff_r)*lam%QVR
                        - vac_dur_r*QVR
                        - omega*QVR
                        + ageing*QVR
                        - mort_col%QVR 
                        - (1.0/quarantine_days)*QVR // <------------------- this term
                        ;

Thanks.

I wanted to enquire about the QVR and QEVR compartments. I understand the QVR compartment represents the group of individuals that are under quarantine (among other things, i.e. vacc and prev infection). Given they're under quarantine, should the lambda term they're subjected to be lamq instead of lam?

The terms in dQVRdt (1 below) and dQEVRdt (2 below) would then balance (from L1029 of comoOde.cpp, here) and would be in line with other transitions (i.e. QV -> QEV, QR -> QER, QS -> QE).

      arma::vec dQVRdt = quarantine_rate*VR
                        - (1.0-vaccine_eff_r)*lam%QVR  // <--------------------- 1
                        - vac_dur_r*QVR
                        - omega*QVR
                        + ageing*QVR
                        - mort_col%QVR 
                        ;

      arma::vec dQEVRdt = quarantine_rate*EVR
                        - gamma*QEVR
                        + ageing*QEVR
                        - mort_col%QEVR
                        - (1.0/quarantine_days)*QEVR
                        + (1.0-vaccine_eff_r)*lamq%QVR // <--------------------- 2
                        ;

Thanks for the advice.

This feature is planned for "School Closure" and "Vaccination" only.
Maximum 2 age groups may have different coverage rates per day.

Progress:

• School Closure

  • R version
  • Cpp version

• Vaccination

• Online