Effects of AIDS-related disability on workforce participation and earned income in Botswana: A quasi-experimental evaluation ()
1. INTRODUCTION
In Botswana, over 350 thousand people are estimated to be living with HIV/AIDS [1]. The vast majority of HIV cases in the country are among working aged adults between 15 and 49, with a prevalence in this age group of 25% [2]. The potential for reduced labor productivity of individuals who make up the majority of the workforce in the country has implications for the economic wellbeing of their households as well as the growth of the economy in Botswana. Several studies have demonstrated the detrimental macroeconomic effects of AIDS in countries with a high prevalence [3-13]. A recent model estimates that the growth rate of the economy in Botswana will slow by up to 2% annually through 2021, rendering it one-fifth smaller than it might have been in the absence of AIDS, even considering widespread antiretroviral therapy [8]. An estimated 40 thousand uninfected people in Botswana have been pulled into poverty due to the slow economic growth caused by the epidemic [13], reflecting similar results from other settings showing the effect of the epidemic on household economies [14-22]. The effects of the epidemic related to absenteeism, workforce participation, and earnings could therefore have both shortand long-term implications for the economy in Botswana.
In recognition of such far-reaching social and economic consequences of the epidemic, the Government of Botswana established a national antiretroviral therapy (ART) program in 2002, known as Masa or “new dawn” in Setswana, and began rapid scale-up across the country in 2004. The Masa program provides antiretroviral drugs, essential laboratory testing, and other care and treatment services free of charge to any individual who has tested positive for HIV and qualifies for treatment based on virologic or clinical criteria. The Masa program has a dual purpose to improve the length and quality of life for individuals with HIV and defray the enormous cost of life-long treatment, easing financial burdens on households and individuals.
Several years after the Masa program has been established and expanded nation-wide, research has shown dramatic reductions in HIV mortality rates and promising clinical outcomes among individuals enrolled in the program [23]. Important questions still remain to be answered regarding the financial well-being of individuals with HIV in the country. Therefore, two important factors motivate the evaluation of the indirect economic effects of HIV-related illness in Botswana. First, as HIV becomes more of a chronic condition and the life expectancy of those with HIV is greatly extended thanks to the successful scale up of treatment and care services, workforce participation and economic productivity may be affected to a greater degree by illness and disability rather than death [24]. Secondly, given that the national ART program is fully subsidized in Botswana, we assume that direct costs for care and treatment are of lesser concern in this context than indirect costs, such as earnings potentially foregone via absenteeism.
There is abundant evidence from developed countries linking HIV-related morbidity to increased absenteeism and reduced productivity and wages along a gradient of disease severity [25-37]. Only a small number of studies examine these effects in developing countries [38-42]. More often, studies conducted in the developing world demonstrate adverse economic effects of the epidemic at the household level, focusing mostly on expenditures for care and treatment rather than potential opportunity costs, and more frequently examine the effects of deaths due to HIV/AIDS rather than illness [7,8,12,14-22,39,43,44].
Our analysis therefore adds to the evidence base in terms of indirect economic impacts of HIV at the individual level in a developing country context. We use a random sample of adults from across the country, and rigorous quasi-experimental methods to evaluate the effects of HIV-related illness on absenteeism and earnings in Botswana, with special attention to differences between men and women.
2. METHODS
2.1. Data
We surveyed 1840 households containing 6480 individuals randomly selected in three urban and three rural sites across Botswana over three months in 2007. Overall HIV prevalence in these six locations ranges from 13.3% to 26.5% [2]. One member of each household, most often the head of household, was asked to respond to the questionnaire and provide information for each household member. Informed consent was obtained from all respondents before the administration of the questionnaire, and ethical review committee approvals were acquired from the Harvard School of Public Health and the Botswana Ministry of Health before the survey.
We included 3999 individuals from 1768 households between the ages of 15 to 49 with complete information in our analysis, representing 62% of the individuals from our randomly selected population. Since the survey instrument was randomly administered to individuals in their home, sometimes in the presence of others in the household, it was not consistent with ethical human research practices to ask respondents directly about their HIV infection status. Instead, they were asked to list all episodes of illness within the last year and last month, and attribute a cause to each illness episode. The “HIVrelated illness” variable is therefore an indicator equal to 1 if an individual attributed any episode of illness in the past year or month to HIV through self-report, and 0 otherwise.
Outcomes of interest include absenteeism and earnings in the last month. Absenteeism was defined as the number of days last month that the individual was unable to perform normal activities, including any inpatient days. Absenteeism questions were asked uniformly for all episodes of illness for each member of every household surveyed. The comparison we make is therefore the number of days of absenteeism reported by those who attribute the illness episode to HIV, versus those who attribute it to any other cause. Earnings include all reported income from paid work, and these skewed data were transformed with the natural logarithm for the analysis. Our indicator of workforce participation is based on the earnings data, and is defined such that anyone who reported a positive amount of earnings in the past month was considered part of the workforce, while those reporting no earnings from paid work in the last month were considered not part of the workforce.
2.2. Propensity Score Matching
We apply a propensity score matching approach to compare self-reports of absenteeism in the last month between individuals reporting HIV-related illness and a group of matched controls [45-47]. Propensity score matching methods allow us to simulate “treatment” and “control” groups by matching on a few observed covariates, resulting in a quasi-randomized design that mimics the result of balance achieved by randomization [46,48,49]. We use three matching methods in order to validate the results, including nearest-neighbor, kernel, and radius matching [45,50]. Once individuals are matched, we can compare days of absenteeism between those who have reported HIV-related illness and those who have not. This is an estimate of the “average treatment effect on the treated” (ATT) of HIV-related illness. We stratified the analysis on sex, and patients in either group who are not matched are eliminated from the comparison. It is important to note that while propensity score methods can reduce bias due to observed differences between treatment groups, they are still subject to biases from unobserved differences.
2.3. Heckman Selection Model
We evaluate whether HIV-related illness has an impact on earnings using a Heckman selection model, which allows us to control for the fact that not all individuals enter the workforce. This conditional analysis is important because we cannot extrapolate the effect of HIVrelated illness on earnings that we observe among workers to non-workers since many of the unobservable characteristics that lead people to forego employment, e.g. worse health status, also lead them to less success once employed [51]. If people do self-select into the sample of workers, we may obtain biased estimates without properly adjusting for that phenomenon, as the combination of factors leading to entrance into the workforce would be similar to an omitted variable in the regression [52, 53].
The Heckman selection model addresses this potential bias by first estimating the probability of entering the workforce, as a function of a series of observable and measured covariates, using a probit regression. The degree to which there is selection in the sample based on workforce participation is estimated in a term that is entered into the main linear regression to evaluate earnings among those who participated in the workforce last month. In this model, if the correlation coefficient rho (ρ) given in the model output is not equal to zero; we have evidence of sample selection, meaning that indeed the observed characteristics in the probit regression have influenced participation in the workforce [54]. In all of our analyses, we specify our use of survey data with enumeration areas as the primary sampling unit, and household weights as the sampling weights. We used STATA 11 software to carry out our analyses [55].
3. RESULTS
3.1. Descriptive Statistics
Sample characteristics and summary statistics are presented in Table 1, stratified by sex. The samples of men and women are similar in many ways. The average age is
Table 1. Sample characteristics and descriptive statistics.

*Urban includes Gaborone, Selebi-Phikwe, and Jwaneng; Rural includes Kanye, Ghanzi, and Hukuntsi; **Earnings rounded to the nearest Pula (7.7 Pula to the US dollar), only presented for those reporting a positive amount of earnings in the last month.
about 28 years old, and two-thirds of both male and female respondents have achieved secondary education. More men have achieved some form of higher education (20%) compared to women (15%). Women live in households with higher dependency ratios (0.50) compared to men (0.33). In our sample, 72% of individuals lived in rural areas, and 28% in urban areas. We also constructed a household wealth index using ownership of household assets in a principal component analysis [56], and partitioned this into five quintiles. Overall l9% of males and 22% of females fell into the lowest quintile, while 21% of males and 19% of females fell into the highest quintile. Approximately 2% of males and females alike reported HIV-related illness. The average age among individuals reporting HIV-related illness in the last year was 35%, 57% were female, 31% were married or living with a partner, 49% were residents of rural areas, 64% were educated at the secondary or tertiary level, and 66% fell into the lowest wealth quintile. In total across the whole sample, men and women reported an average of one half-day of absenteeism in the last month (0.51 and 0.54 respectively). However, excluding those who reported no absenteeism, men report on average 9.79 days and women report 8.24 days of absence in the last month. More men (56%) participated in the workforce based on positive earnings data compared to women (45%). Among these individuals, men made an average of about 2871 Pula (US $372) in the last month, and women made an average of 2194 Pula (US $284) in the last month.
3.2. Outcomes
Absenteeism: The results of the logistic regression (not shown) used to calculate the propensity scores indicate that the probability of reporting HIV-related illness in our sample rises with age at first and then declines. The peak age at which HIV-related illness is reported is 44 for men and 48 for women. Rural residents and people from a low socioeconomic status are more likely to report HIV-related illness than urban residents and individuals of middle and high socioeconomic status, respectively. Table 2 presents our findings from the propensity score matching comparisons. We report the average treatment effect on the treated (ATT), the average difference in absenteeism reported between persons reporting HIVrelated illness and their matched set of controls. Men reporting HIV-related illness are more likely to report longer periods of absenteeism in the previous month by
Table 2. Absenteeism due to HIV-related illness, propensity score matching.

Estimates represent the average treatment effect among the treated (ATT); Standard errors in parentheses, Kernel matching uses bootstrapped standard errors, all other methods use analytical standard errors, *p < 0.10, **p < 0.05, *** p < 0.01.
approximately 5.2 days, while women report about 3.3 additional days, compared to their matched controls.
Earnings: The results of the Heckman Selection Model are presented in Table 3. The first step in comparing earnings was controlling for self-selection into the workforce. The results of this step align with general expectations. HIV-related illness appears to have a negative effect on participation in the workforce for both sexes, and the effect is stronger and highly significant among females. Age is significant and positive for both sexes and the effect drops off at older ages as the age-squared term is negative and significant. Both men and women in rural areas appear to be less likely to par-
Table 3. Impact of HIV-related illness on Earnings, Heckman Selection Model (MLE).

Standard errors in parentheses; *p < 0.10, **p < 0.05, ***p < 0.01; The natural logarithm is used to transform the skewed earnings data. Earnings data includes all reports of positive amounts.
ticipate in the workforce, though the effect is significant only among men. The positive sign of the education variables suggests that higher education is correlated with an increased probability of working. The effects of secondary education and tertiary education are stronger and significant among females. A larger dependency ratio in the household appears to discourage working, which is significant for both men and women but a much stronger effect among men. Among males, being married or living with a partner has a strong positive effect on workforce participation, but a negative and significant effect on women’s workforce participation. Men in the highest wealth quintile are significantly more likely to work than those in the lowest, while the effect is reversed for women, who are much less likely to participate in the workforce in the highest quintile compared to the lowest.
In the main linear regression, we observe that the effect of HIV-related illness on earnings is negative and significant for both sexes. Since we used the natural logarithm to transform the skewed earnings data, we exponentiate the coefficients to quantify the actual effect [57]. The effect of HIV-related illness is an approximately 38% reduction in last month’s earnings among men, and 43% reduction among women. Among males, it seems that earnings increase with age but the effect is not statistically significant. For women, age has a positive and significant effect on earnings until the age of 53. Primary education is not statistically different from “no schooling” with regard to earnings, while earnings increase with secondary and tertiary levels of education for both sexes. The education effects are stronger among women. Interaction terms between HIV-related illness, age, and education level were insignificant in the model.
4. DISCUSSION
We find that working age adults reporting HIV-related illness in Botswana are reporting greater amounts of absenteeism, and that HIV-related illness contributes to reduced earnings. In addition, we find that the difference in absenteeism between those with HIV-related illness and without is smaller among women than among men, but that HIV-related illness is associated with reduced earnings of almost half of the previous month’s earnings for women, and reduced earnings of about one-third for men. Our findings are consistent with studies from both the developed and developing world showing that HIV morbidity is associated with reduced likelihood of workforce participation and earnings [33-38,40,42]. However, our findings differ slightly in magnitude compared with previous research.
Specifically, the findings from this study show rates of absenteeism similar or higher to the upper bound from previous studies, and reductions in earnings at least double that from other studies in sub-Saharan Africa. The rates in absenteeism observed in the present analysis are very similar to reports from other studies specific to individuals in the month before they began initiation of antiretroviral therapy (ART) [38,40,42]. Mahal et al. (2008) use propensity score matching to show that individuals with HIV report about 2 additional days of absenteeism in the last month, compared with those without HIV [41]. Data from employees of the Debswana Diamond mining firm in Botswana have been used to show that absenteeism among workers at the firm with HIV peaks at approximately 5 days in the month prior to antiretroviral treatment (ART) initiation, and these researchers found no difference in absenteeism between men and women [39].
Our study offers several advantages over previous research. The sample we use is a large, random sample of adults from across the country, rather than a convenience sample. We use rigorous quasi-experimental methods, and focus specifically on differences between men and women. The study also has some limitations. Interestingly, the results regarding absenteeism in this study are very similar to estimates from previous studies specific to individuals who have not yet initiated ART. Our measurement of HIV-related illness and absenteeism rely on self-reports. It is possible that some individuals may have been hesitant to reveal their HIV status, and attributed their illness episode to some other cause. However, if this occurred, the result would be an underestimation of the impact of HIV-related illness on absenteeism on earnings.
In addition, we expect individuals to be sick at two distinct periods throughout their disease progression: within one year before and after their initiation on antiretroviral treatment (ART), and later when the disease has progressed substantially. If we were to make the assumption, considering the widespread voluntary testing and counseling campaign in Botswana in past years, that a self-report of HIV-related illness means that the individual has been tested and knows their status, the rate of absenteeism we see in this analysis is surprisingly high, since we might have expected these individuals to have already enrolled in the Masa program and been initiated on ART if they were eligible for treatment. It is possible that individuals could have been tested and not initiated treatment because they were waiting to qualify by virologic or clinical criteria, but this degree of absenteeism suggests that many might be eligible if not already on treatment. It is also possible that these individuals are simply not adhering to their treatment regimen, or that they have been on treatment for longer periods of time and are reaching later phases of disease progression. A limitation of this study is the inability to capture this element of treatment status, duration, and adherence among individuals in the sample.
Stigma may play a role in the results we observe. Discrimination from employers could account for part of the difference in workforce participation and earnings observed among those reporting HIV-related illness. Without an indicator for stigma or discrimination that may have influenced employment or earnings, we could not distinguish whether individuals who did not report any earnings in the last month had been seeking work or were denied employment for any reason related to their illness. Enterprises may refrain from hiring HIV positive workers for various reasons including loss of investment in training, increased costs related to health, or reduced productivity, despite the increase in lifespan afforded by effective treatment. However, unpublished work using data from Botswana National ART Program database suggests that individuals often present late into the course of their disease for testing and treatment [58]. This suggests some delay with respect to testing or enrolling in treatment, which may explain some of the high reports of absenteeism in our data.
5. RECOMMENDATIONS
Our findings suggest at least two clear avenues for further research. It would be useful to continue this analysis over time in Botswana using a random sample from the adult population, as well as measure the enrollment of respondents in the Masa program, since our cross-sectional data does not allow us to make conclusions about trends in absenteeism or earnings over time and control for participation in the national program. Additionally, if absenteeism and earnings are impacted by HIV-related illness, which is of increasing importance as HIV becomes a chronic rather than terminal condition, effects on household coping mechanisms should be studied more intensively. Individuals or households that manage small enterprises, especially those engaged in the informal sector, may bear a heavier burden if the person with HIV or AIDS is the key worker. Some households may borrow, or receive cash or in-kind donations from other households while others may pull children out of school in order for them to work to bring in needed income. Other individuals in the household may be compelled to contribute, either taking up longer work hours or additional household responsibilities [18]. This information would have important implications for the types of policies and interventions that would be designed to ease financial shocks among these households.
6. CONCLUSION
Our findings indicate that HIV-related illness increases absenteeism among working age adults in Botswana, and reduces their earning potential. We find substantial differences in these effects between men and women. Our research is the first to our knowledge highlighting the high level of HIV-related absenteeism and reduced earnings using a random sample of the adult population, and rigorous quasi-experimental methods. These results have potential implications for employer practices such as granting paid sick leave, as well as government policies surrounding welfare and disability grants to people living with HIV or AIDS in Botswana. It also may indicate that individuals who present late into the course of their disease for testing and treatment are dramatically and negatively impacting their prospects for continued employment. Botswana spends more than any country per capita on HIV, an amount totaling about 70 million USD annually [59], and the Masa program represents the most extensive of its kind relative to need [8]. While this program protects against short-term catastrophic spending, the length of AIDS illness demands that attention be paid to the indirect consequences such as loss of earnings as they relate to the wellbeing of individuals supported by working aged adults, as well as the country’s labor market.
7. AUTHOR CONTRIBUTIONS
RM and AM led the conceptualization and design of the study, with the assistance of MH, GM and CM. MF and DR carried out the analyses and interpretation of the data. MF and DR wrote the first draft of the manuscript. RM, AM, MH, GM and CM assisted in editing and revising the manuscript. All authors have approved of the final version of the manuscript.
8. ACKNOWLEDGEMENTS
This research was supported by the Global Infectious Diseases (GID) Program’s AIDS Innovative Research Fund at the Harvard Initiative for Global Health.
NOTES