SwissTPH / openmalaria

A simulator of malaria epidemiology and control
GNU General Public License v2.0
75 stars 30 forks source link

2 or 3 compartment PK Model Concentrations #225

Closed sharma-raman closed 5 years ago

sharma-raman commented 5 years ago

PPQ 2 compartment PK model using Hodel 2013 parameterisation (from the original paper) doesn’t give the same concentration-time plot as standalone 2 compartment PK model. The stand-alone model has been verified using independently implemented 2 compartment model both in the analytic form (Monolix 3 compartment model with k13=k31=0) and ODE form.

Although the Monolix intermediate analytic equations found in LSTMDrugThreeComp::updateCached(double bm) give the correct values when given the correct peripheral rate constants, follow on steps to calculate concentrations seem to differ from Mentre et al 2008 in a significant manner.

sharma-raman commented 5 years ago

After alterations, concentrations now match those outputted by stand-alone models.