if (FALSE) {
data(NeonatesVA5)
fit_insilico <- codeVA(NeonatesVA5)
insilico_prep <- prepCalibration(fit_insilico)
calib_insilico = vacalibration::vacalibration(va_data = insilico_prep,
age_group = "neonate",
country = "Mozambique")
summary(calib_insilico)
fit_interva <- codeVA(NeonatesVA5, model = "InterVA", version = "5", write = FALSE)
interva_prep <- prepCalibration(fit_interva)
calib_interva = vacalibration::vacalibration(va_data = interva_prep,
age_group = "neonate",
country = "Mozambique")
summary(calib_interva, top = 3, rnd = 2)
two_fits <- prepCalibration(fit_insilico, fit_interva)
calib_ensemble = vacalibration::vacalibration(va_data = two_fits,
age_group = "neonate",
country = "Mozambique")
summary(calib_ensemble, algorithm = "ensemble")
summary(calib_ensemble, algorithm = c("ensemble", "insilicova"))
}
Run the code above in your browser using DataLab