forgot to log2 the nabs
@@ -1382,7 +1382,7 @@ class BaseVaccination(Intervention):
if self.p['doses'] == len(self.p['target_eff']):
# determine efficacy of first dose (assume efficacy supplied is against symptomatic disease)
nabs = np.arange(-8, 4, 0.1)
- VE_symp = cvi.calc_VE_symp(nabs, sim.pars['nab_eff'])
+ VE_symp = cvi.calc_VE_symp(2**nabs, sim.pars['nab_eff'])
peak_nab = nabs[np.argmax(VE_symp>self.p['target_eff'][0])]
self.p['nab_init'] = dict(dist='normal', par1=peak_nab, par2=2)
if self.p['doses'] == 2:
