ve {vaxpmx}

Vaccine efficacy estimation

Description

Calculates vaccine efficacy and confidence interval as described in Dudasova et al., 2024, BMC Med Res Methodol

Usage

ve(Fit, Data, nboot = 2000)

Arguments

Fit

an object of class inheriting from "glm" representing the fitted model

Data

a data frame containing the variables in the fitted model; data must include a column called "vaccine" with binary indicator of vaccination status

nboot

a numeric value for number of bootstrap samples for confidence interval construction

Value

a value of vaccine efficacy VE and lower and upper bound of confidence interval CI

Examples

# Load an example dataset
data(data_temp)

# Fit logistic model relating neutralizing titer to disease status
logisticFit <- glm(disease_any ~ nAb1, data = data_temp, family = binomial())

# Estimate vaccine efficacy and 95\% confidence interval based on the fitted model
ve(logisticFit, data_temp, nboot = 500)