Skip to content

Latest commit

 

History

History
82 lines (70 loc) · 3.39 KB

README.md

File metadata and controls

82 lines (70 loc) · 3.39 KB
Raw

Global Sensitivity Analysis (GSA) of Covid-19 Transmission on a University Campus

All code is written in MATLAB R2021b

Directories

  1. contact_matrices/

Folder contains the contact matrices for both campus contacts and dorm living contacts under different classroom and dorm capacity scenarios. The following key indicates the specifics of each contact scenario based on a case number 1 through 24.
Covid-19 model and transmission dynamics
Covid-19 transmission diagram.

  1. model_sols/

Model solutions created by the script fun_compute_model_solutions.m are saved in this folder.
Model solutions.

  1. sobol_indices_dt/

Sobol indices computed using the doubling time metric, with the script compute_sobol_indices_doubling_time.m, are saved to this folder.
GSA of DT.

  1. sobol_indices_in_time/

Sobol indices computed in time, using the cumulative number of cases computed with the script compute_sobol_indices_over_time.m are saved to this folder.

Scripts

  1. compute_sobol_indices_doubling_time.m

Computes the first and total-effect Sobol indices using the disease doubling time metric (see fun_doubling_time.m) and model solutions found in the model_sols/ folder. The sobol indices are saved in the sobol_indices_dt/ folder.

  1. compute_sobol_indices_in_time.m

Computes the first and total-effect Sobol indices using the cumulative number of cases (see fun_sobol_indices_by_solution_in_time.m) and the model solutions found in the model_sols/ folder. The resulting sobol indices are saved in the sobol_indices_in_time/ folder.

  1. RUN_ME.m

Shows examples of how to run the code to produce model solutions, perform the Sobol analysis and produce plots.

Functions

  1. fun_compute_model_solutions.m

Computed model solutions for particular contact scenarios (see contact_matrices) and model parameter ranges (see fun_model_parameter_ranges.m). All solutions are saved in the model_sols/ folder.

  1. fun_covid_model_ode.m

SEIR model of disease transmission dynamics on a university campus. This function is used by the compute_all_model_solutions.m script.

  1. fun_doubling_time.m

Computes the disease doubling time over a desired number of consecutive days since the start of a semester.

  1. fun_initialize_contact_matrices.m

Combines contact matrices found in the contact_matrices/ folder into arrays that can be accessed by case number (see Directories section).

  1. fun_lhs_sampling_of_parameters.m

Performs Latin hypercube sampling (lhs) of desired parameters in the ranges specified in the fun_model_parameter_ranges.m function. Model solutions are then computed using the sets of sampled parameter values.

  1. fun_model_parameter_ranges.m

Stores the user indicated model parameters ranges and has a binary option for whether a parameter is to be varied (true) or not (false).

  1. fun_sobol_indices_by_measure.m

Computes the first and total-effect Sobol indices using a measure function that is applied to model solutions stored in the model_solutions/ folder.

  1. fun_sobol_indices_by_solution_in_time.m

Computes the first and total-effect Sobol indices using the model solutions saved in the model_solutions/ folder.