Skip to content

Commit

Permalink
Fix seeding time
Browse files Browse the repository at this point in the history
  • Loading branch information
romain-ragonnet committed Aug 31, 2023
1 parent 1ffa59e commit 28a11e6
Showing 1 changed file with 11 additions and 6 deletions.
17 changes: 11 additions & 6 deletions autumn/projects/sm_covid2/common_school/project_maker.py
Original file line number Diff line number Diff line change
Expand Up @@ -214,10 +214,15 @@ def get_school_project_parameter_set(iso3, first_date_with_death, sero_age_min,
}
baseline_params = baseline_params.update(sero_age_params)

# Set seeding time 40 days prior first reported death
baseline_params = baseline_params.update(
{"infectious_seed_time": first_date_with_death - 40.}
)

# update using MLE params, if available
mle_path= param_path / "mle_files" / f"mle_{iso3}.yml"
if exists(mle_path):
baseline_params = baseline_params.update(mle_path, calibration_format=True)
# mle_path= param_path / "mle_files" / f"mle_{iso3}.yml"
# if exists(mle_path):
# baseline_params = baseline_params.update(mle_path, calibration_format=True)

# update using potential Sensitivity Analysis params
sa_params_path = param_path / "SA_analyses" / f"{analysis}.yml"
Expand Down Expand Up @@ -360,12 +365,12 @@ def get_school_project_priors(first_date_with_death):
"""

# Work out max infectious seeding time so transmission starts before first observed deaths
min_seed_time = first_date_with_death - 100
max_seed_time = first_date_with_death - 1
# min_seed_time = first_date_with_death - 100
# max_seed_time = first_date_with_death - 1

priors = [
UniformPrior("contact_rate", [0.01, 0.06]),
UniformPrior("infectious_seed_time", [min_seed_time, max_seed_time]),
# UniformPrior("infectious_seed_time", [min_seed_time, max_seed_time]),
UniformPrior("age_stratification.ifr.multiplier", [0.5, 1.5]),

# VOC-related parameters
Expand Down

0 comments on commit 28a11e6

Please sign in to comment.