EfficacyCI {vaxpmx} | R Documentation |
Efficacy summary (mean, median, confidence intervals)
Description
Function summarizes efficacy statistics (mean, median, confidence intervals) based on the set of estimated efficacy values and chosen condfidence interval.
Usage
EfficacyCI(efficacySet, ci = 0.95)
Arguments
efficacySet |
numeric vector - vector of estimated efficacy values |
ci |
numeric - required confidence level |
Details
Confidence intervals are calculated using quantiles of estimated efficacy values.
Value
named list - mean, median, CILow, CIHigh
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, 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)
[Package vaxpmx version 0.0.3 Index]