Skip to contents

Run a specified number of simulations with identical parameters.

Usage

scenario_sim(
  n_sim,
  prop_pep,
  rel_risk,
  hrc_mu,
  hrc_disp,
  lrc_mu,
  lrc_disp,
  p_hrc_case,
  p_lrc_case,
  n_initialcases,
  cap_max_gen,
  cap_max_cases
)

Arguments

prop_pep

Proportion of high-risk contacts that receive PEP

rel_risk

Relative risk of a high-risk contact becoming a case after receiving PEP

hrc_mu

Average number of high-risk contacts

hrc_disp

Dispersion parameter for number of high-risk contacts

lrc_mu

Average number of low-risk contacts

lrc_disp

Dispersion parameter for number of low-risk contacts

p_hrc_case

Probability that a high-risk contact becomes a case

p_lrc_case

Probability that a low-risk contact becomes a case

n_initialcases

Integer number of initial cases

cap_max_gen

Stop the simulation when this many generations is reached

cap_max_cases

Stop the simulation when this many cases is reached

nsim

Number of simulations to run

Value

A data.frame containing the results for multiple simulations.

Author

Joel Hellewell, Sophie Meakin