Childhood stunting in Indonesia: assessing the performance of Bayesian spatial conditional autoregressive models

Stunting continues to be a significant health issue, particularly in developing nations, with Indonesia ranking third in prevalence in Southeast Asia. This research examined the risk of stunting and influencing factors in Indonesia by implementing various Bayesian spatial conditional autoregressive (CAR) models that include covariates. A total of 750 models were run, including five different Bayesian spatial CAR models (Besag-York-Mollie (BYM), CAR Leroux and three forms of localised CAR), with 30 covariate combinations and five different hyperprior combinations for each model. The Poisson distribution was employed to model the counts of stunting cases. After a comprehensive evaluation of all model selection criteria utilized, the Bayesian localised CAR model with three covariates were preferred, either allowing up to 2 clusters with a variance hyperprior of inverse-gamma (1, 0.1) or allowing 3 clusters with a variance hyperprior of inverse-gamma (1, 0.01). Poverty and recent low birth weight (LBW) births are significantly associated with an increased risk of stunting, whereas child diet diversity is inversely related to the risk of stunting. Model results indicated that Sulawesi Barat Province has the highest risk of stunting, with DKI Jakarta Province the lowest. These areas with high stunting require interventions to reduce poverty, LBW births and increase child diet diversity.