Disaggregating proportional multistate lifetables by population heterogeneity to estimate intervention impacts on inequalities

Simulation models can be used to quantify the projected health impact of interventions. Quantifying heterogeneity in these impacts, for example by socioeconomic status, is important to understand impacts on health inequalities. We aim to disaggregate one type of Markov macro-simulation model, the proportional multistate lifetable, ensuring that under business-as-usual (BAU) the sum of deaths across disaggregated strata in each time step returns the same as the initial non-disaggregated model. We developed a disaggregation algorithm that iteratively rescales mortality, incidence and case fatality rates by time-step of the model to ensure correct total population counts were retained at each step. We then demonstrate the application by deprivation quintiles for New Zealand (NZ), for: hypothetical interventions (50% lower all-cause mortality, 50% lower coronary heart disease mortality) and a dietary intervention to substitute 59% of sodium with potassium chloride in the food supply. We found that the intervention scenario life years and health adjusted life years (HALYs) gained differed slightly when summed over the deprivation quintile compared to the aggregated model, due to the stratified model (appropriately) allowing for differential background mortality rates by strata. Modest differences in health gains (health adjusted life years) resulted from rescaling of sub-population mortality and incidence rates to ensure consistency with the aggregate population.