Estimated mortality on HIV treatment among active patients and patients lost to follow-up in 4 provinces of Zambia: Findings from a multistage sampling-based survey

To improve estimates of mortality for patients who initiate HIV treatment, Charles Holmes and colleagues use a multistage sampling approach to find and trace health outcomes for patients lost to follow-up after starting antiretroviral therapy in Zambia.

Published in the journal: . PLoS Med 15(1): e32767. doi:10.1371/journal.pmed.1002489
Category: Research Article
doi: 10.1371/journal.pmed.1002489


To improve estimates of mortality for patients who initiate HIV treatment, Charles Holmes and colleagues use a multistage sampling approach to find and trace health outcomes for patients lost to follow-up after starting antiretroviral therapy in Zambia.


As the global response to HIV moves towards universal treatment, with the intent of ending AIDS as a public health crisis by 2030, a rigorous understanding of contemporary program effectiveness demonstrates progress as well as identifies important opportunities for improvement [1]. Although other metrics also contain important information, mortality captures program effectiveness better than any other single process or clinical outcome. For example, although plasma HIV RNA suppression is a critical measure of HIV control and uncontrolled viremia indicates immunological deterioration, suppressed HIV RNA does not reflect benefits from other aspects of care such as isoniazid preventative therapy, management of side effects, or treatment of opportunistic infections. Of note, even 90-90-90 [2]—a centrally important set of metrics of national program performance—does not include monitoring of a program’s success in saving lives.

Despite the importance of mortality, relatively little data exist about long-term survival after initiation of HIV treatment in lower- and middle-income countries [3]. In South Africa, a new National Population Register for recording vital status has enabled an advanced understanding of mortality, but most high-burden settings lack such registries [4,5]. Routine monitoring and evaluation of treatment programs also fail to capture most deaths because they often happen outside of the health system in patients considered “lost to follow up” [6,7]. A nomogram has been proposed that incorporates the results of some of these tracing studies and uses the fraction of individuals lost to assess mortality, but it is not always highly accurate at site-level estimates [8]. Efforts have been undertaken to understand the “true” mortality rates through tracing random samples of patients lost to follow-up in convenience samples of clinics in multiple studies [911], but this approach has only been examined in the context of a limited number of conveniently sampled facilities. Mathematical models for the Joint United Nations Programme on HIV/AIDS (UNAIDS) and others report on mortality but depend on a number of assumptions and limited underlying epidemiological data [12].

In order to rigorously assess mortality, we apply for the first time, to our knowledge, a multistage sampling design to obtain representative estimates of mortality in 4 provinces in Zambia as well as site-level estimates with enough precision to identify site-to-site variation. This approach offers a rigorous estimate of (1) the magnitude of deaths after starting antiretroviral therapy (ART) among those already on treatment at the start of the study and among new treatment initiates, (2) when deaths occur, (3) which groups are at highest risk of death, and (4) whether these factors differ by region, facility, or other important factors such as prior healthcare utilization patterns. Both the approach and findings have implications for assessment of program effectiveness and program improvement efforts.


The protocol and study were approved by the University of Zambia Biomedical Research Ethics Committee (UNZABREC) (004-06-14) and the IRB of the University of Alabama, Birmingham, School of Medicine (F160122006), and the submitted protocol is available in the supplemental material. In this paper, we used a number of analytic innovations not specified in the analysis protocol, including the use of a modified Lorenz curve to depict the distribution of deaths across facilities and survival analyses of mortality in which deaths are classified by retention histories. In addition, analyses of pretreatment experience as specified in the protocol are also planned but not included in this paper.

Population and sampling

The target population for this analysis is the contemporary population of HIV-infected adults on ART in Zambia. We considered “contemporary” to be patients who had a clinical encounter while on ART (including the ART initiation visit) in the 24 months before the evaluation, and they are referred to as “all ART users” and provide the most complete picture of treatment outcomes. We also designate a subpopulation of this group as “new ART users,” i.e., those who initiated ART during the 24-month period prior to the evaluation—this group is more readily comparable to outcomes from other analyses of those starting ART. Our source population consisted of patients in 4 provinces (Lusaka, Southern, Eastern, and Western) who visited government-operated HIV treatment sites in these provinces, each of which receives technical assistance support from Centre for Infectious Disease Research in Zambia (CIDRZ), a local Zambian organization.

The analysis population was enumerated using a multistage sampling approach to yield estimates representative of the entire population of patients on treatment as well as to enable comparison, through stratification, across 4 provinces and 3 facility types (hospital, urban health center, and rural health center). We selected using probabilities proportional to size a minimum of 2 to 10 facilities from each of 12 strata defined by facility type and province for a total of 32 facilities. In each selected facility, we enumerated all adults on ART who made a visit in the previous 24 months and identified loss to follow-up from this cohort (defined as >90 days late at the time of sampling). From patients lost to follow-up, we selected a simple random sample for tracing inverse proportional to size. This approach seeks to enable both a regionally representative estimate of survival after starting ART and site-specific estimates with approximately 95% confidence estimates from 5% to 15% in a setting where mortality is estimated to be 10% [9]. We additionally used simulation to optimize the sampling strategy to achieve these aims. A diagram detailing the sampling process and tracing outcomes is presented in the supplementary materials (S1 Fig).

Procedures and measurements

As described in previous work, we used data from an electronic medical records system to enumerate the cohort and to obtain sociodemographic and clinical data. Peer health workers traced lost patients intensively in the community to ascertain vital status [13]. Tracing included in-depth review of paper and electronic medical records, phone calls, and in-person tracing in the community. Tracers used bicycles, public transportation, study vehicles, or motorcycles as needed. A minimum of 3 tracing attempts was required. If patients were found to be dead, death dates, cause of death, and location of death were solicited from family and other close contacts. The list of patients to be traced was ranked in random order and released in blocks of 25 to each facility to preserve the benefits of a random sample.


We estimated the cumulative incidence of mortality using the Kaplan Meier (KM) approach as well as simple mortality rates. In KM estimates, we treated time zero as the date of ART initiation. Patients who started ART before the observation period were left-truncated, which yields an estimate of survival after ART initiation during the 2-year period of observation of this study that is contingent on surviving into this period. This approach is analogous to life expectancy estimates. Mortality estimates were carried out for all ART users as well as stratified by new ART initiators during the 2-year observation period and those already on ART at the start of the current period. These analyses are intended to convey the totality of the contemporary experience of the patient population, including both that of those on treatment for extended periods of time and that of important subpopulations. We carried out mortality estimates both using only deaths known to the program before tracing and after incorporating findings from tracing lost patients to “correct” estimates through use of probability weights. These individual weights were combined with facility-level weights for estimates in the entire patient population.

In additional survival analyses, we used the Aalen Johansen method to estimate the cumulative incidence of death events classified by their current and previous engagement history using the reweighted data. We defined patients as having “died in care” if they had died within 30 days of their last visit; the remaining patients were classified as having “died out of care.” In addition, patients were categorized as “previously lost” if they had had at least 1 episode of being out of care for more than 90 days after an appointment in the past and returned to care prior to our study sampling and “previously late” if they had at least 1 episode of being late for an appointment by more than 30 days but less than 90 days prior to our study sampling. These categories yielded 6 engagement states for the analysis [14,15]. We examined the sociodemographic (e.g., sex and age), clinical (e.g., CD4 level at ART initiation), and health systems factors (e.g., facility size) associated with mortality through simple stratification and multivariate Poisson regression. Analytic weights inverse to the probability of missingness were used to address missing predictor side variables under the missing-at-random assumption [16]. We quantified site-to-site variability in mortality by examining the distribution of site-specific mortality estimates. In order to summarize the unequal contribution to total mortality across facilities, we used a modification of the Lorenz curve in which the y-axis displays the cumulative proportion of excess deaths and the x-axis the cumulative population proportion after adjustment for patient (e.g., sex and CD4 level at ART initiation) as well as facility factors (e.g., facility volume) [17]. Excess deaths are defined as deaths above a threshold of 1% at 2 years by facility—the approximate death rate in cohorts from Europe of stable patients used to represent “ideal” outcomes [18]. The Lorenz analysis was restricted to new ART initiators in order to minimize the effect of local, clinic-specific distribution of patient time of treatment.

We used 3 approaches to characterize deaths over time. First, we graphically displayed the retention history with all patients who died represented by a single horizontal bar. The length of the bar represents time since ART initiation, and colors represent time spent in care, late for care (>30 days and <90 days late for a visit), or lost (>90 days late for last appointment). Second, we examined mortality rates in each year after ART initiation among all patients as well as only among those patients who had made a visit in the last 2 years (and therefore survived into the current era) and tested for significance in the change over time. Third, we added time-varying covariates—(1) engaged, never late or lost; (2) engaged, but with a history of being late or lost; (3) late without a history of previous episodes of being late or lost; (4) late with a history of being lost; and (5) being lost (defined as being >90 days since last appointment)—in multivariable regression to capture the association between lapses in retention states and mortality. Based on the revised estimates, we quantified the fraction of deaths that fell into each category of visit histories and in each of the segments of time after ART initiation. We calculated the “population attributable fraction” of all deaths that occurred in each of these states—an estimate of how much mortality could be avoided if all measured characteristics were set at the lowest level of risk.


Patient characteristics

In a network of 64 HIV delivery facilities operated by the Zambian Ministry of Health and supported by CIDRZ across 4 provinces, we identified a total of 165,464 individuals on ART and who had made a visit in the 2 years before the project period (i.e., 1 August 2013 to 31 July 2015), which included 49,129 individuals newly initiating ART during that period and 116,335 individuals already on treatment at the opening of the observation period (Fig 1).

Flowchart depicting tracing outcomes among all patients on antiretroviral therapy (ART) and the subset of new ART initiators during the 2-year study period at 64 facilities.
Fig. 1. Flowchart depicting tracing outcomes among all patients on antiretroviral therapy (ART) and the subset of new ART initiators during the 2-year study period at 64 facilities.

In the entire cohort (Table 1), 64% were female, the median age was 39 years (IQR 33–46), 56% attended urban facilities, about half were from Lusaka Provence (52%), median CD4 at ART initiation was 201 cells/mm3 (IQR 111–312), and 58% were WHO stage 1 or 2 at ART initiation. New ART users were similar demographically but had a higher median CD4 level at ART initiation of 262 cells/mm3 (IQR 138–288), and 65% were WHO stage 1 or 2 at enrollment. The median time on ART was 1,142 days (IQR 390–2,139) for all ART users in the cohort and 255 days (IQR 63–407) in new ART starters. The median date of ART initiation for all users was 14 February 2012 (IQR 03 June 2009–23 January 2014) and for new ART users was 24 July 2014 (IQR 3 February 2014–23 December 2014).

Tab. 1. Descriptive characteristics of the study cohort.
Descriptive characteristics of the study cohort.
ART, antiretroviral therapy.

Naïve and corrected mortality estimates

Routine program data before the tracing exercise suggested that 1% of all patients had died. A total of 17% of all patients were lost to follow-up, and this fraction lost was highly variable across clinic sites (IQR 10%–31%). Among the 28,111 lost patients from the entire cohort, we made an attempt to trace 2,892 (10%). Tracing efforts led to updated vital status in 2,163/2,892 (75%) of patients, but ascertainment of vital status varied from 54% to 93% across all sites (S1 Table). Among patients lost to follow-up from the entire cohort, 17% (95% CI 15%–19%) had died. Among patients lost from the group newly starting ART during the 2-year observation period, 21% had died (95% CI 17%–25%). After incorporating outcomes from tracing into underlying data known to the program, we found the 2-year cumulative incidence of mortality among all ART users increased from a naive estimate of 1.9% (95% CI 1.7%–2.0%) to a sample-revised estimate of 7.0% (95% CI 5.7%–8.4%) (Fig 2). Among new ART users, sampling changed the cumulative incidence of mortality at 2 years from a naïve estimate of 2.1% (95% CI 1.8%–2.4%) to a revised estimate of 8.3% (95% CI 6.1%–10.7%) (Fig 2). Facility-level revised mortality estimates ranged widely, with an IQR from 1.6 to 2.6 deaths per 100 person-years among all patients and an IQR of 3.5 to 7.5 deaths per 100 person-years and a high of 13.4 deaths per 100 person-years among new ART initiators (Fig 2). Most facilities showed marked differences between naïve and revised estimates, with 1 site showing a 23-fold difference in mortality among all patients and another site showing a 14-fold difference for new ART users.

Fig. 2.
<b>Naïve and revised mortality estimates among all antiretroviral therapy (ART) users (<i>N</i> = 165,464) and new ART users (<i>N</i> = 49, 129) in (A) cumulative incidence of mortality among all ART users.</b> (B) Cumulative incidence of mortality among new ART users. (C) Mortality rate among all ART users. (D) Mortality rate among new ART users. PY, person-year.

The deaths during the first year on treatment occur mostly in persons who have not missed appointments. There is an increase over time in the fraction of deaths that occur in patients who are either late for a visit or who have a history of being late or lost from care (Fig 3). The cumulative incidence estimates and confidence intervals for this analysis are presented in the supplementary materials (S2 Table). Furthermore, general categories of reported causes of death were available for about half of the patients who died (S3 Table).

Proportion of patients who died after starting antiretroviral therapy (ART) by engagement status, with left truncation for individuals on ART before the study observation period (<i>N</i> = 165,464).
Fig. 3. Proportion of patients who died after starting antiretroviral therapy (ART) by engagement status, with left truncation for individuals on ART before the study observation period (N = 165,464).
Lost was defined as having been more than 90 days late for a visit, and late was defined as being 15–90 days late for a visit.

Regional and facility-level mortality rates

Mortality varied markedly across geographical areas and facilities. At the provincial level, revised estimates of mortality were higher than naïve estimates, ranging in all ART users from 3-fold higher in Eastern and Southern Province to 8-times higher in Lusaka Province and in new ART users from 3-fold higher in Eastern Province to 9-fold higher in Lusaka Province (Fig 4, S4 Table). The highest absolute mortality rates, prior to revision, were in Eastern Province (all ART users) and Southern Province (new ART users), whereas the highest revised mortality rates were in Western Province, followed by Southern, Eastern, and Lusaka Provinces, for both all ART users and new ART users. Facility-level mortality estimates were highly variable even within province and across facilities with similar median CD4 levels at ART initiation (Fig 5).

Fig. 4.
<b>Naïve and revised mortality rates by province in (A) all antiretroviral therapy (ART) users (<i>N</i> = 165,464) and (B) new ART users (<i>N</i> = 49,129).</b> PY, person-year.

Fig. 5.
<b>(A) Facility-level mortality rate by median CD4 at treatment initiation and (B) facility-level mortality estimates by province among all antiretroviral therapy (ART) users (<i>N</i> = 165,464); (C) facility-level mortality rate by median CD4 and (B) facility-level mortality rate by province in new ART users (<i>N</i> = 49,129).</b> PY, person-year.

Multivariable regression models using patient and facility characteristics at ART initiation found that male sex, CD4 level (≥200 cells/μL) at ART initiation, and higher WHO stage were associated with higher mortality in new ART users, whereas facility size and type (health center) were protective. For those already on ART at the start of observation, male sex; province; facility size, type, and setting; and disclosure of HIV status at treatment initiation were associated with mortality (Table 2).

Tab. 2. Predictors of mortality among all antiretroviral therapy (ART) users and new ART users.
Predictors of mortality among all antiretroviral therapy (ART) users and new ART users.

A modified Lorenz curve to capture how “concentrated” mortality was across facilities showed that approximately half (57%) of excess deaths occurred in 15% of the population (Fig 6). The Gini coefficient was 0.714, suggesting that 71.4% of deaths were not evenly distributed across the population.

Modified Lorenz curve depicting excess mortality in new antiretroviral therapy (ART) users (&lt;i&gt;N&lt;/i&gt; = 49,129).
Fig. 6. Modified Lorenz curve depicting excess mortality in new antiretroviral therapy (ART) users (<i>N</i> = 49,129).

Mortality in all ART users was highest during the first year of treatment at 4.6 per 100 person-years (95% CI 4.0–5.5), fell to 2.9 per 100 person-years (95% CI 2.1–3.9) in the second year, and then remained steady at approximately 1.6% per year for up to 8 years on treatment. The overall mortality rate in this group was 1.7 per 100 person-years (95% CI 1.5–2.0), and a test of difference in mortality rates by year on treatment for up to 10 years was not statistically significant (p-value = 0.46). When examined alone, patients who started ART during the observation time for this study (new ART users) displayed the highest mortality in the first 90 days of treatment, at 7.9 per 100 person-years (95% CI 6.3–10.0); this fell to 4.6 per 100 person-years from 91–180 days and then fell to 3.3 per 100 person-years (95% CI 2.0–5.8). Regression models with time-varying factors capturing changing retention status found that being late and being lost were both associated with mortality, but a history of lateness or loss among those in care did not predict mortality after adjustment for other sociodemographic and clinical factors. These associations were more pronounced among all ART users than among new ART users (Table 2).

Proportion of current deaths in windows of time on treatment

We estimated that 38.1% (95% CI 33.1%–43.4%) of mortality occurred among individuals in their first year on treatment, while the remaining 61.8% (95% CI 56.6%–66. 9%) was distributed among patients who had been on treatment for longer periods of time (Table 3). For patients newly starting treatment, 50.3% (95% CI 42.4%–58.2%) of all deaths occurred among those who had no history of missed visits, while 14.0% (95% CI 9.2%–20.7%) occurred in individuals while they were late, and 31.7% (95% CI 24.5%–39.9%) occurred among individuals while they were lost (S2 Fig). For those already on ART in the current period, 8.8% (95% CI 6.4%–12.1%) of deaths occurred among those who were always in care at the time of death, 13.7% (95% CI 9.7%–18.9%) occurred in those with a history of lateness (but who were in care at the time of death), 29.1% (95% CI 23.1%–35.9%) occurred in those who were in care but had a history of being lost, 11.7% (95% CI 8.1%–16.7%) occurred in individuals while they were late, and another 36.7% (95% CI 30.1%–43.8%) occurred among individuals while they were lost (S2 Fig).

Tab. 3. Duration on antiretroviral therapy (ART) and contribution to overall mortality among all ART users.
Duration on antiretroviral therapy (ART) and contribution to overall mortality among all ART users.


Despite the success of implementing HIV treatment in high-burden settings in Africa, we have identified high rates of mortality among those on treatment in Zambia, substantial under-reporting of deaths, and marked heterogeneity among provinces and sites. We applied for the first time, to our knowledge, a multistage sampling approach that included a sample of both sites and patients lost to follow-up in order to develop both provincially representative and site-level revised mortality estimates for an HIV program. In a sample representing over 160,000 patients across 4 provinces in Zambia, we found that 7% of all patients and 8% of new ART starters died within 2 years—from 10- to 20-fold higher than treated patients in Europe [18]. Provincial mortality rates developed from routine program data under-estimated true mortality rates in the provinces for new ART users by 3- to 9-fold, and revised mortality rates revealed the highest mortality rates for the HIV program in Western Province, followed by Southern, Eastern, and Lusaka. We identified important patient characteristics associated with short-term (e.g., CD4 level at initiation and male sex) and long-term mortality (e.g., male sex and disclosure), as others have described [9,19]. At the site level, care in a hospital and care in a nonurban setting were associated with a protective effect on mortality risk for all ART users, and the size of the facility had a statistically significant but modest protective effect. Importantly, we also found that the facility itself (regardless of province) is a critical driver of survival: even after adjustment for patients and facility characteristics, mortality across facilities ranged from 0 per 100 person years to 13.4 per 100 person years, even though the facilities in theory receive similar support (e.g., human resources and infrastructure). We also found that although mortality rates were highest in the first months of treatment, death rates among those on treatment for longer periods of time remain steady and high. None of these epidemiological observations would be possible using routinely collected death data, which differed markedly from sample-revised estimates.

Although mortality from childbirth, road accidents, and other factors are indeed higher in Africa, over 95% of the deaths ascertained in this cohort were due to “illness” rather than death from these competing causes. We believe the potential drivers of deaths include unmet needs for more advanced medical care, especially given the protective effect in our multivariable model of being treated in a hospital versus a health center, which may indicate a better capacity for advanced diagnostics and care in hospital-based clinics. Indeed, approximately half of deaths occur relatively shortly after a clinic visit. In a previous analysis, we found that the majority of patients who died after a recent visit had some opportunity for medical intervention that was missed—a more detailed assessment of the clinical encounters before death in these patients should be urgently undertaken [13]. At a time when public health approaches are focused on deintensifying care through differentiated service delivery [20,21], it should be recalled that in many cases, improving outcomes means identifying appropriate patients in whom to escalate care.

The provincially representative nature of the study enabled a new and more nuanced understanding of regional mortality among those newly starting ART and among the cross-section of all ART users in a province. First, the ratio of the revised mortality rate versus the naïve rate was highest in Lusaka Province—8-fold higher for all ART users and 9-fold higher for new ART users—compared with the other 3 provinces (ratios between 3–5). One interpretation of this is that populations in dense urban population centers such as Lusaka may be more likely to have their mortality under-reported. Lusaka Province is dominated by the capital city and transport hub of Lusaka, which has an urban population of 2.4 million people. It is conceivable that high mobility and/or fewer community linkages to health facilities result in less passive reporting of mortality than in the less urban provinces. The highest absolute revised mortality rates were in Western Province, which is a largely rural floodplain of the Zambezi River with the lowest population density of the 4 provinces, an HIV prevalence equal to that in Lusaka Province (16%), and limited economic activity other than cattle farming. Although our overall multivariable analysis suggested an association of nonurban facility setting with reduced mortality, health services in much of Western Province region are notoriously difficult to deliver and up-referrals for advanced care can be particularly challenging in the rainy season because of flooding and long distances to higher-level facilities. We believe these and potentially other factors may account for the high revised mortality rates in Western Province, although further investigation is needed at the provincial and site level. Indeed, the differences between the revised provincial mortality rates and between the ratios of the naïve and corrected morality rates in the provinces serve as starting points for both quality improvement work by provincial health offices and for further study. For example, future investigation could examine the system- and human-level causes of low levels of passive mortality reporting that we found across all 4 provinces, but particularly in Lusaka Province, and could investigate differential utilization of up-referral for sick patients by province and facility type and the association of provincial mortality rates with the quality and availability of advanced care services in the region.

Our results, drawn from a sample of patients visiting HIV treatment clinics between 2013 and 2015, indicate that even in an era of expanded eligibility for therapy and rising CD4 cell count thresholds, CD4 cell counts remain low at initiation and are a risk factor for mortality [22,23]. With clinical guidelines in Zambia now recommending universal treatment, we expect CD4 cell counts at treatment initiation will increase and mortality rates will decline overall. However, we anticipate that inconsistent healthcare utilization will continue to be a threat that could undermine potential gains. Although mortality amongst samples of lost patients has declined in recent years [11], perhaps reflecting an overall increase in CD4 cell counts at treatment initiation, our study also demonstrated that death rates remain unacceptably high even after 1 or more years of therapy, a time period when many might assume patients to be stable. An elevated risk of death, even years after initiating ART, was also seen in South African cohorts with high ascertainment of mortality through their national vital status registry. Standardized mortality ratios (SMRs) (observed deaths divided by the number of deaths expected if all patients were HIV negative) demonstrated that men in particular had elevated SMRs even after 48 months of therapy and having started therapy with a CD4 cell count >200 cells/mm3; women’s rates approached those of HIV-uninfected individuals by 48 months [4], a gender differential that has been widely noted [19]. These observations problematize a number of assumptions. First, it has been assumed that most patients on treatment for longer periods of time will be stable and can be placed in deintensified, demedicalized groups that emphasize social support and community-based treatment. Yet, high mortality rates in this group suggest that time on treatment itself cannot be a reliable marker of stability. Furthermore, it is assumed that intensive adherence support is critical only at initiation for most patients. Persistent deaths after 2 years on treatment imply that a subset of patients is poorly adherent or retained and that progression of HIV disease on treatment appears to be common. Therefore, effective support for adherence that is attuned to mitigating threats to adherence over time remains a crucial issue [24]. Previous work has suggested that adherence, as measured by medical possession ratio, is not optimal [25], and the high death rates on treatment may reflect disease progression. Finally, treatment does not fully protect individuals against tuberculosis; thus, rolling out isoniazid preventative therapy is a priority [26].

The marked heterogeneity of mortality across facilities is a crucial and novel observation made possible by our sampling scheme. The fact that mortality differs across facilities, even when higher-level systems support (e.g., supply chain, information systems, and central laboratory services) is relatively uniform, implies that the next generation of improvement efforts should be focused at the facility level and with access to improved data for better decision-making. In South Africa, where facilities are now linked to the national vital status registry, this is already possible, but in most settings in the region, it is not. Another implication of the site-to-site variability is that attention should be targeted to the high-volume, high-mortality-rate facilities because the absolute reductions in mortality are likely to be the biggest through improvements at these sites. Third, the presence of high-performing facilities suggests that strategies such as exemplars by positive deviants—in which high-performing facilities are brought in to assess and guide lower-performing facilities—should be explored [27]. Unexplained variability in small geographical areas supported by largely uniform systems spurred a generation of improvement science in the United States through the Dartmouth Atlas and other investigators and should be a crucial driver of improvement in the AIDS response as well [28,29].

Overall, our findings also imply that in Zambia and other similar settings, even if the regular use of site and subnational region-level mortality data for improvement were attempted, inaccuracy in the underlying data would undermine the utility of the approach. As we have shown, under-recording of mortality is the norm, and variability is masked under typical approaches to monitoring and evaluation, as evidenced by the differences in the ordering of mortality burden by province with use of naïve versus revised mortality rates. Although Demographic and Health Surveys (DHSs) and other door-to-door surveys provide representative data that can be used to evaluate trends in mortality, they are often not linked to sites and are not conducted with the frequency required to monitor for program performance and to measure the proximal effects of health systems interventions [30]. We propose 2 potential approaches to this issue. The first is to renew existing calls for nationwide death and birth registries, an approach successfully taken to scale in South Africa [31] and now being piloted in Zambia. Where this is not immediately possible, a scalable version of the representative multistage sampling performed in this study would be a feasible approach to incorporate into regional surveillance. Indeed, this approach is currently being adapted for broader use in Zambia, and we will shortly be publishing a cost analysis and a toolkit to encourage its broader use in the region. If efforts to expand this approach to improving our understanding of mortality are successful and supported by governments and funders, we would ideally be able to equip every site, district, and country with up-to-date access to accurate mortality data. In addition, incorporation of corrected provincial mortality rates into projection models could further improve their ability to inform resource needs and epidemic trends.


The methodology is predicated on tracing a random sample of individuals lost to care. We were unable to ascertain the vital status of about one-quarter of those lost and traced, which could have biased our results. However, we did not find differences in the sociodemographic or clinical characteristics among those who were traced and found compared with those we were unable to find. We were also not able to adjudicate causes of death and therefore are unable to distinguish clearly all-cause mortality versus HIV-related mortality. Further studies are needed to better understand clinical causes of death among HIV-infected individuals on ART in Zambia. Lastly, although the study included a comprehensive sample of Ministry of Health (MOH) ART sites in the 4 provinces, faith-based and other sites were not included in the estimates.


To our knowledge, this study is the first to apply a multistage sampling-based approach to generating provincially representative mortality rates for HIV programs in low- and middle-income countries (LMICs). In addition to demonstrating that provincial mortality is substantially under-reported and variable among HIV-infected individuals on ART, we found unexpectedly high site-to-site variability in mortality after starting HIV treatment under routine health service delivery conditions. This heterogeneity among sites suggests that service delivery has much room to improve and that greater attention to systematically generating and using corrected mortality rates is a logical next step for ongoing quality improvement for national programs and major donors. In addition, given the associations of late/missed visits with mortality, greater attention to patient-centered and differentiated care approaches that are also sensitive to the need for care intensification may have a role in not only improving retention but addressing the markedly high death rates, despite the widespread use and availability of ART.

Supporting Information

Attachment 3

Attachment 4

Attachment 5

Attachment 6


1. UNAIDS. Fast-track: ending the AIDS epidemic by 2030. Geneva: Joint United Nations Programme on HIV. AIDS. 2014.

2. UNAIDS. 90-90-90: An ambitious treatment target to help end the AIDS epidemic. UNAIDS, 2014.

3. Farahani M, Vable A, Lebelonyane R, Seipone K, Anderson M, Avalos A, et al. Outcomes of the Botswana national HIV/AIDS treatment programme from 2002 to 2010: a longitudinal analysis. Lancet Glob Health. 2014;2(1):e44–50. Epub 2014/08/12. doi: 10.1016/S2214-109X(13)70149-9 25104635.

4. Cornell M, Johnson LF, Wood R, Tanser F, Fox MP, Prozesky H, et al. Twelve-year mortality in adults initiating antiretroviral therapy in South Africa. J Int AIDS Soc. 2017;20(1):21902. Epub 2017/09/28. doi: 10.7448/IAS.20.1.21902 28953328.

5. Schomaker M, Gsponer T, Estill J, Fox M, Boulle A. Non-ignorable loss to follow-up: correcting mortality estimates based on additional outcome ascertainment. Stat Med. 2014;33(1):129–42. Epub 2013/07/23. doi: 10.1002/sim.5912 23873614; PubMed Central PMCID: PMCPMC3859810.

6. Brinkhof MW, Pujades-Rodriguez M, Egger M. Mortality of patients lost to follow-up in antiretroviral treatment programmes in resource-limited settings: systematic review and meta-analysis. PLoS ONE. 2009;4(6):e5790. doi: 10.1371/journal.pone.0005790 19495419.

7. PEPFAR. PEPFAR Monitoring, Evaluation, and Reporting (MER 2.0) Indicator Reference Guide: Version 2.1. January 2017.

8. Stover J, Johnson P, Zaba B, Zwahlen M, Dabis F, Ekpini RE. The Spectrum projection package: improvements in estimating mortality, ART needs, PMTCT impact and uncertainty bounds. Sexually Transmitted Infections. 2008;84:I24–I30. doi: 10.1136/sti.2008.029868 ISI:000257843600005. 18647862

9. Geng EH, Odeny TA, Lyamuya RE, Nakiwogga-Muwanga A, Diero L, Bwana M, et al. Estimation of mortality among HIV-infected people on antiretroviral treatment in east Africa: a sampling based approach in an observational, multisite, cohort study. The Lancet HIV. 2015. Epub 2015/10/02. doi: 10.1016/S2352-3018(15)00002-8 26424542; PubMed Central PMCID: PMCPMC4480204.

10. Yiannoutsos CT, An MW, Frangakis CE, Musick BS, Braitstein P, Wools-Kaloustian K, et al. Sampling-based approaches to improve estimation of mortality among patient dropouts: experience from a large PEPFAR-funded program in Western Kenya. PLoS ONE. 2008;3(12):e3843. doi: 10.1371/journal.pone.0003843 19048109.

11. Zurcher K, Mooser A, Anderegg N, Tymejczyk O, Couvillon MJ, Nash D, et al. Outcomes of HIV-positive patients lost to follow-up in African treatment programmes. Trop Med Int Health. 2017;22(4):375–87. Epub 2017/01/20. doi: 10.1111/tmi.12843 28102610; PubMed Central PMCID: PMCPMC5580236.

12. Egger M, Spycher BD, Sidle J, Weigel R, Geng EH, Fox MP, et al. Correcting mortality for loss to follow-up: a nomogram applied to antiretroviral treatment programmes in sub-Saharan Africa. PLoS Med. 2011;8(1):e1000390.

13. Geng EH, Bangsberg DR, Musinguzi N, Emenyonu N, Bwana MB, Yiannoutsos CT, et al. Understanding reasons for and outcomes of patients lost to follow-up in antiretroviral therapy programs in Africa through a sampling-based approach. J Acquir Immune Defic Syndr. 2010;53(3):405–11. Epub 2009/09/12. doi: 10.1097/QAI.0b013e3181b843f0 19745753.

14. Aalen OO, Johansen S. An empirical transition matrix for non-homogeneous Markov chains based on censored observations. Scandinavian Journal of Statistics. 1978:141–50.

15. Satagopan JM, Ben-Porat L, Berwick M, Robson M, Kutler D, Auerbach AD. A note on competing risks in survival data analysis. Br J Cancer. 2004;91(7):1229–35. doi: 10.1038/sj.bjc.6602102 15305188.

16. Allison PD. Missing Data. Thousand Oaks, CA: Sage Publications, Inc.; 2002.

17. Lorenz MO. Methods of measuring the concentration of wealth. Publications of the American statistical association. 1905;9(70):209–19.

18. Helleberg M, Afzal S, Kronborg G, Larsen CS, Pedersen G, Gerstoft J, et al. Mortality attributable to smoking among HIV-1–infected individuals: a nationwide, population-based cohort study. Clinical Infectious Diseases. 2013;56(5):727–34. doi: 10.1093/cid/cis933 23254417

19. Bor J, Rosen S, Chimbindi N, Haber N, Herbst K, Mutevedzi T, et al. Mass HIV Treatment and Sex Disparities in Life Expectancy: Demographic Surveillance in Rural South Africa. PLoS Med. 2015;12(11):e1001905; discussion e. Epub 2015/11/26. doi: 10.1371/journal.pmed.1001905 26599699; PubMed Central PMCID: PMCPMC4658174.

20. Grimsrud A, Bygrave H, Doherty M, Ehrenkranz P, Ellman T, Ferris R, et al. Reimagining HIV service delivery: the role of differentiated care from prevention to suppression. Journal of the International AIDS Society. 2016;19(1).

21. Duncombe C, Rosenblum S, Hellmann N, Holmes C, Wilkinson L, Biot M, et al. Reframing HIV care: putting people at the centre of antiretroviral delivery. Tropical Medicine & International Health. 2015;20(4):430–47.

22. Siedner MJ, Ng CK, Bassett IV, Katz IT, Bangsberg DR, Tsai AC. Trends in CD4 count at presentation to care and treatment initiation in sub-Saharan Africa, 2002–2013: a meta-analysis. Clinical Infectious Diseases. 2014;60(7):1120–7. doi: 10.1093/cid/ciu1137 25516189

23. Van Cutsem G, Ford N, Hildebrand K, Goemaere E, Mathee S, Abrahams M, et al. Correcting for mortality among patients lost to follow up on antiretroviral therapy in South Africa: a cohort analysis. PLoS ONE. 2011;6(2):e14684. Epub 2011/03/08. doi: 10.1371/journal.pone.0014684 21379378; PubMed Central PMCID: PMCPMC3040750.

24. Heestermans T, Browne JL, Aitken SC, Vervoort SC, Klipstein-Grobusch K. Determinants of adherence to antiretroviral therapy among HIV-positive adults in sub-Saharan Africa: a systematic review. BMJ global health. 2016;1(4):e000125. doi: 10.1136/bmjgh-2016-000125 28588979

25. Czaicki NL, Holmes CB, Sikazwe I, Bolton C, Savory T, wa Mwanza M, et al. Nonadherence to antiretroviral therapy among HIV-infected patients in Zambia is concentrated among a minority of patients and is highly variable across clinics. AIDS. 2017;31(5):689–96. doi: 10.1097/QAD.0000000000001347 28225707

26. Van Rie A, Westreich D, Sanne I. Tuberculosis in patients receiving antiretroviral treatment: incidence, risk factors and prevention strategies. Journal of acquired immune deficiency syndromes (1999). 2011;56(4):349. doi: 10.1097/QAI.0b013e3181f9fb39 20926954

27. Marsh DR, Schroeder DG, Dearden KA, Sternin J, Sternin M. The power of positive deviance. BMJ: British Medical Journal. 2004;329(7475):1177. doi: 10.1136/bmj.329.7475.1177 15539680

28. Schoenfeld AJ, Davies JM, Marafino BJ, Dean M, DeJong C, Bardach NS, et al. Variation in quality of urgent health care provided during commercial virtual visits. JAMA internal medicine. 2016;176(5):635–42. doi: 10.1001/jamainternmed.2015.8248 27042813

29. Munson J, Morden N, Goodman D, Valle L, Wennberg J. The Dartmouth atlas of Medicare prescription drug use. Lebanon, NH: The Dartmouth Institute for Health Policy and Clinical Practice; 2013. Available from:

30. Bendavid E, Holmes CB, Bhattacharya J, Miller G. HIV development assistance and adult mortality in Africa. Jama. 2012;307(19):2060–7. doi: 10.1001/jama.2012.2001 22665105

31. Boulle A, Schomaker M, May MT, Hogg RS, Shepherd BE, Monge S, et al. Mortality in patients with HIV-1 infection starting antiretroviral therapy in South Africa, Europe, or North America: a collaborative analysis of prospective studies. PLoS Med. 2014;11(9):e1001718. doi: 10.1371/journal.pmed.1001718 25203931

Interní lékařství

Článek vyšel v časopise

PLOS Medicine

2018 Číslo 1

Nejčtenější v tomto čísle

Tomuto tématu se dále věnují…


Zvyšte si kvalifikaci online z pohodlí domova

Antiseptika a prevence ve stomatologii
nový kurz
Autoři: MUDr. Ladislav Korábek, CSc., MBA

Citikolin v neuroprotekci a neuroregeneraci: od výzkumu do klinické praxe nejen očních lékařů
Autoři: MUDr. Petr Výborný, CSc., FEBO

Zánětlivá bolest zad a axiální spondylartritida – Diagnostika a referenční strategie
Autoři: MUDr. Monika Gregová, Ph.D., MUDr. Kristýna Bubová

Diagnostika a léčba deprese pro ambulantní praxi
Autoři: MUDr. Jan Hubeňák, Ph.D

Význam nemocničního alert systému v době SARS-CoV-2
Autoři: doc. MUDr. Helena Lahoda Brodská, Ph.D., prim. MUDr. Václava Adámková

Všechny kurzy
Kurzy Doporučená témata Časopisy
Zapomenuté heslo

Nemáte účet?  Registrujte se

Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.


Nemáte účet?  Registrujte se