EfficacyCI function

Efficacy summary (mean, median, confidence intervals)

Function summarizes efficacy statistics (mean, median, confidence intervals) based on the set of estimated efficacy values and chosen condfidence interval.

EfficacyCI(efficacySet, ci = 0.95)

Arguments

• efficacySet: numeric vector - vector of estimated efficacy values
• ci: numeric - required confidence level

Returns

named list - mean, median, CILow, CIHigh

Details

Confidence intervals are calculated using quantiles of estimated efficacy values.

Examples

# Load required packages
library(dplyr)

# Load an example dataset
data(data_temp)
Data.vaccinated <- filter(data_temp, vaccine == 1)
Data.control <- filter(data_temp, vaccine == 0)

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

# Estimate 95\% confidence interval of vaccine efficacy based on the fitted model
efficacySet <- glmParametricSampling(logisticFit, nboot = 500, Data.vaccinated, Data.control)
EfficacyCI(efficacySet)