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.

Multilevel model with random effects for clustered survival data with multiple failure outcomes

We present a multilevel frailty model for handling serial dependence and simultaneous heterogeneity in survival data with a multilevel structure attributed to clustering of subjects and the presence of multiple failure outcomes. One commonly observes such data, for example, in multi-institutional, randomized placebo-controlled trials in which patients suffer repeated episodes (eg, recurrent migraines) of the disease outcome being measured. The model extends the proportional hazards model by incorporating a random covariate and unobservable random institution effect to respectively account for treatment-by-institution interaction and institutional variation in the baseline risk. Moreover, a random effect term with correlation structure driven by a first-order autoregressive process is attached to the model to facilitate estimation of between patient heterogeneity and serial dependence. By means of the generalized linear mixed model methodology, the random effects distribution is assumed normal and the residual maximum likelihood and the maximum likelihood methods are extended for estimation of model parameters. Simulation studies are carried out to evaluate the performance of the residual maximum likelihood and the maximum likelihood estimators and to assess the impact of misspecifying random effects distribution on the proposed inference. We demonstrate the practical feasibility of the modeling methodology by analyzing real data from a double-blind randomized multi-institutional clinical trial, designed to examine the effect of rhDNase on the occurrence of respiratory exacerbations among patients with cystic fibrosis.