Run a specified number of simulations with identical parameters

scenario_sim(
  n,
  initial_cases,
  offspring,
  delays,
  event_probs,
  interventions,
  sim
)

Arguments

n: a positive integer scalar: number of simulations to run
initial_cases: a non-negative integer scalar: number of initial or starting cases which are all assumed to be missed.
offspring: a list with class <ringbp_offspring_opts>: the offspring distribution functions for the ringbp model, returned by offspring_opts(). Contains three elements: community, isolated, and asymptomatic
delays: a list with class <ringbp_delay_opts>: the delay distribution functions for the ringbp model, returned by delay_opts(). Contains two elements: incubation_period and onset_to_isolation
event_probs: a list with class <ringbp_event_prob_opts>: the event probabilities for the ringbp model, returned by event_prob_opts(). Contains three elements: asymptomatic, presymptomatic_transmission and symptomatic_ascertained
interventions: a list with class <ringbp_intervention_opts>: the intervention settings for the ringbp model, returned by intervention_opts(). Contains one element: quarantine
sim: a list with class <ringbp_sim_opts>: the simulation control options for the ringbp model, returned by sim_opts()

Value

A data.table object returning the results for multiple simulations using the same set of parameters. The table has columns

week: The week in the simulation.
weekly_cases: The number of new cases that week.
cumulative: The cumulative cases.
effective_r0: The effective reproduction rate for the whole simulation
cases_per_gen: A list column with the cases per generation. This is repeated each row.
sim: Index column for which simulation.

Author

Joel Hellewell

Examples

offspring <- offspring_opts(
  community = \(n) rnbinom(n = n, mu = 2.5, size = 0.16),
  isolated = \(n) rnbinom(n = n, mu = 0, size = 1),
  asymptomatic = \(n) rnbinom(n = n, mu = 2.5, size = 0.16)
)
delays <- delay_opts(
  incubation_period = \(n) rweibull(n = n, shape = 2.32, scale = 6.49),
  onset_to_isolation = \(n) rweibull(n = n, shape = 2.5, scale = 5)
)
event_probs <- event_prob_opts(
  asymptomatic = 0,
  presymptomatic_transmission = 0.3,
  symptomatic_ascertained = 0
)
interventions <- intervention_opts(quarantine = TRUE)
sim <- sim_opts(
  cap_max_days = 365,
  cap_cases = 2000
)
res <- scenario_sim(
  n = 5,
  initial_cases = 5,
  offspring = offspring,
  delays = delays,
  event_probs = event_probs,
  interventions = interventions,
  sim = sim
)
res
#>        sim  week weekly_cases cumulative effective_r0
#>      <int> <num>        <num>      <num>        <num>
#>   1:     1     0            5          5     0.000000
#>   2:     1     1            0          5     0.000000
#>   3:     1     2            0          5     0.000000
#>   4:     1     3            0          5     0.000000
#>   5:     1     4            0          5     0.000000
#>  ---                                                 
#> 261:     5    48            0       2171     2.900912
#> 262:     5    49            0       2171     2.900912
#> 263:     5    50            0       2171     2.900912
#> 264:     5    51            0       2171     2.900912
#> 265:     5    52            0       2171     2.900912
#>                    cases_per_gen
#>                           <list>
#>   1:                           0
#>   2:                           0
#>   3:                           0
#>   4:                           0
#>   5:                           0
#>  ---                            
#> 261:   4, 46, 35, 69,154,215,...
#> 262:   4, 46, 35, 69,154,215,...
#> 263:   4, 46, 35, 69,154,215,...
#> 264:   4, 46, 35, 69,154,215,...
#> 265:   4, 46, 35, 69,154,215,...