# \donttest{
## using the data
data(comsamoz_public_broad)
head(comsamoz_public_broad$data) # head of the data
comsamoz_public_broad$data[1,] # binary vector indicating cause of death for individual 1
## mapped to national death counts
comsamoz_public_asdeathcount = colSums(comsamoz_public_broad$data)
## VA-calibration for the "neonate" age group and InSilicoVA algorithm
## input as broad cause
calib_out_asbroad = vacalibration(va_data = setNames(list(comsamoz_public_broad$data),
list(comsamoz_public_broad$va_algo)),
age_group = comsamoz_public_broad$age_group,
country = "Mozambique")
## input as specific cause
calib_out_asdeathcount = vacalibration(va_data = setNames(list(comsamoz_public_asdeathcount),
list(comsamoz_public_broad$va_algo)),
age_group = comsamoz_public_broad$age_group,
country = "Mozambique")
## comparing uncalibrated CSMF estimates and posterior summary of calibrated CSMF estimates
## all are the same
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