# \donttest{
## using the data (as output by crossVA function in openVA package for InSilicoVA algorithm)
data(comsamoz_public_openVAout)
head(comsamoz_public_openVAout$data) # head of the data
comsamoz_public_openVAout$data[1,] # ID and specific cause of death for individual 1
## mapped to broad cause
### same as comsamoz_public_broad$data
comsamoz_public_asbroad = cause_map(df = comsamoz_public_openVAout$data, age_group = "neonate")
head(comsamoz_public_asbroad)
### store broad cause map of the data
data(comsamoz_public_broad)
head(comsamoz_public_broad$data) # identical to head(comsamoz_public_asbroad)
## mapped to national death counts
comsamoz_public_asdeathcount = colSums(comsamoz_public_asbroad)
## VA-calibration for the "neonate" age group and InSilicoVA algorithm
## input as specific cause
calib_out_asspecific = vacalibration(va_data = setNames(list(comsamoz_public_openVAout$data),
list(comsamoz_public_openVAout$va_algo)),
age_group = comsamoz_public_openVAout$age_group,
country = "Mozambique")
## input as broad cause
calib_out_asbroad = vacalibration(va_data = setNames(list(comsamoz_public_asbroad),
list(comsamoz_public_openVAout$va_algo)),
age_group = comsamoz_public_openVAout$age_group,
country = "Mozambique")
## input as specific cause
calib_out_asdeathcount = vacalibration(va_data = setNames(list(comsamoz_public_asdeathcount),
list(comsamoz_public_openVAout$va_algo)),
age_group = comsamoz_public_openVAout$age_group,
country = "Mozambique")
## comparing uncalibrated CSMF estimates and posterior summary of calibrated CSMF estimates
calib_out_asspecific$p_uncalib
calib_out_asspecific$pcalib_postsumm[1,,]
calib_out_asbroad$p_uncalib
calib_out_asbroad$pcalib_postsumm[1,,]
calib_out_asdeathcount$p_uncalib
calib_out_asdeathcount$pcalib_postsumm[1,,]
# }
Run the code above in your browser using DataLab