The global burden of falls: Global, regional and national estimates of morbidity and mortality from the Global Burden of Disease Study 2017

Background: Falls can lead to severe health loss including death. Past research has shown that falls are an important cause of death and disability worldwide. The Global Burden of Disease Study 2017 (GBD 2017) provides a comprehensive assessment of morbidity and mortality from falls. Methods: Estimates for mortality, years of life lost (YLLs), incidence, prevalence, years lived with disability (YLDs) and disability-adjusted life years (DALYs) were produced for 195 countries and territories from 1990 to 2017 for all ages using the GBD 2017 framework. Distributions of the bodily injury (eg, hip fracture) were estimated using hospital records. Results: Globally, the age-standardised incidence of falls was 2238 (1990-2532) per 100 000 in 2017, representing a decline of 3.7% (7.4 to 0.3) from 1990 to 2017. Age-standardised prevalence w

Mapping 123 million neonatal, infant and child deaths between 2000 and 2017

Since 2000, many countries have achieved considerable success in improving child survival, but localized progress remains unclear. To inform efforts towards United Nations Sustainable Development Goal 3.2—to end preventable child deaths by 2030—we need consistently estimated data at the subnational level regarding child mortality rates and trends. Here we quantified, for the period 2000–2017, the subnational variation in mortality rates and number of deaths of neonates, infants and children under 5 years of age within 99 low- and middle-income countries using a geostatistical survival model. We estimated that 32% of children under 5 in these countries lived in districts that had attained rates of 25 or fewer child deaths per 1,000 live births by 2017, and that 58% of child deaths between 2000 and 2017 in these countries could have been averted in the absence of geographical inequality. This study enables the identification of high-mortality clusters, patterns of progress and geographical inequalities to inform appropriate investments and implementations that will help to improve the health of all populations

Global age-sex-specific fertility, mortality, healthy life expectancy (HALE), and population estimates in 204 countries and territories, 1950-2019 : a comprehensive demographic analysis for the Global Burden of Disease Study 2019

Background: Accurate and up-to-date assessment of demographic metrics is crucial for understanding a wide range of social, economic, and public health issues that affect populations worldwide. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019 produced updated and comprehensive demographic assessments of the key indicators of fertility, mortality, migration, and population for 204 countries and territories and selected subnational locations from 1950 to 2019. Methods: 8078 country-years of vital registration and sample registration data, 938 surveys, 349 censuses, and 238 other sources were identified and used to estimate age-specific fertility. Spatiotemporal Gaussian process regression (ST-GPR) was used to generate age-specific fertility rates for 5-year age groups between ages 15 and 49 years. With extensions to age groups 10–14 and 50–54 years, the total fertility rate (TFR) was then aggregated using the estimated age-specific fertility between ages 10 and 54 years. 7417 sources were used for under-5 mortality estimation and 7355 for adult mortality. ST-GPR was used to synthesise data sources after correction for known biases. Adult mortality was measured as the probability of death between ages 15 and 60 years based on vital registration, sample registration, and sibling histories, and was also estimated using ST-GPR. HIV-free life tables were then estimated using estimates of under-5 and adult mortality rates using a relational model life table system created for GBD, which closely tracks observed age-specific mortality rates from complete vital registration when available. Independent estimates of HIV-specific mortality generated by an epidemiological analysis of HIV prevalence surveys and antenatal clinic serosurveillance and other sources were incorporated into the estimates in countries with large epidemics. Annual and single-year age estimates of net migration and population for each country and territory were generated using a Bayesian hierarchical cohort component model that analysed estimated age-specific fertility and mortality rates along with 1250 censuses and 747 population registry years. We classified location-years into seven categories on the basis of the natural rate of increase in population (calculated by subtracting the crude death rate from the crude birth rate) and the net migration rate. We computed healthy life expectancy (HALE) using years lived with disability (YLDs) per capita, life tables, and standard demographic methods. Uncertainty was propagated throughout the demographic estimation process, including fertility, mortality, and population, with 1000 draw-level estimates produced for each metric. Findings: The global TFR decreased from 2·72 (95% uncertainty interval [UI] 2·66–2·79) in 2000 to 2·31 (2·17–2·46) in 2019. Global annual livebirths increased from 134·5 million (131·5–137·8) in 2000 to a peak of 139·6 million (133·0–146·9) in 2016. Global livebirths then declined to 135·3 million (127·2–144·1) in 2019. Of the 204 countries and territories included in this study, in 2019, 102 had a TFR lower than 2·1, which is considered a good approximation of replacement-level fertility. All countries in sub-Saharan Africa had TFRs above replacement level in 2019 and accounted for 27·1% (95% UI 26·4–27·8) of global livebirths. Global life expectancy at birth increased from 67·2 years (95% UI 66·8–67·6) in 2000 to 73·5 years (72·8–74·3) in 2019. The total number of deaths increased from 50·7 million (49·5–51·9) in 2000 to 56·5 million (53·7–59·2) in 2019. Under-5 deaths declined from 9·6 million (9·1–10·3) in 2000 to 5·0 million (4·3–6·0) in 2019. Global population increased by 25·7%, from 6·2 billion (6·0–6·3) in 2000 to 7·7 billion (7·5–8·0) in 2019. In 2019, 34 countries had negative natural rates of increase; in 17 of these, the population declined because immigration was not sufficient to counteract the negative rate of decline. Globally, HALE increased from 58·6 years (56·1–60·8) in 2000 to 63·5 years (60·8–66·1) in 2019. HALE increased in 202 of 204 countries and territories between 2000 and 2019

Global burden of chronic respiratory diseases and risk factors, 1990–2019: an update from the Global Burden of Disease Study 2019

Background: Updated data on chronic respiratory diseases (CRDs) are vital in their prevention, control, and treatment in the path to achieving the third UN Sustainable Development Goals (SDGs), a one-third reduction in premature mortality from non-communicable diseases by 2030. We provided global, regional, and national estimates of the burden of CRDs and their attributable risks from 1990 to 2019. Methods: Using data from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019, we estimated mortality, years lived with disability, years of life lost, disability-adjusted life years (DALYs), prevalence, and incidence of CRDs, i.e. chronic obstructive pulmonary disease (COPD), asthma, pneumoconiosis, interstitial lung disease and pulmonary sarcoidosis, and other CRDs, from 1990 to 2019 by sex, age, region, and Socio-demographic Index (SDI) in 204 countries and territories. Deaths and DALYs from CRDs attributable to each risk factor were estimated according to relative risks, risk exposure, and the theoretical minimum risk exposure level input. Findings: In 2019, CRDs were the third leading cause of death responsible for 4.0 million deaths (95% uncertainty interval 3.6–4.3) with a prevalence of 454.6 million cases (417.4–499.1) globally. While the total deaths and prevalence of CRDs have increased by 28.5% and 39.8%, the age-standardised rates have dropped by 41.7% and 16.9% from 1990 to 2019, respectively. COPD, with 212.3 million (200.4–225.1) prevalent cases, was the primary cause of deaths from CRDs, accounting for 3.3 million (2.9–3.6) deaths. With 262.4 million (224.1–309.5) prevalent cases, asthma had the highest prevalence among CRDs. The age-standardised rates of all burden measures of COPD, asthma, and pneumoconiosis have reduced globally from 1990 to 2019. Nevertheless, the age-standardised rates of incidence and prevalence of interstitial lung disease and pulmonary sarcoidosis have increased throughout this period. Low- and low-middle SDI countries had the highest age-standardised death and DALYs rates while the high SDI quintile had the highest prevalence rate of CRDs. The highest deaths and DALYs from CRDs were attributed to smoking globally, followed by air pollution and occupational risks. Non-optimal temperature and high body-mass index were additional risk factors for COPD and asthma, respectively. Interpretation: Albeit the age-standardised prevalence, death, and DALYs rates of CRDs have decreased, they still cause a substantial burden and deaths worldwide. The high death and DALYs rates in low and low-middle SDI countries highlights the urgent need for improved preventive, diagnostic, and therapeutic measures. Global strategies for tobacco control, enhancing air quality, reducing occupational hazards, and fostering clean cooking fuels are crucial steps in reducing the burden of CRDs, especially in low- and lower-middle income countries

Global burden of 369 diseases and injuries in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019

Background: In an era of shifting global agendas and expanded emphasis on non-communicable diseases and injuries along with communicable diseases, sound evidence on trends by cause at the national level is essential. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) provides a systematic scientific assessment of published, publicly available, and contributed data on incidence, prevalence, and mortality for a mutually exclusive and collectively exhaustive list of diseases and injuries. Methods: GBD estimates incidence, prevalence, mortality, years of life lost (YLLs), years lived with disability (YLDs), and disability-adjusted life-years (DALYs) due to 369 diseases and injuries, for two sexes, and for 204 countries and territories. Input data were extracted from censuses, household surveys, civil registration and vital statistics, disease registries, health service use, air pollution monitors, satellite imaging, disease notifications, and other sources. Cause-specific death rates and cause fractions were calculated using the Cause of Death Ensemble model and spatiotemporal Gaussian process regression. Cause-specific deaths were adjusted to match the total all-cause deaths calculated as part of the GBD population, fertility, and mortality estimates. Deaths were multiplied by standard life expectancy at each age to calculate YLLs. A Bayesian meta-regression modelling tool, DisMod-MR 2.1, was used to ensure consistency between incidence, prevalence, remission, excess mortality, and cause-specific mortality for most causes. Prevalence estimates were multiplied by disability weights for mutually exclusive sequelae of diseases and injuries to calculate YLDs. We considered results in the context of the Socio-demographic Index (SDI), a composite indicator of income per capita, years of schooling, and fertility rate in females younger than 25 years. Uncertainty intervals (UIs) were generated for every metric using the 25th and 975th ordered 1000 draw values of the posterior distribution. Findings: Global health has steadily improved over the past 30 years as measured by age-standardised DALY rates. After taking into account population growth and ageing, the absolute number of DALYs has remained stable. Since 2010, the pace of decline in global age-standardised DALY rates has accelerated in age groups younger than 50 years compared with the 1990–2010 time period, with the greatest annualised rate of decline occurring in the 0–9-year age group. Six infectious diseases were among the top ten causes of DALYs in children younger than 10 years in 2019: lower respiratory infections (ranked second), diarrhoeal diseases (third), malaria (fifth), meningitis (sixth), whooping cough (ninth), and sexually transmitted infections (which, in this age group, is fully accounted for by congenital syphilis; ranked tenth). In adolescents aged 10–24 years, three injury causes were among the top causes of DALYs: road injuries (ranked first), self-harm (third), and interpersonal violence (fifth). Five of the causes that were in the top ten for ages 10–24 years were also in the top ten in the 25–49-year age group: road injuries (ranked first), HIV/AIDS (second), low back pain (fourth), headache disorders (fifth), and depressive disorders (sixth). In 2019, ischaemic heart disease and stroke were the top-ranked causes of DALYs in both the 50–74-year and 75-years-and-older age groups. Since 1990, there has been a marked shift towards a greater proportion of burden due to YLDs from non-communicable diseases and injuries. In 2019, there were 11 countries where non-communicable disease and injury YLDs constituted more than half of all disease burden. Decreases in age-standardised DALY rates have accelerated over the past decade in countries at the lower end of the SDI range, while improvements have started to stagnate or even reverse in countries with higher SDI. Interpretation: As disability becomes an increasingly large component of disease burden and a larger component of health expenditure, greater research and developm nt investment is needed to identify new, more effective intervention strategies. With a rapidly ageing global population, the demands on health services to deal with disabling outcomes, which increase with age, will require policy makers to anticipate these changes. The mix of universal and more geographically specific influences on health reinforces the need for regular reporting on population health in detail and by underlying cause to help decision makers to identify success stories of disease control to emulate, as well as opportunities to improve. Funding: Bill & Melinda Gates Foundation. © 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 licens

Measuring universal health coverage based on an index of effective coverage of health services in 204 countries and territories, 1990–2019 : A systematic analysis for the Global Burden of Disease Study 2019

Background Achieving universal health coverage (UHC) involves all people receiving the health services they need, of high quality, without experiencing financial hardship. Making progress towards UHC is a policy priority for both countries and global institutions, as highlighted by the agenda of the UN Sustainable Development Goals (SDGs) and WHO's Thirteenth General Programme of Work (GPW13). Measuring effective coverage at the health-system level is important for understanding whether health services are aligned with countries' health profiles and are of sufficient quality to produce health gains for populations of all ages. Methods Based on the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019, we assessed UHC effective coverage for 204 countries and territories from 1990 to 2019. Drawing from a measurement framework developed through WHO's GPW13 consultation, we mapped 23 effective coverage indicators to a matrix representing health service types (eg, promotion, prevention, and treatment) and five population-age groups spanning from reproductive and newborn to older adults (≥65 years). Effective coverage indicators were based on intervention coverage or outcome-based measures such as mortality-to-incidence ratios to approximate access to quality care; outcome-based measures were transformed to values on a scale of 0–100 based on the 2·5th and 97·5th percentile of location-year values. We constructed the UHC effective coverage index by weighting each effective coverage indicator relative to its associated potential health gains, as measured by disability-adjusted life-years for each location-year and population-age group. For three tests of validity (content, known-groups, and convergent), UHC effective coverage index performance was generally better than that of other UHC service coverage indices from WHO (ie, the current metric for SDG indicator 3.8.1 on UHC service coverage), the World Bank, and GBD 2017. We quantified frontiers of UHC effective coverage performance on the basis of pooled health spending per capita, representing UHC effective coverage index levels achieved in 2019 relative to country-level government health spending, prepaid private expenditures, and development assistance for health. To assess current trajectories towards the GPW13 UHC billion target—1 billion more people benefiting from UHC by 2023—we estimated additional population equivalents with UHC effective coverage from 2018 to 2023. Findings Globally, performance on the UHC effective coverage index improved from 45·8 (95% uncertainty interval 44·2–47·5) in 1990 to 60·3 (58·7–61·9) in 2019, yet country-level UHC effective coverage in 2019 still spanned from 95 or higher in Japan and Iceland to lower than 25 in Somalia and the Central African Republic. Since 2010, sub-Saharan Africa showed accelerated gains on the UHC effective coverage index (at an average increase of 2·6% [1·9–3·3] per year up to 2019); by contrast, most other GBD super-regions had slowed rates of progress in 2010–2019 relative to 1990–2010. Many countries showed lagging performance on effective coverage indicators for non-communicable diseases relative to those for communicable diseases and maternal and child health, despite non-communicable diseases accounting for a greater proportion of potential health gains in 2019, suggesting that many health systems are not keeping pace with the rising non-communicable disease burden and associated population health needs. In 2019, the UHC effective coverage index was associated with pooled health spending per capita (r=0·79), although countries across the development spectrum had much lower UHC effective coverage than is potentially achievable relative to their health spending. Under maximum efficiency of translating health spending into UHC effective coverage performance, countries would need to reach $1398 pooled health spending per capita (US$ adjusted for purchasing power parity) in order to achieve 80 on the UHC effective coverage index. From 2018 to 2023, an estimated 388·9 million (358·6–421·3) more population equivalents would have UHC effective coverage, falling well short of the GPW13 target of 1 billion more people benefiting from UHC during this time. Current projections point to an estimated 3·1 billion (3·0–3·2) population equivalents still lacking UHC effective coverage in 2023, with nearly a third (968·1 million [903·5–1040·3]) residing in south Asia. Interpretation The present study demonstrates the utility of measuring effective coverage and its role in supporting improved health outcomes for all people—the ultimate goal of UHC and its achievement. Global ambitions to accelerate progress on UHC service coverage are increasingly unlikely unless concerted action on non-communicable diseases occurs and countries can better translate health spending into improved performance. Focusing on effective coverage and accounting for the world's evolving health needs lays the groundwork for better understanding how close—or how far—all populations are in benefiting from UHC

Mapping local patterns of childhood overweight and wasting in low- and middle-income countries between 2000 and 2017

A double burden of malnutrition occurs when individuals, household members or communities experience both undernutrition and overweight. Here, we show geospatial estimates of overweight and wasting prevalence among children under 5 years of age in 105 low- and middle-income countries (LMICs) from 2000 to 2017 and aggregate these to policy-relevant administrative units. Wasting decreased overall across LMICs between 2000 and 2017, from 8.4% (62.3 (55.1–70.8) million) to 6.4% (58.3 (47.6–70.7) million), but is predicted to remain above the World Health Organization’s Global Nutrition Target of <5% in over half of LMICs by 2025. Prevalence of overweight increased from 5.2% (30 (22.8–38.5) million) in 2000 to 6.0% (55.5 (44.8–67.9) million) children aged under 5 years in 2017. Areas most affected by double burden of malnutrition were located in Indonesia, Thailand, southeastern China, Botswana, Cameroon and central Nigeria. Our estimates provide a new perspective to researchers, policy makers and public health agencies in their efforts to address this global childhood syndemic

Measuring universal health coverage based on an index of effective coverage of health services in 204 countries and territories, 1990–2019: a systematic analysis for the Global Burden of Disease Study 2019

Background Achieving universal health coverage (UHC) involves all people receiving the health services they need, of high quality, without experiencing financial hardship. Making progress towards UHC is a policy priority for both countries and global institutions, as highlighted by the agenda of the UN Sustainable Development Goals (SDGs) and WHO's Thirteenth General Programme of Work (GPW13). Measuring effective coverage at the health-system level is important for understanding whether health services are aligned with countries' health profiles and are of sufficient quality to produce health gains for populations of all ages. Methods Based on the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2019, we assessed UHC effective coverage for 204 countries and territories from 1990 to 2019. Drawing from a measurement framework developed through WHO's GPW13 consultation, we mapped 23 effective coverage indicators to a matrix representing health service types (eg, promotion, prevention, and treatment) and five population-age groups spanning from reproductive and newborn to older adults (≥65 years). Effective coverage indicators were based on intervention coverage or outcome-based measures such as mortality-to-incidence ratios to approximate access to quality care; outcome-based measures were transformed to values on a scale of 0–100 based on the 2·5th and 97·5th percentile of location-year values. We constructed the UHC effective coverage index by weighting each effective coverage indicator relative to its associated potential health gains, as measured by disability-adjusted life-years for each location-year and population-age group. For three tests of validity (content, known-groups, and convergent), UHC effective coverage index performance was generally better than that of other UHC service coverage indices from WHO (ie, the current metric for SDG indicator 3.8.1 on UHC service coverage), the World Bank, and GBD 2017. We quantified frontiers of UHC effective coverage performance on the basis of pooled health spending per capita, representing UHC effective coverage index levels achieved in 2019 relative to country-level government health spending, prepaid private expenditures, and development assistance for health. To assess current trajectories towards the GPW13 UHC billion target—1 billion more people benefiting from UHC by 2023—we estimated additional population equivalents with UHC effective coverage from 2018 to 2023. Findings Globally, performance on the UHC effective coverage index improved from 45·8 (95% uncertainty interval 44·2–47·5) in 1990 to 60·3 (58·7–61·9) in 2019, yet country-level UHC effective coverage in 2019 still spanned from 95 or higher in Japan and Iceland to lower than 25 in Somalia and the Central African Republic. Since 2010, sub-Saharan Africa showed accelerated gains on the UHC effective coverage index (at an average increase of 2·6% [1·9–3·3] per year up to 2019); by contrast, most other GBD super-regions had slowed rates of progress in 2010–2019 relative to 1990–2010. Many countries showed lagging performance on effective coverage indicators for non-communicable diseases relative to those for communicable diseases and maternal and child health, despite non-communicable diseases accounting for a greater proportion of potential health gains in 2019, suggesting that many health systems are not keeping pace with the rising non-communicable disease burden and associated population health needs. In 2019, the UHC effective coverage index was associated with pooled health spending per capita (r=0·79), although countries across the development spectrum had much lower UHC effective coverage than is potentially achievable relative to their health spending. Under maximum efficiency of translating health spending into UHC effective coverage performance, countries would need to reach $1398 pooled health spending per capita (US$ adjusted for purchasing power parity) in order to achieve 80 on the UHC effective coverage index. From 2018 to 2023, an estimated 388·9 million (358·6–421·3) more population equivalents would have UHC effective coverage, falling well short of the GPW13 target of 1 billion more people benefiting from UHC during this time. Current projections point to an estimated 3·1 billion (3·0–3·2) population equivalents still lacking UHC effective coverage in 2023, with nearly a third (968·1 million [903·5–1040·3]) residing in south Asia. Interpretation The present study demonstrates the utility of measuring effective coverage and its role in supporting improved health outcomes for all people—the ultimate goal of UHC and its achievement. Global ambitions to accelerate progress on UHC service coverage are increasingly unlikely unless concerted action on non-communicable diseases occurs and countries can better translate health spending into improved performance. Focusing on effective coverage and accounting for the world's evolving health needs lays the groundwork for better understanding how close—or how far—all populations are in benefiting from UHC

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 0.71 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 50.2% exceed this threshold for suitability in at least one 5×5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify

Predicting the environmental suitability for onchocerciasis in Africa as an aid to elimination planning

Recent evidence suggests that, in some foci, elimination of onchocerciasis from Africa may be feasible with mass drug administration (MDA) of ivermectin. To achieve continental elimination of transmission, mapping surveys will need to be conducted across all implementation units (IUs) for which endemicity status is currently unknown. Using boosted regression tree models with optimised hyperparameter selection, we estimated environmental suitability for onchocerciasis at the 5 × 5-km resolution across Africa. In order to classify IUs that include locations that are environmentally suitable, we used receiver operating characteristic (ROC) analysis to identify an optimal threshold for suitability concordant with locations where onchocerciasis has been previously detected. This threshold value was then used to classify IUs (more suitable or less suitable) based on the location within the IU with the largest mean prediction. Mean estimates of environmental suitability suggest large areas across West and Central Africa, as well as focal areas of East Africa, are suitable for onchocerciasis transmission, consistent with the presence of current control and elimination of transmission efforts. The ROC analysis identified a mean environmental suitability index of 071 as a threshold to classify based on the location with the largest mean prediction within the IU. Of the IUs considered for mapping surveys, 502% exceed this threshold for suitability in at least one 5 × 5-km location. The formidable scale of data collection required to map onchocerciasis endemicity across the African continent presents an opportunity to use spatial data to identify areas likely to be suitable for onchocerciasis transmission. National onchocerciasis elimination programmes may wish to consider prioritising these IUs for mapping surveys as human resources, laboratory capacity, and programmatic schedules may constrain survey implementation, and possibly delaying MDA initiation in areas that would ultimately qualify.SUPPORTING INFORMATION : FIGURE S1. Data coverage by year. Here we visualise the volume of data used in the analysis by country and year. Larger circles indicate more data inputs. ‘NA’ indicates records for which no year was reported (eg, ‘pre-2000’). https://doi.org/10.1371/journal.pntd.0008824.s001FIGURE S2. Illustration of covariate values for year 2000. Maps were produced using ArcGIS Desktop 10.6. https://doi.org/10.1371/journal.pntd.0008824.s002FIGURE S3. Environmental suitability of onchocerciasis including locations that have received MDA for which no pre-intervention data are available. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s003FIGURE S4. Environmental suitability prediction uncertainty including locations that have received MDA for which no pre-intervention data are available. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s004FIGURE S5. Environmental suitability of onchocerciasis excluding morbidity data. This plot shows suitability predictions from green (low = 0%) to pink (high = 100%), representing those areas where environmental conditions are most similar to prior pathogen detections. Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s005FIGURE S6. Environmental suitability prediction uncertainty excluding morbidity data. This plot shows uncertainty associated with environmental suitability predictions colored from blue to red (least to most uncertain). Countries in grey with hatch marks were excluded from the analysis based on a review of national endemicity status. Areas in grey only represent locations masked due to sparse population. https://doi.org/10.1371/journal.pntd.0008824.s006FIGURE S7. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. The figures show the fit per covariate relative to the data that correspond to specific values of the covariate. https://doi.org/10.1371/journal.pntd.0008824.s007FIGURE S8. Covariate Effect Curves for all onchocerciasis occurrences (measures of infection prevalence and disability). On the right set of axes we show the frequency density of the occurrences taking covariate values over 20 bins of the horizontal axis. The left set of axes shows the effect of each on the model, where the mean effect is plotted on the black line and its uncertainty is represented by the upper and lower confidence interval bounds plotted in dark grey. https://doi.org/10.1371/journal.pntd.0008824.s008FIGURE S9. ROC analysis for threshold. Results of the area under the receiver operating characteristic (ROC) curve analysis are presented below, with false positive rate (FPR) on the x-axis and true positive rate (TPR) on the y-axis. The red dot on the curve represents the location on the curve that corresponds to a threshold that most closely agreed with the input data. For each of the 100 BRT models, we estimated the optimal threshold that maximised agreement between occurrence inputs (considered true positives) and the mean model predictions as 0·71. https://doi.org/10.1371/journal.pntd.0008824.s009TABLE S1. Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) checklist. https://doi.org/10.1371/journal.pntd.0008824.s010TABLE S2. Total number of occurrence data classified as point and polygon inputs by diagnostic. We present the total number of occurrence points extracted from the input data sources by diagnostic type. ‘Other diagnostics’ include: DEC Patch test; Knott’s Method (Mazotti Test); 2 types of LAMP; blood smears; and urine tests. https://doi.org/10.1371/journal.pntd.0008824.s011TABLE S3. Total number of occurrence data classified as point and polygon inputs by location. https://doi.org/10.1371/journal.pntd.0008824.s012TABLE S4. Covariate information. https://doi.org/10.1371/journal.pntd.0008824.s013TEXT S1. Details outlining construction of occurrence dataset. https://doi.org/10.1371/journal.pntd.0008824.s014TEXT S2. Covariate rationale. https://doi.org/10.1371/journal.pntd.0008824.s015TEXT S3. Boosted regression tree methodology additional details. https://doi.org/10.1371/journal.pntd.0008824.s016APPENDIX S1. Country-level maps and data results. Maps were produced using ArcGIS Desktop 10.6 and shapefiles to visualize administrative units are available at https://espen.afro.who.int/tools-resources/cartography-database. https://doi.org/10.1371/journal.pntd.0008824.s017This work was primarily supported by a grant from the Bill & Melinda Gates Foundation OPP1132415 (SIH). Financial support from the Neglected Tropical Disease Modelling Consortium (https://www.ntdmodelling.org/), which is funded by the Bill & Melinda Gates Foundation (grants No. OPP1184344 and OPP1186851), and joint centre funding (grant No. MR/R015600/1) by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).The Neglected Tropical Disease Modelling Consortium which is funded by the Bill & Melinda Gates Foundation, the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement which is also part of the EDCTP2 programme supported by the European Union (MGB).http://www.plosNTDS.orgam2022Medical Microbiolog