The goal of vaxpmx is to provide the functions for pharmacometricmodeling in vaccines, specifically analyses related to correlates ofprotection (CoPs) and vaccine efficacy/effectiveness predictions using aCoP.
You can install the released version of vaxpmx fromCRAN with:
install.packages("vaxpmx")This is a basic example which shows you how to solve a commonproblem:
library(vaxpmx)library(survival)# Load an example datasetdata(data_temp)# Fit logistic model relating neutralizing titer to disease statuslogisticFit<-glm(disease_any~ nAb1,data = data_temp,family =binomial())# Fit Cox proportional hazards model relating neutralizing titer to time to disease or end of follow-up# coxFit <- coxph(Surv(time_event, disease_any) ~ nAb1, data = data_temp)# Estimate vaccine efficacy and 95\% confidence interval based on the fitted modelsve(logisticFit, data_temp,nboot =500)ve(coxFit, data_temp,nboot =500)