Approximation of the Cox survival regression model by MCMC Bayesian Hierarchical Poisson modelling of factors associated with childhood mortality in Nigeria

The need for more pragmatic approaches to achieve sustainable development goal on childhood mortality reduction necessitated this study. Simultaneous study of the influence of where the children live and the censoring nature of children survival data is scarce. We identified the compositional and contextual factors associated with under-five (U5M) and infant (INM) mortality in Nigeria from 5 MCMC Bayesian hierarchical Poisson regression models as approximations of the Cox survival regression model. The 2018 DHS data of 33,924 under-five children were used. Life table techniques and the Mlwin 3.05 module for the analysis of hierarchical data were implemented in Stata Version 16. The overall INM rate (INMR) was 70 per 1000 livebirths compared with U5M rate (U5MR) of 131 per 1000 livebirth. The INMR was lowest in Ogun (17 per 1000 live births) and highest in Kaduna (106), Gombe (112) and Kebbi (116) while the lowest U5MR was found in Ogun (29) and highest in Jigawa (212) and Kebbi (248). The risks of INM and U5M were highest among children with none/low maternal education, multiple births, low birthweight, short birth interval, poorer households, when spouses decide on healthcare access, having a big problem getting to a healthcare facility, high community illiteracy level, and from states with a high proportion of the rural population in the fully adjusted model. Compared with the null model, 81% vs 13% and 59% vs 35% of the total variation in INM and U5M were explained by the state- and neighbourhood-level factors respectively. Infant- and under-five mortality in Nigeria is influenced by compositional and contextual factors. The Bayesian hierarchical Poisson regression model used in estimating the factors associated with childhood deaths in Nigeria fitted the survival data.