Do ethnic and socio-economic inequalities in mortality vary by region in New Zealand? An application of hierarchical Bayesian modelling

We hypothesised that ethnic and socio-economic inequality in mortality might vary by region in New Zealand. Linked 2001-2004 census-mortality data were stratified by region (District Health Boards or DHBs), sex, age and ethnic groups, and income quintiles. To accommodate data sparseness, and to achieve accurate estimates of DHB-specific mortality rates and rate ratios by ethnicity and income, we used hierarchical Bayesian methods. To aid presentation of results, we used posterior mortality rates from the models to calculate directly standardised rates and rate ratios, with credible intervals. Māori-European/Other mortality rate ratios were often similar across DHBs, but Waitemata and Canterbury DHBs (both predominantly urban areas with low Māori population) had significantly lower rate ratios. In contrast, Bay of Plenty and Waikato DHBs (heterogeneous by both ethnicity and socio-economic position) had significantly higher rate ratios. There was little variation in mortality inequalities by income across DHBs. Examining the underlying rates for ethnic and income groups separately, there were significant variations across DHBs, but these were often correlated such that the ethnic or income rate ratio was similar across DHBs. The application of hierarchical Bayesian allowed more definitive conclusions than routine empirical methods when comparing small populations such as social groups across regions. The range of hierarchical Bayesian estimates of Māori mortality and Māori:European rate ratios across regions was considerably narrower than empirical standardisation estimates.