A difference-in-difference, quasi-experimental design was used to estimate a “treatment effect” (PEPFAR), based on comparison to a control group (the counterfactual). The difference-in-difference design compares the before and after change in outcomes for the treatment group to the before and after change in outcomes for the comparison group. The outcomes of interest, their definitions and sources are listed in Table 1. Baseline variables are listed in Table 2. The panel data set of 157 low- and middle- income countries used in the prior analysis, covering 1990 to 2018, was updated to include data through 2022 (due to data limitations, the time trend for maternal mortality began in 2000). COVID-related funding was not included. All values were adjusted to constant 2022 dollars.
The PEPFAR group included 90 countries that had received PEPFAR support over the period. The comparison group included 67 low- and middle- income countries that had not received any PEPFAR support or had received minimal PEPFAR support (<$1M over the period or <$.05 per capita) between 2004 and 2022. Data on PEPFAR spending by country were obtained from the U.S. government’s https://foreignassistance.gov/ database and represent U.S. fiscal year disbursements. Data on mortality rates were obtained from the World Bank’s World Development Indicators (WDI) (https://datatopics.worldbank.org/world-development-indicators/. Several different model specifications were explored. Each specification controlled for numerous baseline variables, compared to an unadjusted model, variables which may be expected to influence the outcomes of interest and which help make the comparison group more comparable to the PEPFAR group.
Table 3 provides the model specifications tested in the updated analysis. Each model specification produced similar, statistically significant results. All models were also run with and without China and India, the two most populous countries in the world, to assess whether they were influencing the results. In both cases, PEPFAR’s impact was still significant and results were similar. The final reported results are from model specification 3. The pre-intervention period for this model started in 2002. Almost all results were significant at the p<0.001 level; one result was significant at the p<0.01 level and one at the p<0.05 level. Table 4 provides the mean values for baseline outcomes and Tables 5-6 provide model results. The Table 5 difference-in-difference estimates should be interpreted as the unit change (e.g., in number of deaths per 1,000) in the outcome associated with PEPFAR. The Table 6 estimates should be interpreted as the percent change in the outcome, relative to the baseline, associated with PEPFAR.
Despite the strengths of the difference-in-difference design, there are limitations to this approach. While the models adjusted for numerous baseline factors that could be correlated with the outcomes of interest, there may be other, unobservable factors not captured. Similarly, while baseline factors are also intended to adjust for selection bias, and make the PEPFAR and comparison groups more similar, there may be other ways in which comparison countries differed from PEPFAR countries (and factors which influenced which countries received PEPFAR support), which could bias the estimates. Recent, published research, based on the earlier 2004-2018 period, examined several other model specifications as well as double log transformations for each model. The results were similar across all models, adding to the confidence of the analytic approach used here.4 Another analysis, under review for publication, conducted additional robustness tests and found consistent results.5
| Table 1: Outcome Variables |
| Variable |
Data Source |
| 1. Crude death rate |
All causes per 1,000 population |
| 2. Child mortality rate |
Probability of a child dying between birth and 5 years of age per 1,000 live births |
| 3. Maternal mortality ratio |
Number of women who die from pregnancy-related causes while pregnant or within 42 days of pregnancy termination per 100,000 live births |
| 4. Measles immunization |
Percent of children ages 12-23 months who received the measles vaccination |
| 5. DPT immunizations |
Percent of children ages 12-23 months who received DPT vaccinations (3 doses) |
| 6. Hepatitis B immunizations |
Percent of children ages 12-23 months who received hepatitis B vaccinations (3 doses) |
| 7. Newborns protected against tetanus |
Percentage of births by women of child-bearing age who are immunized against tetanus |
| Source: World Bank, WDI, https://datatopics.worldbank.org/world-development-indicators/ |
| Table 2: Baseline Variables |
| Variable |
Data Source |
| 1. Gross Domestic Product (GDP) per capita (current USD) |
World Bank Development Indicators |
| 2. Recipient of U.S. HIV funding prior to 2004 (dummy variable) |
https://foreignassistance.gov/ |
| 3. Total population |
United Nations, Department of Economic and Social Affairs, Population Division |
| 4. Life expectancy at birth (years) |
World Bank Development Indicators |
| 5. Total fertility rate (births per woman) |
World Bank Development Indicators |
| 6. Percent urban population (of total population) |
World Bank Development Indicators |
| 7. School enrollment, secondary (% gross) |
World Bank Development Indicators |
| 8. World Bank country income classification |
World Bank |
| 9. HIV prevalence (% of population ages 15-49) |
World Bank Development Indicators. To address missing values in some cases, additional data were obtained from the Global Burden of Disease Collaborative Network |
| 10. Per capita donor spending on health (non-PEPFAR) (constant $) |
OECD Creditor Reporting System database |
| 11. Per capita domestic health spending, government and private, PPP (current $) |
World Bank Development Indicators |
| 12. Measles prevalence in under 5 population (measles immunization models only) |
IHME, http://ghdx.healthdata.org/gbd-results-tool |
| 13. Diphtheria prevalence in under 5 population (DPT immunization models only) |
IHME, http://ghdx.healthdata.org/gbd-results-tool |
| 14. Whooping cough prevalence 15. in under 5 population (DPT immunization models only) |
IHME, http://ghdx.healthdata.org/gbd-results-tool |
| 15. Tetanus prevalence in under 5 population (DPT immunization models only) |
IHME, http://ghdx.healthdata.org/gbd-results-tool |
| 16. Hepatitis B prevalence in under 5 population (Hepatitis B immunization models only) |
IHME, http://ghdx.healthdata.org/gbd-results-tool |
| Table 3: Model Specifications |
| 1. Unadjusted model |
| 2. Includes baseline variables 1-9 (and an additional baseline variable for disease incidence, 12-16, depending on outcome measure) |
| 3. Includes baseline variables 1-11 (and an additional baseline variable for disease incidence, 12-16, depending on outcome measure) |
| Table 4: Baseline Mean Outcome Values, 2004 |
| Outcome |
All PEPFAR Countries |
PEPFAR COP Countries |
| All-cause mortality rate |
10.5 |
12.6 |
| Maternal mortality ratio |
409.8 |
497.5 |
| Child mortality rate |
78.9 |
97.9 |
| Measles immunization percent |
77.5 |
74.5 |
| DPT immunization percent |
78.0 |
74.2 |
| Hepatitis B immunization percent |
79.1 |
74.4 |
| Prevalence of newborns protected against tetanus |
75.7 |
75.6 |
| Table 5: Difference-in-Difference Estimates Associated with PEPFAR, 2004-2022 (standard errors in parentheses) |
| Outcome |
All PEPFAR Countries |
PEPFAR COP Countries |
| All-cause mortality rate |
-2.452*** |
-3.523*** |
| |
(0.171) |
(0.235) |
| Maternal mortality ratio |
-119.1*** |
-155.9*** |
| |
(14.618) |
(19.101) |
| Child mortality rate |
-30.02*** |
-39.48*** |
| |
(1.329) |
(1.979) |
| Measles immunization percent |
7.004*** |
7.795*** |
| |
(0.678) |
(0.956) |
| DPT immunization percent |
8.871*** |
9.022*** |
| |
(0.626) |
(0.880) |
| Hepatitis B immunization percent |
6.435** |
13.67*** |
| |
(2.133) |
(2.915) |
| Prevalence of newborns protected against tetanus |
4.962*** |
3.303* |
| |
(1.299) |
(1.495) |
| ***p < 0.001 **p < 0.01 *p < 0.05 |
| Table 6: Estimated Percent Change Associated with PEPFAR, 2004-2022 (Relative to 2004 Baseline) |
| Outcome |
All PEPFAR Countries |
PEPFAR COP Countries |
| All-cause mortality rate |
-23.4%*** |
-28.0%*** |
| Maternal mortality ratio |
-29.1%*** |
-31.3%*** |
| Child mortality rate |
-38.0%*** |
-40.3%*** |
| Measles immunization percent |
9.0%*** |
10.5%*** |
| DPT immunization percent |
11.4%*** |
12.2%*** |
| Hepatitis B immunization percent |
8.1%** |
18.4%*** |
| Prevalence of newborns protected against tetanus |
6.6%*** |
4.4%* |
| ***p < 0.001 **p < 0.01 *p < 0.05 |
| Table 7: Estimated Incremental Percent Change Associated with PEPFAR by Time Period (Relative to 2004 Baseline) |
| Outcome |
2004 – 2008 |
2009 – 2013 |
2014 – 2018 |
2019 – 2022 |
| All-cause mortality rate |
-9.8% |
-5.7% |
-5.0% |
-2.8% |
| Maternal mortality ratio |
-14.2% |
-6.3% |
-4.1% |
-4.5% |
| Child mortality rate |
-22.0% |
-7.0% |
-5.7% |
-3.3% |
| Measles immunization percent |
5.3% |
2.4% |
1.2% |
0.1% |
| DPT immunization percent |
7.9% |
2.3% |
1.1% |
0.1% |
| Hepatitis B immunization percent |
6.8% |
0.6% |
1.0% |
-0.3% |
| Prevalence of newborns protected against tetanus |
9.0% |
-1.3% |
-0.8% |
-0.4% |