Female genital mutilation/cutting (FGM/C) is considered a public health and human rights concern, mainly concentrated in Africa, and has been targeted for elimination under the sustainable development goals. Interventions aimed at ending the practice often rely on data from household surveys which employ complex designs leading to outcomes that are not totally independent, thus requiring advanced statistical techniques. Combining data from multiple surveys within robust statistical framework holds promise to provide more precise estimates due to increased sample size, and accurately identify ‘hotspots’ and allow for assessment of changes over time. In this study, rich datasets from six (6) successive waves of the Nigeria Demographic and Health Surveys and Multiple Indicator Cluster Surveys undertaken between 2003 and 2016/17, were combined and analyzed in order to better assess changes in the likelihood and prevalence of FGM/C among 0-14-year old girls in Nigeria. We used Bayesian hierarchical regression models which explicitly accounted for the inherent spatial and temporal autocorrelations within the data while simultaneously adjusting for variations due to different survey methods and the effects of linear and non-linear covariates. Parameters were estimated using Markov chain Mote Carlo techniques and model fit assessments were based on Deviance Information Criterion. Results show that prevalence of FGM/C among 0–14 years old girls in Nigeria varied over time and across geographical locations and peaked in 2008 with a shift from South to North. A girl was more likely to be cut if her mother was cut, supported FGM/C continuation, or had no higher education. The effects of mother’s age, wealth and type of residence (urban-rural) were no longer significant in 2016. These results reflect the gains of interventions over the years, but also echo the belief that FGM/C is a social norm thus requiring tailored all-inclusive interventions for the total abandonment of FGM/C in Nigeria.
