sus_ORs no longer averaged for 10-20 and 60-70
@@ -194,7 +194,7 @@ def get_prognoses(by_age=True):
else:
prognoses = dict(
age_cutoffs = np.array([0, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80]), # Age cutoffs (lower limits)
- sus_ORs = np.array([0.34, 0.34, 0.34, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1.47, 1.47, 1.47, 1.47]), # Odds ratios for relative susceptibility -- from https://science.sciencemag.org/content/early/2020/05/04/science.abb8001; 10-20 and 60-70 bins are the average across the ORs
+ sus_ORs = np.array([0.34, 0.34, 0.34, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1.47, 1.47, 1.47, 1.47]), # Odds ratios for relative susceptibility -- from https://science.sciencemag.org/content/early/2020/05/04/science.abb8001;
trans_ORs = np.array([1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]), # Odds ratios for relative transmissibility -- no evidence of differences
comorbidities = np.array([1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]), # Comorbidities by age -- set to 1 by default since already included in disease progression rates
symp_probs = np.array([0.5, 0.5, 0.55, 0.55, 0.6, 0.6, 0.65, 0.65, 0.7, 0.7, 0.75, 0.75, 0.8, 0.8, 0.85, 0.85, 0.9]), # Overall probability of developing symptoms (based on https://www.medrxiv.org/content/10.1101/2020.03.24.20043018v1.full.pdf, scaled for overall symptomaticity)
