Spatio-temporal modelling of routine health facility data for malaria risk micro-stratification in mainland Tanzania.

As malaria transmission declines, the need to monitor the heterogeneity of malaria risk at finer scales becomes critical to guide community-based targeted interventions. Although routine health facility (HF) data can provide epidemiological evidence at high spatial and temporal resolution, its incomplete nature of information can result in lower administrative units without empirical data. To overcome geographic sparsity of data and its representativeness, geo-spatial models can leverage routine information to predict risk in un-represented areas as well as estimate uncertainty of predictions. Here, a Bayesian spatio-temporal model was applied on malaria test positivity rate (TPR) data for the period 2017-2019 to predict risks at the ward level, the lowest decision-making unit in mainland Tanzania. To quantify the associated uncertainty, the probability of malaria TPR exceeding programmatic threshold was estimated. Results showed a marked spatial heterogeneity in malaria TPR across wards. 17.7 million people resided in areas where malaria TPR was high (≥ 30; 90% certainty) in the North-West and South-East parts of Tanzania. Approximately 11.7 million people lived in areas where malaria TPR was very low (< 5%; 90% certainty). HF data can be used to identify different epidemiological strata and guide malaria interventions at micro-planning units in Tanzania. These data, however, are imperfect in many settings in Africa and often require application of geo-spatial modelling techniques for estimation.

Geographic accessibility and hospital competition for emergency blood transfusion services in Bungoma, Western Kenya.

BACKGROUND: Estimating accessibility gaps to essential health interventions helps to allocate and prioritize health resources. Access to blood transfusion represents an important emergency health requirement. Here, we develop geo-spatial models of accessibility and competition to blood transfusion services in Bungoma County, Western Kenya. METHODS: Hospitals providing blood transfusion services in Bungoma were identified from an up-dated geo-coded facility database. AccessMod was used to define care-seeker's travel times to the nearest blood transfusion service. A spatial accessibility index for each enumeration area (EA) was defined using modelled travel time, population demand, and supply available at the hospital, assuming a uniform risk of emergency occurrence in the county. To identify populations marginalized from transfusion services, the number of people outside 1-h travel time and those residing in EAs with low accessibility indexes were computed at the sub-county level. Competition between the transfusing hospitals was estimated using a spatial competition index which provided a measure of the level of attractiveness of each hospital. To understand whether highly competitive facilities had better capacity for blood transfusion services, a correlation test between the computed competition metric and the blood units received and transfused at the hospital was done. RESULTS: 15 hospitals in Bungoma county provide transfusion services, however these are unevenly distributed across the sub-counties. Average travel time to a blood transfusion centre in the county was 33 min and 5% of the population resided outside 1-h travel time. Based on the accessibility index, 38% of the EAs were classified to have low accessibility, representing 34% of the population, with one sub-county having the highest marginalized population. The computed competition index showed that hospitals in the urban areas had a spatial competitive advantage over those in rural areas. CONCLUSION: The modelled spatial accessibility has provided an improved understanding of health care gaps essential for health planning. Hospital competition has been illustrated to have some degree of influence in provision of health services hence should be considered as a significant external factor impacting the delivery, and re-design of available services.

Quality of maternal healthcare and travel time influence birthing service utilisation in Ghanaian health facilities: a geographical analysis of routine health data.

OBJECTIVES: To investigate how the quality of maternal health services and travel times to health facilities affect birthing service utilisation in Eastern Region, Ghana. DESIGN: The study is a cross-sectional spatial interaction analysis of birth service utilisation patterns. Routine birth data were spatially linked to quality care, service demand and travel time data. SETTING: 131 Health facilities (public, private and faith-based) in 33 districts in Eastern Region, Ghana. PARTICIPANTS: Women who gave birth in health facilities in the Eastern Region, Ghana in 2017. OUTCOME MEASURES: The count of women giving birth, the quality of birthing care services and the geographic coverage of birthing care services. RESULTS: As travel time from women's place of residence to the health facility increased up to two2 hours, the utilisation rate markedly decreased. Higher quality of maternal health services haves a larger, positive effect on utilisation rates than service proximity. The quality of maternal health services was higher in hospitals than in primary care facilities. Most women (88.6%) travelling via mechanised transport were within two2 hours of any birthing service. The majority (56.2%) of women were beyond the two2 -hour threshold of critical comprehensive emergency obstetric and newborn care (CEmONC) services. Few CEmONC services were in urban centres, disadvantaging rural populations. CONCLUSIONS: To increase birthing service utilisation in Ghana, higher quality health facilities should be located closer to women, particularly in rural areas. Beyond Ghana, routinely collected birth records could be used to understand the interaction of service proximity and quality.

Unmet need for COVID-19 vaccination coverage in Kenya.

COVID-19 has impacted the health and livelihoods of billions of people since it emerged in 2019. Vaccination for COVID-19 is a critical intervention that is being rolled out globally to end the pandemic. Understanding the spatial inequalities in vaccination coverage and access to vaccination centres is important for planning this intervention nationally. Here, COVID-19 vaccination data, representing the number of people given at least one dose of vaccine, a list of the approved vaccination sites, population data and ancillary GIS data were used to assess vaccination coverage, using Kenya as an example. Firstly, physical access was modelled using travel time to estimate the proportion of population within 1 hour of a vaccination site. Secondly, a Bayesian conditional autoregressive (CAR) model was used to estimate the COVID-19 vaccination coverage and the same framework used to forecast coverage rates for the first quarter of 2022. Nationally, the average travel time to a designated COVID-19 vaccination site (n = 622) was 75.5 min (Range: 62.9 - 94.5 min) and over 87% of the population >18 years reside within 1 hour to a vaccination site. The COVID-19 vaccination coverage in December 2021 was 16.70% (95% CI: 16.66 - 16.74) - 4.4 million people and was forecasted to be 30.75% (95% CI: 25.04 - 36.96) - 8.1 million people by the end of March 2022. Approximately 21 million adults were still unvaccinated in December 2021 and, in the absence of accelerated vaccine uptake, over 17.2 million adults may not be vaccinated by end March 2022 nationally. Our results highlight geographic inequalities at sub-national level and are important in targeting and improving vaccination coverage in hard-to-reach populations. Similar mapping efforts could help other countries identify and increase vaccination coverage for such populations.

Understanding factors associated with attending secondary school in Tanzania using household survey data.

BACKGROUND: Sustainable Development Goal (SDG) 4 aims to ensure inclusive and equitable access for all by 2030, leaving no one behind. One indicator selected to measure progress towards achievement is the participation rate of youth in education (SDG 4.3.1). Here we aim to understand drivers of school attendance using one country in East Africa as an example. METHODS: Nationally representative household survey data (2015-16 Tanzania Demographic and Health Survey) were used to explore individual, household and contextual factors associated with secondary school attendance in Tanzania. These included, age, head of household's levels of education, gender, household wealth index and total number of children under five. Contextual factors such as average pupil to qualified teacher ratio and geographic access to school were also tested at cluster level. A two-level random intercept logistic regression model was used in exploring association of these factors with attendance in a multi-level framework. RESULTS: Age of household head, educational attainments of either of the head of the household or parent, child characteristics such as gender, were important predictors of secondary school attendance. Being in a richer household and with fewer siblings of lower age (under the age of 5) were associated with increased odds of attendance (OR = 0.91, CI 95%: 0.86; 0.96). Contextual factors were less likely to be associated with secondary school attendance. CONCLUSIONS: Individual and household level factors are likely to impact secondary school attendance rates more compared to contextual factors, suggesting an increased focus of interventions at these levels is needed. Future studies should explore the impact of interventions targeting these levels. Policies should ideally promote gender equality in accessing secondary school as well as support those families where the dependency ratio is high. Strategies to reduce poverty will also increase the likelihood of attending school.

Malaria hospitalisation in East Africa: age, phenotype and transmission intensity.

BACKGROUND: Understanding the age patterns of disease is necessary to target interventions to maximise cost-effective impact. New malaria chemoprevention and vaccine initiatives target young children attending routine immunisation services. Here we explore the relationships between age and severity of malaria hospitalisation versus malaria transmission intensity. METHODS: Clinical data from 21 surveillance hospitals in East Africa were reviewed. Malaria admissions aged 1 month to 14 years from discrete administrative areas since 2006 were identified. Each site-time period was matched to a model estimated community-based age-corrected parasite prevalence to provide predictions of prevalence in childhood (PfPR2-10). Admission with all-cause malaria, severe malaria anaemia (SMA), respiratory distress (RD) and cerebral malaria (CM) were analysed as means and predicted probabilities from Bayesian generalised mixed models. RESULTS: 52,684 malaria admissions aged 1 month to 14 years were described at 21 hospitals from 49 site-time locations where PfPR2-10 varied from < 1 to 48.7%. Twelve site-time periods were described as low transmission (PfPR2-10 < 5%), five low-moderate transmission (PfPR2-10 5-9%), 20 moderate transmission (PfPR2-10 10-29%) and 12 high transmission (PfPR2-10 ≥ 30%). The majority of malaria admissions were below 5 years of age (69-85%) and rare among children aged 10-14 years (0.7-5.4%) across all transmission settings. The mean age of all-cause malaria hospitalisation was 49.5 months (95% CI 45.1, 55.4) under low transmission compared with 34.1 months (95% CI 30.4, 38.3) at high transmission, with similar trends for each severe malaria phenotype. CM presented among older children at a mean of 48.7 months compared with 39.0 months and 33.7 months for SMA and RD, respectively. In moderate and high transmission settings, 34% and 42% of the children were aged between 2 and 23 months and so within the age range targeted by chemoprevention or vaccines. CONCLUSIONS: Targeting chemoprevention or vaccination programmes to areas where community-based parasite prevalence is ≥10% is likely to match the age ranges covered by interventions (e.g. intermittent presumptive treatment in infancy to children aged 2-23 months and current vaccine age eligibility and duration of efficacy) and the age ranges of highest disease burden.

Insecticide-treated net distribution in Western Kenya: impacts related to COVID-19 and health worker strikes.

We examined the impact of coronavirus disease (COVID) mitigation, supply and distribution interruptions on the delivery of long-lasting insecticide-treated nets (LLINs) in Western Kenya. The median monthly distribution of LLINs declined during COVID mitigation strategies (March-July 2020) and during the health worker strikes (December 2020-February 2021). Recovery periods followed initial declines, indicative of a 'catching up' on missed routine distribution. Mass community campaigns were delayed by 10 months. These observations offer encouragement for routine net distribution resilience, but complete interruptions of planned mass distributions require alternate strategies during pandemics.

Model building and assessment of the impact of covariates for disease prevalence mapping in low-resource settings: to explain and to predict.

This paper provides statistical guidance on the development and application of model-based geostatistical methods for disease prevalence mapping. We illustrate the different stages of the analysis, from exploratory analysis to spatial prediction of prevalence, through a case study on malaria mapping in Tanzania. Throughout the paper, we distinguish between predictive modelling, whose main focus is on maximizing the predictive accuracy of the model, and explanatory modelling, where greater emphasis is placed on understanding the relationships between the health outcome and risk factors. We demonstrate that these two paradigms can result in different modelling choices. We also propose a simple approach for detecting over-fitting based on inspection of the correlation matrix of the estimators of the regression coefficients. To enhance the interpretability of geostatistical models, we introduce the concept of domain effects in order to assist variable selection and model validation. The statistical ideas and principles illustrated here in the specific context of disease prevalence mapping are more widely applicable to any regression model for the analysis of epidemiological outcomes but are particularly relevant to geostatistical models, for which the separation between fixed and random effects can be ambiguous.

Malaria micro-stratification using routine surveillance data in Western Kenya.

BACKGROUND: There is an increasing need for finer spatial resolution data on malaria risk to provide micro-stratification to guide sub-national strategic plans. Here, spatial-statistical techniques are used to exploit routine data to depict sub-national heterogeneities in test positivity rate (TPR) for malaria among patients attending health facilities in Kenya. METHODS: Routine data from health facilities (n = 1804) representing all ages over 24 months (2018-2019) were assembled across 8 counties (62 sub-counties) in Western Kenya. Statistical model-based approaches were used to quantify heterogeneities in TPR and uncertainty at fine spatial resolution adjusting for missingness, population distribution, spatial data structure, month, and type of health facility. RESULTS: The overall monthly reporting rate was 78.7% (IQR 75.0-100.0) and public-based health facilities were more likely than private facilities to report ≥ 12 months (OR 5.7, 95% CI 4.3-7.5). There was marked heterogeneity in population-weighted TPR with sub-counties in the north of the lake-endemic region exhibiting the highest rates (exceedance probability > 70% with 90% certainty) where approximately 2.7 million (28.5%) people reside. At micro-level the lowest rates were in 14 sub-counties (exceedance probability < 30% with 90% certainty) where approximately 2.2 million (23.1%) people lived and indoor residual spraying had been conducted since 2017. CONCLUSION: The value of routine health data on TPR can be enhanced when adjusting for underlying population and spatial structures of the data, highlighting small-scale heterogeneities in malaria risk often masked in broad national stratifications. Future research should aim at relating these heterogeneities in TPR with traditional community-level prevalence to improve tailoring malaria control activities at sub-national levels.

Plasmodium falciparum parasite prevalence in East Africa: Updating data for malaria stratification.

The High Burden High Impact (HBHI) strategy for malaria encourages countries to use multiple sources of available data to define the sub-national vulnerabilities to malaria risk, including parasite prevalence. Here, a modelled estimate of Plasmodium falciparum from an updated assembly of community parasite survey data in Kenya, mainland Tanzania, and Uganda is presented and used to provide a more contemporary understanding of the sub-national malaria prevalence stratification across the sub-region for 2019. Malaria prevalence data from surveys undertaken between January 2010 and June 2020 were assembled form each of the three countries. Bayesian spatiotemporal model-based approaches were used to interpolate space-time data at fine spatial resolution adjusting for population, environmental and ecological covariates across the three countries. A total of 18,940 time-space age-standardised and microscopy-converted surveys were assembled of which 14,170 (74.8%) were identified after 2017. The estimated national population-adjusted posterior mean parasite prevalence was 4.7% (95% Bayesian Credible Interval 2.6-36.9) in Kenya, 10.6% (3.4-39.2) in mainland Tanzania, and 9.5% (4.0-48.3) in Uganda. In 2019, more than 12.7 million people resided in communities where parasite prevalence was predicted ≥ 30%, including 6.4%, 12.1% and 6.3% of Kenya, mainland Tanzania and Uganda populations, respectively. Conversely, areas that supported very low parasite prevalence (<1%) were inhabited by approximately 46.2 million people across the sub-region, or 52.2%, 26.7% and 10.4% of Kenya, mainland Tanzania and Uganda populations, respectively. In conclusion, parasite prevalence represents one of several data metrics for disease stratification at national and sub-national levels. To increase the use of this metric for decision making, there is a need to integrate other data layers on mortality related to malaria, malaria vector composition, insecticide resistance and bionomic, malaria care-seeking behaviour and current levels of unmet need of malaria interventions.

Routine data for malaria morbidity estimation in Africa: challenges and prospects.

BACKGROUND: The burden of malaria in sub-Saharan Africa remains challenging to measure relying on epidemiological modelling to evaluate the impact of investments and providing an in-depth analysis of progress and trends in malaria response globally. In malaria-endemic countries of Africa, there is increasing use of routine surveillance data to define national strategic targets, estimate malaria case burdens and measure control progress to identify financing priorities. Existing research focuses mainly on the strengths of these data with less emphasis on existing challenges and opportunities presented. CONCLUSION: Here we define the current imperfections common to routine malaria morbidity data at national levels and offer prospects into their future use to reflect changing disease burdens.

Severe-malaria infection and its outcomes among pregnant women in Burkina Faso health-districts: Hierarchical Bayesian space-time models applied to routinely-collected data from 2013 to 2018.

Fine-scale hotspots detection is crucial for optimum delivery of essential health-services for reducing severe malaria in pregnancy (MiP) and death cases in Burkina Faso. This study used hierarchical-Bayesian Spatio-temporal modeling to explore space-time patterns and pinpoint health-districts with an exceedance probability of severe MiP incidence and fatality rate. Study also assessed effect of health-district service delivery (readiness) on severe-MiP outcomes. Severe-MiP fatality rate declined considerably while its incidence rate remained unchanged between January-2013 and December-2018. Severe-MiP cases persisted throughout the year with peaks between August and November. These peaks increased 2.5-fold the fatality rate. Furthermore, severe-MiP fatality was higher in health-districts classified as low-readiness (IRR = 2.469, 95%CrI: 1.632-3.738). However, the fatality rate decreased significantly with proper coverage with three doses for intermittent-preventive-treatment with sulphadoxine-pyrimethamine. Severe-MiP burden was heterogeneous spatially and temporally. The study suggested that health-programs should increase health-districts readiness and optimize resource allocation in high burden areas and months.

Bayesian Spatiotemporal Modeling of Routinely Collected Data to Assess the Effect of Health Programs in Malaria Incidence During Pregnancy in Burkina Faso.

Control of malaria in pregnancy (MiP) remains a major challenge in Burkina Faso. Surveillance of the burden due to MiP based on routinely collected data at a fine-scale level, followed by an appropriate analysis and interpretation, may be crucial for evaluating and improving the effectiveness of existing control measures. We described the spatio-temporal dynamics of MiP at the community-level and assessed health program effects, mainly community-based health promotion, results-based financing, and intermittent-preventive-treatment with sulphadoxine-pyrimethamine (IPTp-SP). Community-aggregated monthly MiP cases were downloaded from Health Management Information System and combined with covariates from other sources. The MiP spatio-temporal pattern was decomposed into three components: overall spatial and temporal trends and space-time interaction. Bayesian hierarchical spatio-temporal Poisson models were used to fit the MiP incidence rate and assess health program effects. The overall annual incidence increased between 2015 and 2017. The findings reveal spatio-temporal heterogenicity throughout the year, which peaked during rainy season. From the model without covariates, 96 communities located mainly in the Cascades, South-West, Center-West, Center-East, and Eastern regions, exhibited significant relative-risk levels. The combined effect (significant reducing effect) of RBF, health promotion and IPTp-SP strategies was greatest in 17.7% (17/96) of high burden malaria communities. Despite intensification of control efforts, MiP remains high at the community-scale. The provided risk maps are useful tools for highlighting areas where interventions should be optimized, particularly in high-risk communities.

Estimating hospital catchments from in-patient admission records: a spatial statistical approach applied to malaria.

Admission records are seldom used in sub-Saharan Africa to delineate hospital catchments for the spatial description of hospitalised disease events. We set out to investigate spatial hospital accessibility for severe malarial anaemia (SMA) and cerebral malaria (CM). Malaria admissions for children between 1 month and 14 years old were identified from prospective clinical surveillance data recorded routinely at four referral hospitals covering two complete years between December 2015 to November 2016 and November 2017 to October 2018. These were linked to census enumeration areas (EAs) with an age-structured population. A novel mathematical-statistical framework that included EAs with zero observations was used to predict hospital catchment for malaria admissions adjusting for spatial distance. From 5766 malaria admissions, 5486 (95.14%) were linked to specific EA address, of which 272 (5%) were classified as cerebral malaria while 1001 (10%) were severe malaria anaemia. Further, results suggest a marked geographic catchment of malaria admission around the four sentinel hospitals although the extent varied. The relative rate-ratio of hospitalisation was highest at <1-hour travel time for SMA and CM although this was lower outside the predicted hospital catchments. Delineation of catchments is important for planning emergency care delivery and in the use of hospital data to define epidemiological disease burdens. Further hospital and community-based studies on treatment-seeking pathways to hospitals for severe disease would improve our understanding of catchments.

A spatial database of health facilities managed by the public health sector in sub Saharan Africa.

Health facilities form a central component of health systems, providing curative and preventative services and structured to allow referral through a pyramid of increasingly complex service provision. Access to health care is a complex and multidimensional concept, however, in its most narrow sense, it refers to geographic availability. Linking health facilities to populations has been a traditional per capita index of heath care coverage, however, with locations of health facilities and higher resolution population data, Geographic Information Systems allow for a more refined metric of health access, define geographic inequalities in service provision and inform planning. Maximizing the value of spatial heath access requires a complete census of providers and their locations. To-date there has not been a single, geo-referenced and comprehensive public health facility database for sub-Saharan Africa. We have assembled national master health facility lists from a variety of government and non-government sources from 50 countries and islands in sub Saharan Africa and used multiple geocoding methods to provide a comprehensive spatial inventory of 98,745 public health facilities.

Using parasite genetic and human mobility data to infer local and cross-border malaria connectivity in Southern Africa.

Local and cross-border importation remain major challenges to malaria elimination and are difficult to measure using traditional surveillance data. To address this challenge, we systematically collected parasite genetic data and travel history from thousands of malaria cases across northeastern Namibia and estimated human mobility from mobile phone data. We observed strong fine-scale spatial structure in local parasite populations, providing positive evidence that the majority of cases were due to local transmission. This result was largely consistent with estimates from mobile phone and travel history data. However, genetic data identified more detailed and extensive evidence of parasite connectivity over hundreds of kilometers than the other data, within Namibia and across the Angolan and Zambian borders. Our results provide a framework for incorporating genetic data into malaria surveillance and provide evidence that both strengthening of local interventions and regional coordination are likely necessary to eliminate malaria in this region of Southern Africa.

Mapping vaccination coverage to explore the effects of delivery mechanisms and inform vaccination strategies.

The success of vaccination programs depends largely on the mechanisms used in vaccine delivery. National immunization programs offer childhood vaccines through fixed and outreach services within the health system and often, additional supplementary immunization activities (SIAs) are undertaken to fill gaps and boost coverage. Here, we map predicted coverage at 1 × 1 km spatial resolution in five low- and middle-income countries to identify areas that are under-vaccinated via each delivery method using Demographic and Health Surveys data. We compare estimates of the coverage of the third dose of diphtheria-tetanus-pertussis-containing vaccine (DTP3), which is typically delivered through routine immunization (RI), with those of measles-containing vaccine (MCV) for which SIAs are also undertaken. We find that SIAs have boosted MCV coverage in some places, but not in others, particularly where RI had been deficient, as depicted by DTP coverage. The modelling approaches outlined here can help to guide geographical prioritization and strategy design.

Spatial distribution and determinants of asymptomatic malaria risk among children under 5 years in 24 districts in Burkina Faso.

BACKGROUND: In malaria endemic countries, asymptomatic cases constitute an important reservoir of infections sustaining transmission. Estimating the burden of the asymptomatic population and identifying areas with elevated risk is important for malaria control in Burkina Faso. This study analysed the spatial distribution of asymptomatic malaria infection among children under 5 in 24 health districts in Burkina Faso and identified the determinants of this distribution. METHODS: The data used in this study were collected in a baseline survey on "evaluation of the impact of pay for performance on the quality of care" conducted in 24 health districts in Burkina Faso, between October 2013 and March 2014. This survey involved 7844 households and 1387 community health workers. A Bayesian hierarchical logistic model that included spatial dependence and covariates was implemented to identify the determinants of asymptomatic malaria infection. The posterior probability distribution of a parameter from the model was summarized using odds ratio (OR) and 95% credible interval (95% CI). RESULTS: The overall prevalence of asymptomatic malaria infection in children under 5 years of age was estimated at 38.2%. However, significant variation was observed between districts ranging from 11.1% in the district of Barsalgho to 77.8% in the district of Gaoua. Older children (48-59 vs < 6 months: OR: 6.79 [5.62, 8.22]), children from very poor households (Richest vs poorest: OR: 0.85 [0.74-0.96]), households located more than 5 km from a health facility (< 5 km vs  ≥ 5 km: OR: 1.14 [1.04-1.25]), in localities with inadequate number of nurses (< 3 vs ≥ 3: 0.72 [0.62, 0.82], from rural areas (OR: 1.67 [1.39-2.01]) and those surveyed in high transmission period of asymptomatic malaria (OR: 1.27 [1.10-1.46]) were most at risk for asymptomatic malaria infection. In addition, the spatial analysis identified the following nine districts that reported significantly higher risks: Batié, Boromo, Dano, Diébougou, Gaoua, Ouahigouya, Ouargaye, Sapouy and Toma. The district of Zabré reported the lowest risk. CONCLUSION: The analysis of spatial distribution of infectious reservoir allowed the identification of risk areas as well as the identification of individual and contextual factors. Such national spatial analysis should help to prioritize areas for increased malaria control activities.

Temporal trends in spatial inequalities of maternal and newborn health services among four east African countries, 1999-2015.

BACKGROUND: Sub-Saharan Africa continues to account for the highest regional maternal mortality ratio (MMR) in the world, at just under 550 maternal deaths per 100,000 live births in 2015, compared to a global rate of 216 deaths. Spatial inequalities in access to life-saving maternal and newborn health (MNH) services persist within sub-Saharan Africa, however, with varied improvement over the past two decades. While previous research within the East African Community (EAC) region has examined utilisation of MNH care as an emergent property of geographic accessibility, no research has examined how these spatial inequalities have evolved over time at similar spatial scales. METHODS: Here, we analysed temporal trends of spatial inequalities in utilisation of antenatal care (ANC), skilled birth attendance (SBA), and postnatal care (PNC) among four East African countries. Specifically, we used Bayesian spatial statistics to generate district-level estimates of these services for several time points using Demographic and Health Surveys data in Kenya, Tanzania, Rwanda, and Uganda. We examined temporal trends of both absolute and relative indices over time, including the absolute difference between estimates, as well as change in performance ratios of the best-to-worst performing districts per country. RESULTS: Across all countries, we found the greatest spatial equality in ANC, while SBA and PNC tended to have greater spatial variability. In particular, Rwanda represented the only country to consistently increase coverage and reduce spatial inequalities across all services. Conversely, Tanzania had noticeable reductions in ANC coverage throughout most of the country, with some areas experiencing as much as a 55% reduction. Encouragingly, however, we found that performance gaps between districts have generally decreased or remained stably low across all countries, suggesting countries are making improvements to reduce spatial inequalities in these services. CONCLUSIONS: We found that while the region is generally making progress in reducing spatial gaps across districts, improvement in PNC coverage has stagnated, and should be monitored closely over the coming decades. This study is the first to report temporal trends in district-level estimates in MNH services across the EAC region, and these findings establish an important baseline of evidence for the Sustainable Development Goal era.

National and sub-national variation in patterns of febrile case management in sub-Saharan Africa.

Given national healthcare coverage gaps, understanding treatment-seeking behaviour for fever is crucial for the management of childhood illness and to reduce deaths. Here, we conduct a modelling study triangulating household survey data for fever in children under the age of five years with georeferenced public health facility databases (n = 86,442 facilities) in 29 countries across sub-Saharan Africa, to estimate the probability of seeking treatment for fever at public facilities. A Bayesian item response theory framework is used to estimate this probability based on reported fever episodes, treatment choice, residence, and estimated travel-time to the nearest public-sector health facility. Findings show inter- and intra-country variation, with the likelihood of seeking treatment for fever less than 50% in 16 countries. Results highlight the need to invest in public healthcare and related databases. The variation in public sector use illustrates the need to include such modelling in future infectious disease burden estimation.