A Function to Generate Scenario Results from the Sutherland Model

A Function to Generate Scenario Results from the Sutherland Model

gen_scenario_results(scenarios, parameter_samples, cores = 1, ...)

Arguments

scenarios	A dataframe of scenarios, see the examples for a structural example.
parameter_samples	A dataframe of parameter samples, see the examples for an example.
cores	The number of cores to use, defaults to one
...	Additional arguments passed to sutherland_model.

Value

A dataframe combining the scenarios, parameter samples, and model results

Examples

scenarios <- gen_annual_change_scenarios(TB_decrease = c(0.09, 0.05, 0.01),
cohorts_to_sim = 11, samples = 10)

parameter_samples <- sample_model_parameters(incidence_rates = sutherland_incidence_rates,
population_data = sutherland_data, samples = 10)

gen_scenario_results(scenarios, parameter_samples, cores = 1)
#> # A tibble: 30 x 20
#>    scenario sample annual_TB_decre… incidence_rates per_year_one sym_lag
#>    <chr>     <int> <list>           <list>                 <dbl>   <dbl>
#>  1 9% decr…      1 <dbl [55 × 3]>   <list [5]>             0.764    1.32
#>  2 9% decr…      2 <dbl [55 × 3]>   <list [5]>             0.764    1.46
#>  3 9% decr…      3 <dbl [55 × 3]>   <list [5]>             0.764    1.23
#>  4 9% decr…      4 <dbl [55 × 3]>   <list [5]>             0.764    1.44
#>  5 9% decr…      5 <dbl [55 × 3]>   <list [5]>             0.764    1.49
#>  6 9% decr…      6 <dbl [55 × 3]>   <list [5]>             0.764    1.54
#>  7 9% decr…      7 <dbl [55 × 3]>   <list [5]>             0.764    1.28
#>  8 9% decr…      8 <dbl [55 × 3]>   <list [5]>             0.764    1.42
#>  9 9% decr…      9 <dbl [55 × 3]>   <list [5]>             0.764    1.42
#> 10 9% decr…     10 <dbl [55 × 3]>   <list [5]>             0.764    1.42
#> # … with 20 more rows, and 14 more variables: `Table 2 - Estimated risk of
#> #   developing notified TB` <list>, `Table 3 - Estimated of no. of TB
#> #   notifications prevented` <list>, `Table 4 - Estimated no. of TB
#> #   notifications prevented by Schools BCG scheme` <list>, `Table 4 - Total
#> #   notifcations prevented each year` <list>, `Table 5 - Primary additional
#> #   notifications` <list>, `Table 5 - Secondary additional notifcations reduced
#> #   to Sutherland for clarity` <list>, `Table 5 - Secondary additional
#> #   notifications` <list>, `Table 5 - Total Effects of ending the schools BCG
#> #   scheme` <list>, `Table 5 - Total notifications if the scheme
#> #   continues` <list>, `Table 5 - Total secondary notifications all
#> #   time` <list>, `Total additional notifications` <list>, `Total additional
#> #   notifications all time` <list>, total_primary_additional <list>,
#> #   total_secondary_additional <list>