The Weibull model. — weibull_model • serosv

Refers to section 6.1.2.

Usage

weibull_model(t, status = NULL, pos = NULL, tot = NULL)

Arguments

t: the time vector.
status: the serostatus vector (optional if pos & tot are provided).
pos: the positive count vector (optional if status is provided).
tot: the total count vector (optional if status is provided).

Value

list of class weibull_model with the following items

datatype: type of datatype used for model fitting (aggregated or linelisting)
df: the dataframe used for fitting the model
info: fitted "glm" object
sp: seroprevalence
foi: force of infection

See also

[stats::glm()] for more information on the fitted "glm" object

Examples

df <- hcv_be_2006[order(hcv_be_2006$dur), ]
model <- weibull_model(
  t=df$dur,
  status=df$seropositive
  )
plot(model)