Get parameters for vaccine efficacy and antibody titre decay rate.

get_vaccine_ab_titre_parameters(
  vaccine,
  max_dose = 2,
  correlated = FALSE,
  hl_s = 108,
  hl_l = 3650,
  period_s = 250,
  t_period_l = 365,
  ab_50 = 0.2,
  ab_50_severe = 0.03,
  std10 = 0.44,
  k = 2.94,
  nt_transmission_factor = 12,
  nt_efficacy_transmission = FALSE,
  max_ab = 8,
  mu_ab_list = data.frame(name = c("Pfizer", "AstraZeneca", "Sinovac", "Moderna"),
    mu_ab_d1 = c(13/94, 1/59, 28/164, ((185 + 273)/2)/321), mu_ab_d2 = c(223/94, 32/59,
    28/164, 654/158))
)

Arguments

vaccine

which vaccine? should be one of: "Pfizer", "AstraZeneca", "Sinovac", "Moderna" (controls the mean titre associated with each dose)

max_dose

maximum number of doses

correlated

are doses correlated?

hl_s

Half life of antibody decay - short

hl_l

Half life of antibody decay - long

period_s

length of the initial decay rate (shortest half-life)

t_period_l

Time point at which to have switched to longest half-life

ab_50

titre relative to convalescent required to provide 50% protection from infection, on linear scale

ab_50_severe

titre relative to convalescent required to provide 50% protection from severe disease, on linear scale

std10

Pooled standard deviation of antibody level on log10 data

k

shape parameter of efficacy curve

nt_transmission_factor

used by vaccine_efficacy_transmission to compute the effect of antibody titre on onward transmission

nt_efficacy_transmission

a logical variable, specify if antibody titre should affect onward transmission

max_ab

maximum allowable antibody titre draw (on natural log scale)

mu_ab_list

a data.frame