Impact of Health System Inputs on Health Outcome: A Multilevel Longitudinal Analysis of Botswana National Antiretroviral Program (2002-2013)

Objective: To measure the association between the number of doctors, nurses and hospital beds per 10,000 people and individual HIV-infected patient outcomes in Botswana. Design: Analysis of routinely collected longitudinal data from 97,627 patients who received ART through the Botswana National HIV/AIDS Treatment Program across all 24 health districts from 2002 to 2013. Doctors, nurses, and hospital bed density data at district-level were collected from various sources. Methods: A multilevel, longitudinal analysis method was used to analyze the data at both patient- and district-level simultaneously to measure the impact of the health system input at district-level on probability of death or loss-to-follow-up (LTFU) at the individual level. A marginal structural model was used to account for LTFU over time. Results: Increasing doctor density from one doctor to two doctors per 10,000 population decreased the predicted probability of death for each patient by 27%. Nurse density changes from 20 nurses to 25 nurses decreased the predicted probability of death by 28%. Nine percent decrease was noted in predicted mortality of an individual in the Masa program for every five hospital bed density increase. Conclusion: Considerable variation was observed in doctors, nurses, and hospital bed density across health districts. Predictive margins of mortality and LTFU were inversely correlated with doctor, nurse and hospital bed density. The doctor density had much greater impact than nurse or bed density on mortality or LTFU of individual patients. While long-term investment in training more healthcare professionals should be made, redistribution of available doctors and nurses can be a feasible solution in the short term.