EfficacyCI {vaxpmx} | R Documentation |
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)
efficacySet |
numeric vector - vector of estimated efficacy values |
ci |
numeric - required confidence level |
Confidence intervals are calculated using quantiles of estimated efficacy values.
named list - mean, median, CILow, CIHigh
# 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, specific to serotype 2
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)