Abstract
Engaging the private sector in contraceptive method supply is critical for creating equitable, sustainable, and accessible healthcare systems. To achieve this, it is essential to understand where women obtain their modern contraceptives. While national-level estimates provide valuable insights into overall trends in contraceptive supply, they often obscure variation within and across subnational regions. Addressing localised needs has become increasingly important as countries adopt decentralised models for family planning services. Decentralization has also underscored the need for reliable subnational estimates of key family planning indicators. The absence of regularly collected subnational data has hindered effective monitoring and decision-making. To bridge this gap, we propose a novel approach that leverages latent attributes in Demographic and Health Survey (DHS) data to produce Bayesian probabilistic projections of contraceptive method supply shares (the proportions of modern contraceptive methods supplied by public and private sectors) with limited data. Our modelling framework is built on Bayesian hierarchical models. Using penalized splines to track public and private supply shares over time, we leverage the spatial nature of the data and incorporate a correlation structure between recent supply share observations at national and subnational levels. This framework contributes to the domain of subnational estimation of proportions in data-sparse settings, outperforming comparable and previous approaches. As decentralization continues to reshape family planning services, producing reliable subnational estimates of key indicators is increasingly vital for researchers and policymakers.