Geographic inequalities in health intervention coverage – mapping the composite coverage index in Peru using geospatial modelling

Background: The composite coverage index (CCI) provides an integrated perspective towards universal health coverage in the context of reproductive, maternal, newborn and child health. Given the sample design of most household surveys does not provide coverage estimates below the first administrative level, approaches for achieving more granular estimates are needed. We used a model-based geostatistical approach to estimate the CCI at multiple resolutions in Peru. Methods: We generated estimates for the eight indicators on which the CCI is based for the departments, provinces, and areas of 5 × 5 km of Peru using data from two national household surveys carried out in 2018 and 2019 plus geospatial covariates. Bayesian geostatistical models were fit using the INLA-SPDE approach. We assessed model fit using cross-validation at the survey cluster level and by comparing modelled and direct survey estimates at the department-level. Results: CCI coverage in the provinces along the coast was consistently higher than in the remainder of the country. Jungle areas in the north and east presented the lowest coverage levels and the largest gaps between and within provinces. The greatest inequalities were found, unsurprisingly, in the largest provinces where populations are scattered in jungle territory and are difficult to reach. Conclusions: Our study highlighted provinces with high levels of inequality in CCI coverage indicating areas, mostly low-populated jungle areas, where more attention is needed. We also uncovered other areas, such as the border with Bolivia, where coverage is lower than the coastal provinces and should receive increased efforts. More generally, our results make the case for high-resolution estimates to unveil geographic inequities otherwise hidden by the usual levels of survey representativeness.