Skip to contents

Move branching process model forward from generation g to generation g+1.

Usage

outbreak_step(
  case_data,
  g,
  hrc_mu,
  hrc_disp,
  lrc_mu,
  lrc_disp,
  prop_pep,
  rel_risk,
  p_hrc_case,
  p_lrc_case
)

Arguments

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

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

p_hrc_case

Probability that a high-risk contact becomes a case

p_lrc_case

Probability that a low-risk contact becomes a case

outbreak_data

A data.frame of cases up to generation g+1.

Value

A data.frame updated with one new generation of cases.

Author

Sophie Meakin, Joel Hellewell