Skip to contents

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)