INTRODUCTION
Antiretroviral therapy (ART) has become more widely available in resource-limited settings in recent years.1 Similar virological response and relative reductions in mortality compared with industrialized settings have been observed in adult patients.2,3 In sub-Saharan Africa, mortality is, however, high in the first months of ART because many patients start therapy with advanced infection. Patients who start ART with less advanced disease have mortality rates that are comparable to the general population in the second year of ART.3,4 Loss to follow-up (LTFU) is a common unfavorable outcome in ART programs in resource-limited settings.5,6 LTFU introduces bias in estimates of mortality at the level of the treatment program, that is, mortality in all patients starting ART at a given site. Efforts to improve program retention by tracing all or a sample of patients LTFU have been introduced in some programs, and mortality in patients LTFU was found to be substantially higher than in patients retained in the program.7
African children are disproportionately affected by the HIV epidemic: an estimated 2.1 million children worldwide lived with HIV at the end of 2007; 90% of them in sub-Saharan Africa.8 Mortality and program retention of children starting ART are not well defined at present, with substantial heterogeneity between pediatric treatment programs.9-15 We examined early mortality and LTFU of HIV-infected children starting ART in 10 sites participating in the International epidemiologic Databases to Evaluate AIDS in Southern Africa (IeDEA-SA). We estimated mortality and risk factors for death in children retained in care, and mortality at the program level, taking into account mortality among children LTFU.
PATIENTS AND METHODS
Study Population and Cohort Characteristics
The IeDEA-SA collaboration includes programs from Botswana, Malawi, Mozambique, South Africa, Zimbabwe, and Zambia (www.iedea-sa.org ). IeDEA-SA is part of IeDEA, which includes similar networks in other regions of Africa, Latin America and the Caribbean, Asia, and North America (www.hiv-iedea.org ).16 Data are collected at each site as part of routine monitoring at baseline (ART initiation) and each follow-up visit, using standardized definitions. Data from the different sites are transferred to data centers at the Universities of Cape Town, South Africa, or Bern, Switzerland, in a standardized format and merged at regular intervals. Information on each program is collected using standardized site assessments, which are completed each year using the web-based World Health Organization (WHO) Data Collector tool.17 Descriptions on program characteristics are based on the most recent site assessment in 2008 and mainly reflect the situation toward the end of the period of data collection for this study.
Eligibility Criteria and Outcomes
All treatment-naive children (age ≤16 years) who started ART after 1998 were included in the analysis. Enrollment of patients and selection into care vary between sites (see Table, Supplemental Digital Content 1, https://links.lww.com/QAI/A62 ). Small sites with less than 100 eligible children were excluded. Children were categorized into alive, dead, transferred out, and LTFU. The primary outcome was death from all causes. A child was considered LTFU if the last visit was more than 6 months before the closure date for that site, with the closure date defined as the most recent visit date of each site recorded in the database. Only children who commenced ART at least 6 months before the closure of the site databases could therefore be LTFU. Weight and height measurements were converted to age-standardized and sex-standardized z scores using WHO growth reference standards.18,19 Clinical stage and immunological category were defined according to WHO guidelines.20
Statistical Analysis
We used Weibull random effects regression models to identify risk factors for death in children remaining in care. The following variables were considered: age (<18 months, 18-59 months, 60-119 months, ≥120 months), gender, CD4 percentage (<10%, 10% to < 20%, ≥20%), WHO clinical stage, and year of starting ART. Missing values for CD4 percentage and WHO clinical stage were imputed using chained estimation with 10 imputation cycles.21 The same models were used to estimate cumulative mortality at 1 year for each program. Program-specific mortality estimates were then combined in random effects meta-analysis.
We calculated a weighted average of the mortality observed in children remaining in care and mortality estimated for children LTFU.22 We obtained estimates of mortality in children lost to programs from a meta-analysis of studies that traced patients LTFU to ascertain their vital status.7 This study showed an inverse relation between mortality among those lost to programs and the rate of LTFU (see Figure, Supplemental Digital Content 2, https://links.lww.com/QAI/A63 ). We used the mortality estimate that corresponded to the rate of LTFU observed in a program to estimate overall program-level mortality. For example, 0.8% of children were LTFU at year 1 in program C , whereas 12.3% were lost in program F . We, therefore, assumed that 63.2% and 51.9%, respectively, of children LTFU had died (see Table, Supplemental Digital Content 3, https://links.lww.com/QAI/A64 ). Ninety-five percent confidence intervals (CIs) for estimates of program-level mortality were obtained from Monte Carlo simulations.
We examined to what extent heterogeneity in program-level mortality might be explained by differences in prognostic factors at baseline. We adjusted analyses of children remaining in care for age, gender, WHO clinical stage, CD4 percentage, and starting year of ART and applied a cohort-specific correction factor (adjusted mortality estimate divided by crude estimate) to estimates and 95% CI of program-level mortality. Adjusted estimates were centered and correspond to a hypothetical population that is typical for the primary care settings included in this study (age 59.7 months, 50.6% male, 53.4% in WHO stage III, CD4 cell percentage 12.6%, and start of ART in 2005). Considering the young median age of the cohort, we used CD4 percentage in all age groups and thus avoided using 2 separate measures for immune deficiency (percentage and count). All data were analyzed using Stata statistical software version 10 (College Station, TX).
Ethical Approval
All study sites had local institutional review board or ethics committee approval to collect data and participate in IeDEA-SA.
RESULTS
Characteristics of Treatment Programs
Ten treatment programs in 4 countries (South Africa, Malawi, Mozambique, and Zimbabwe) were included in the present analysis. The figure in (Supplemental Digital Content 4, https://links.lww.com/QAI/A65 ) shows the location of sites, and the table in (Supplemental Digital Content 1, https://links.lww.com/QAI/A62 ) shows the characteristics of programs. Four sites are based in Cape Town (South Africa), 1 each in Soweto, Johannesburg, and Durban (South Africa), 1 in Lilongwe (Malawi), 1 in Maputo (Mozambique), and 1 in Harare (Zimbabwe). Two of the 10 programs were tertiary care sites, 8 were public programs, and 3 treated children only. All programs were located in urban areas and started ART according to clinical stage and CD4 percentage. At 4 sites, ART in children younger than 12 months was initiated according to revised WHO guidelines23 toward the end of the period of data collection.
Patient Characteristics
Table 1 summarizes patient characteristics at the start of ART: a total of 8225 HIV-infected children were analyzed. The 3 largest programs (sites B, D, and G) contributed 4413 (53.7%) children; 3 sites with a total of 2153 (26.2%) patients were from outside South Africa. Median age was 49 months (range across sites 16-120 months), 2045 (24.9%) were younger than 18 months (range 2.7% to 53.4%), 4159 (50.6%), were male (range 45.8% to 55.3%), and median follow-up was 17.3 months (range 8.7 to 25.7 months). Median CD4 percent at start of ART ranged from 8% to 15% and the proportion of children with WHO clinical stage III/IV from 63.7% to 96.9%. The most frequent initial ART regimen was stavudine/lamivudine in combination with efavirenz (37.4%), nevirapine (22.0%), or lopinavir/ritonavir (18.1%).
TABLE 1: Baseline Characteristics and Outcomes of Children Starting ART in Each of 10 Treatment Programs in Southern Africa
TABLE 1: (continued ) Baseline Characteristics and Outcomes of Children Starting ART in Each of 10 Treatment Programs in Southern Africa
Children treated in tertiary care programs were younger (median age 19 months versus 57 months, P < 0.0001) and more likely to be in clinical stage III/IV (88.7% versus 75.3%, P < 0.0001). Rates of transfer-out were higher for tertiary care programs: 42.7% versus 7.7% in other sites, P < 0.0001. Overall, 1176 (14.3%) had been transferred to another site (range 0%-45.6%). A total of 7482 eligible children with at least 6 months of follow-up could be included in calculations of the rate of LTFU: 717 children (9.6%) were LTFU (range across programs 1.9% to 14.2%).
Mortality at One Year in Children Remaining in Care and at the Program Level
A total of 391 children (4.8%) died in the first year of ART although in care. Figure 1 shows cumulative mortality at 1 year by treatment program. Mortality was substantially higher in the first 3 months of ART compared with later months but did not significantly change between 4-6 and 7-12 months after start of ART (Table 2 ). There was substantial heterogeneity between programs as follows: mortality ranged from 0.8% (95% CI: 0.2% to 3.1%) in site C to 16.8% (95% CI: 14.2% to 19.9%) in site 1. The combined mortality estimate for all programs from random effects meta-analysis was 4.5% (95% CI: 2.8% to 7.4%) among children remaining in care (Fig. 1 ). Mortality was higher in the tertiary care programs compared with the other sites (13.5% versus 4.1%, P = 0.014). Table 2 shows baseline variables associated with the probability of death in children remaining in care from crude and adjusted Weibull models. Younger age, lower CD4 cell, more advanced clinical stage, and starting ART in earlier calendar years were associated with higher mortality.
TABLE 2: HRs for Death in the First Year of ART in Children Remaining in Care in 10 Treatment Programs in Southern Africa
FIGURE 1: Cumulative mortality at one year in children starting ART and remaining in care in ten treatment programs in Southern Africa. Circles show estimates from individual programs, horizontal lines indicate 95% confidence intervals, and the diamond shows the combined estimate from random-effects meta-analysis with 95% confidence intervals.
Figure 2 shows program-level mortality at 1 year, after taking mortality among children LTFU into account. Mortality ranged from 1.3% (95% CI: 0.20% to 2.7%) in program C to 18.9% (95% CI: 15.9% to 22.1%) in program I. The combined estimate was 8.7% (95% CI: 5.4% to 12.1%). Finally, Figure 3 shows program-level mortality at 1 year after adjusting for age, gender, CD4 percentage, WHO clinical stage, and year of starting ART. Mortality ranged from 1.0% (95% CI: 0.2% to 2.1%) in program C to 10.2% (95% CI: 6.6% to 14.2%) in program A. The higher mortality in tertiary care clinics was clearly explained by differences in prognostic factors, but some heterogeneity persisted due to lower mortality in children treated in 2 community programs in South Africa.
FIGURE 2: Cumulative mortality at one year at the level of the program in ten treatment programs in Southern Africa. Circles show estimates from individual programs, horizontal lines indicate 95% confidence intervals, and the diamond shows the combined estimate from random-effects meta-analysis with 95% confidence intervals.
FIGURE 3: Cumulative mortality at one year at the level of the program in ten treatment programs in Southern Africa, adjusted to correspond to a hypothetical primary care population. Circles show estimates from individual programs, horizontal lines indicate 95% confidence intervals, and the diamond shows the combined estimate from random-effects meta-analysis with 95% confidence intervals.
Sensitivity Analysis
Figure 4 shows estimates of program-level mortality assuming different mortality rates among children LTFU. Taking LTFU up into account increased mortality estimates compared with mortality in children retained in care, with substantial changes being observed in programs with higher rates of LTFU. The combined estimates from random effects meta-analysis were less sensitive to mortality assumptions: assuming that 25%, 50%, or 75% of children LTFU have died resulted in combined estimates of mortality at 1 year of 7.2% (95% CI: 4.2% to 10.2%), 8.4% (95% CI: 5.1% to 11.8%), and 9.7% (95% CI: 6.1% to 13.3%), respectively.
FIGURE 4: Sensitivity analysis of mortality at one year after starting ART by treatment program, assuming mortality rates of 25%, 50% and 75% among children loss to follow-up (LTFU). Sites are ordered by increasing rates of LTFU at one year (from left to right).
DISCUSSION
This collaborative analysis of 10 pediatric treatment programs in 4 countries in Southern Africa found that 1-year mortality in children starting ART varied across treatment programs, with heterogeneity partly explained by differences in prognostic factors at the start of ART, and differences in LTFU. Our study illustrates that ignoring LTFU may lead to substantial underestimation of mortality: mortality in children remaining in care overall was below 5% but was estimated to be about twice as high at the program level, when considering deaths in children LTFU.
We examined mortality and risk factors for death in children remaining in care and then estimated program-level mortality. Previous analyses of mortality in patients starting ART have generally truncated (censored) follow-up time in patients LTFU, thus assuming that their mortality experience is comparable to similar patients remaining in care.12,15,24 This assumption is problematic: many patients who are LTFU stop taking ART and their mortality is high. In addition, patients may not return for a follow-up appointment because they have died. A meta-analysis of studies that traced adult patients LTFU found that in sub-Saharan Africa mortality in patients whose vital status could be ascertained was 46% (95% CI: 39% to 54%).7 Interestingly, there was an inverse relation between mortality among those LTFU and the rate of LTFU. We took this into account and used the mortality estimate that corresponded to the rate of LTFU observed in a given program. A recent analysis of patients LTFU in the rural Médecins sans Frontières program in Chiradzulu, Malawi, found that mortality among patients LTFU was similar for adults and children,25 thus supporting our use of estimates from studies in adults. Another study from Lilongwe, Malawi, described a mortality estimate of 33% (CI: 20.4% to 48.4%) among children LTFU (Weigel R, et al, unpublished data, 2010). The high early mortality in infants and young children means that adult estimates of mortality in those LTFU are not necessarily applicable to cohorts which include a large proportion of very young children.26,27 Our sensitivity analyses showed, however, that combined results were fairly robust to different mortality assumptions for children LTFU.
The context of LTFU in children is different and more complex compared with adults: children are dependent on caregivers, and child outcomes may be related to those of their parents. Even though caregivers are motivated to support children to adhere to ART, the caregiver may change frequently due to ill health, death, or adult migration. In addition, unlike adult ART provision, most pediatric ART services remain located in larger hospitals in urban centers, and the cost of transport and inconvenience of having to attend separate clinics for caregiver and child may contribute to LTFU.28-30 Interestingly, at one of the few family centered ART clinics, it has been shown that children whose caregiver is in HIV care at the same clinic experience better outcomes.31 The cost of pediatric drugs and difficulty of giving children large amounts of medicine which may be unpalatable are further obstacles to effectively remaining in care.28,30,32 Furthermore, many caregivers are elderly grandparents for whom these obstacles may become a substantial barrier to continuing on treatment.27,30 Some of the children LTFU may have continued treatment at another clinic, perhaps a primary health care center closer to their home, without notifying the program where ART was started (“silent” transfer-out).
Our study included children who were treated in 10 ART programs in 4 countries in Southern Africa, including pediatric programs and programs mainly treating adults, and programs at tertiary and lower levels of care. These sites may not necessarily reflect the situation of this region as a whole: all sites were located in urban areas and cohorts from South Africa with a link to a research program predominated. However, this is one of the largest pooled analyses of children on ART ever done and from one of the regions of the world most heavily affected by the HIV epidemic. Results should therefore be applicable to many other children on ART.
Few previous studies reported mortality among children starting ART: a pooled analysis from the Kids-ART-LINC collaboration on over 2000 children on ART in sub-Saharan Africa reported a probability of LTFU of 5.0% and a mortality of 6.0% at 1 year,9 which is slightly higher than the mortality estimate in children remaining in care in the present study. The number of children included was smaller than in the present study and, more importantly, the analysis did not take mortality among children LTFU into account. The estimate from the Kids-ART-LINC analysis is therefore probably an underestimate. A large program from Lusaka, Zambia, with 2398 children on ART showed a mortality rate of 6.6 deaths per 100 child-years, with an estimated cumulative mortality at 1 year (from the Kaplan-Meier curve) of 7.9%. Children LTFU were excluded from this analysis.10 A pediatric ART treatment program in Haiti indicated that 9% of children have died, and 10% were lost after a median follow-up time of 20 months.33
The definition of LTFU was uniform across sites, but the period of 6 months chosen for this definition may be too short for some programs and too long for others. In the absence of dedicated studies on the appropriate definition of LTFU, this remains speculative. Follow-up of children who were transferred out was censored, thus assuming that their prognosis is the same as that of comparable children remaining in care at the site. Some children who were transferred out might, however, have a better or worse prognosis than the children remaining in a program. The proportion of children transferred out was particularly high in tertiary care, and these children will have been transferred to a lower level of care.
We observed differences between tertiary and other care programs in baseline characteristics and a higher mortality in tertiary care sites. The 2 tertiary programs are located in a region in the Republic of South Africa with a large pediatric ART network and a well functioning referral system; hence patients who are younger and represent the severe end of the pediatric HIV disease spectrum predominate at these sites. Statistical adjustment for prognostic factors showed that if children treated in tertiary centers were similar to the children followed in the other programs, particularly with regard to age and CD4 percentage, their mortality would be similar to that observed in the other programs. Indeed, young age is an important prognostic factor as shown in a recent collaborative analysis of 33 European cohorts, which demonstrated that mortality was almost double in children aged less than 2 years, compared with older children.26 Children treated at 2 township programs in South Africa3 had lower mortality, and this was not explained by more favorable prognostic factors at baseline. These programs mainly treat adult patients, and the children followed there may represent a selected group of children with a good prognosis. Other sources of heterogeneity might include different patterns of transfer-out, different parental outcomes related to children's outcome, incorrect assumptions on mortality in children LTFU, and incomplete or absent adjustment for other prognostic variables, for example, exposure to antiretroviral drugs for preventing mother-to-child transmission. Furthermore, where programs do not have the capacity to commence ART quickly, mortality on ART may be lower due to survival bias among children who lived long enough to access treatment. The (non-HIV-related) background mortality in the different populations may also have played a role.4
Treatment programs should attempt to actively trace patients considered as LTFU where this has not yet been done and record outcomes in these patients LTFU. Where such mechanisms are already in place as in most of the sites included in this study, capacities should be increased to intensify these activities. A modeling study based on a program in Côte d'Ivoire found that interventions that enhance program retention would both improve survival and be cost-effective by international criteria.34 An additional benefit is that program-level mortality can then be estimated by using double-sampling designs to account for mortality among children LTFU.35,36 For example, data from an outreach program have been used to adjust mortality estimates from an adult treatment program in Kenya.37,38 An important assumption made by this method is that a random sample of patients LTFU was selected and that the vital status of everyone from that sample was ascertained. These assumptions are not generally met: patients to be traced tend to be selected according to their assumed place of residence, excluding patients who live far from the treatment site, and typically about a third of patients cannot be located.7 Our approach is useful for sensitivity analysis in programs where no dedicated tracing studies are available and to test the robustness of estimates based on the results of tracing activities.
In conclusion, our study indicates that mortality of HIV-infected children starting ART in Southern Africa is substantially higher than the mortality directly observed in treatment programs. A better understanding of outcomes in children LTFU and children transferred out of programs is required, and the determinants of LTFU and transfer-out need to be defined. Dedicated field studies tracing children LTFU and studies linking data from treatment programs with routine mortality data from national death registries are of high priority. It is important both from a clinical and programmatic point of view to ensure capacity to trace this vulnerable patient group as soon as possible after a missed clinic visit.
ACKNOWLEDGMENTS
We thank all children whose data were used in this study and their caregivers. We also would like to thank all who contributed to recording and entering data and preparing and sending it to the IeDEA Southern Africa collaboration.
REFERENCES
1. World Health Organization.
Towards Universal Access. Scaling Up Priority HIV/AIDS Interventions in the Health Sector. Progress Report 2008 . Available at:
http://www.who.int/hiv/pub/towards_universal_access_report_2008.pdf . Geneva, Switzerland: World Health Organization, UNAIDS, Unicef; 2008.
2. Braitstein P, Brinkhof MW, Dabis F, et al. Mortality of HIV-1-infected patients in the first year of antiretroviral therapy: comparison between low-income and high-income countries.
Lancet . 2006;367:817-824.
3. Keiser O, Orrell C, Egger M, et al. Public-health and individual approaches to antiretroviral therapy: township South Africa and Switzerland compared.
PLoS Med . 2008;5:e148.
4. Brinkhof MW, Boulle A, Weigel R, et al. Mortality of HIV-infected patients starting antiretroviral therapy in sub-Saharan Africa: comparison with HIV-unrelated mortality.
PLoS Med . 2009;6:e1000066.
5. Brinkhof MW, Dabis F, Myer L, et al. Early loss of HIV-infected patients on potent antiretroviral therapy programmes in lower-income countries.
Bull World Health Organ . 2008;86:559-567.
6. Rosen S, Fox MP, Gill CJ. Patient retention in antiretroviral therapy programs in sub-Saharan Africa: a systematic review.
PLoS Med . 2007;4:e298.
7. 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:e5790.
8. World Health Organization.Towards universal access by 2010: How WHO is working with countries to scale-up HIV prevention, treatment, care and support. Available at:
http://www.who.int/hiv/mediacenter/universal_access_progres_report_en.pdf . Geneva, Switzerland: WHO; 2007.
9. KIDS-ART-LINC Collaboration. Low risk of death, but substantial program attrition, in pediatric HIV treatment cohorts in Sub-Saharan Africa.
J Acquir Immune Defic Syndr . 2008;49:523-531.
10. Bolton-Moore C, Mubiana-Mbewe M, Cantrell RA, et al. Clinical outcomes and CD4 cell response in children receiving antiretroviral therapy at primary health care facilities in Zambia.
JAMA . 2007;298:1888-1899.
11. Kekitiinwa A, Lee KJ, Walker AS, et al. Differences in factors associated with initial growth, CD4, and viral load responses to ART in HIV-infected children in Kampala, Uganda, and the United Kingdom/Ireland.
J Acquir Immune Defic Syndr . 2008;49:384-392.
12. Sutcliffe CG, van Dijk JH, Bolton C, et al. Effectiveness of antiretroviral therapy among HIV-infected children in sub-Saharan Africa.
Lancet Infect Dis . 2008;8:477-489.
13. Bolton-Moore C, Mubiana-Mbewe M, Cantrell RA, et al. Clinical outcomes and CD4 cell response in children receiving antiretroviral therapy at primary health care facilities in Zambia.
JAMA . 2007;298:1888-1899.
14. Kekitiinwa A, Lee KJ, Walker AS, et al. Differences in factors associated with initial growth, CD4, and viral load responses to ART in HIV-infected children in Kampala, Uganda, and the United Kingdom/Ireland.
J Acquir Immune Defic Syndr . 2008;49:384-392.
15. Sutcliffe CG, van Dijk JH, Bolton C, et al. Effectiveness of antiretroviral therapy among HIV-infected children in sub-Saharan Africa.
Lancet Infect Dis . 2008;8:477-489.
16. Gange SJ, Kitahata MM, Saag MS, et al. Cohort profile: the North American AIDS Cohort Collaboration on Research and Design (NA-ACCORD).
Int J Epidemiol . 2007;36:294-301.
17. World Health Organization. Data collector. Available at:
https://extranet.who.int/datacol/WebHelp/DataCollector_OLH.htm . Accessed February 2008.
18. World Health Organization.Growth reference data for 5-19 years. Available at:
http://www.who.int/growthref/en/ . Accessed February 2008.
19. World Health Organization. WHO child growth standards. Available at:
http://www.who.int/childgrowth/standards/en/ . Accessed February 2008.
20. World Health Organization. Antiretroviral therapy of HIV infection in infants and children: towards universal access. Available at:
www.who.int/hiv/pub/guidelines/WHOpaediatric.pdf . Accessed February 2008.
21. Royston P. Multiple imputation of missing values: an update.
Stata J . 2005;5:1-14.
22. Egger M, Spycher B, Sidle J, et al. Correcting mortality for loss to follow up: a graphical approach applied to ART programmes in resource-limited settings. Presented at: 5th IAS Conference on HIV Pathogenesis, Treatment and Prevention; July 19-22, 2009; Cape Town, South Africa. Abstract WEPED173.
23. World Health Organization.WHO antiretroviral therapy for infants and children 2008. Report of the WHO technical reference group, paediatric HIV/ART care guideline group meeting. Geneva, Switzerland: WHO headquarters. Available at:
http://www.who.int/hiv/pub/meetingreports/art_meeting_april2008/en/index.html . Accessed May 2009.
24. Davies MA, Egger M, Keiser O, et al. Paediatric antiretroviral treatment programmes in sub-Saharan Africa: a review of published clinical studies.
Afr J AIDS Res . 2009;8:329-338.
25. Mukhuna B, McGuire M, Taremba K, et al. Outcomes of patient active tracing by lay community workers in a decentralized HIV care program of rural Malawi. Presented at: 5th IAS Conference on HIV Pathogenesis, Treatment, and Prevention; July 19-22, 2009; Cape Town, South Africa. Abstract WEPED195.
26. Collaboration of Observational HIV Epidemiological Research Europe (COHERE) Study Group. Response to combination antiretroviral therapy: variation by age.
AIDS . 2008;22:1463-1473.
27. Davies MA, Keiser O, Technau K et al. Outcomes of the South African National Antiretroviral Treatment Programme for children: the IeDEA Southern Africa collaboration.
S Afr Med J . 2009;99:730-737.
28. Vreeman RC, Wiehe SE, Ayaya SO, et al. Association of antiretroviral and clinic adherence with orphan status among HIV-infected children in Western Kenya.
J Acquir Immune Defic Syndr . 2008;49:163-170.
29. Eley B, Nuttall J. Antiretroviral therapy for children: challenges and opportunities.
Ann Trop Paediatr . 2007;27:1-10.
30 Meyers T, Moultrie H, Naidoo K, et al. Challenges to pediatric HIV care and treatment in South Africa.
J Infect Dis . 2007;196(Suppl 3):S474-S481.
31. Reddi A, Leeper SC, Grobler AC et al. Preliminary outcomes of a paediatric highly active antiretroviral therapy cohort from KwaZulu-Natal, South Africa.
BMC Pediatr . 2007;7:13.
32. Davies MA, Boulle A, Fakir T, et al. Adherence to antiretroviral therapy in young children in Cape Town, South Africa, measured by medication return and caregiver self-report: a prospective cohort study.
BMC Pediatr . 2008;8:34.
33. George E, Noel F, Bois G et al. Antiretroviral therapy for HIV-1-infected children in Haiti.
J Infect Dis . 2007;195:1411-1418.
34. Losina E, Toure H, Uhler LM, et al. Cost-effectiveness of preventing loss to follow-up in HIV treatment programs: a Cote d'Ivoire appraisal.
PLoS Med . 2009;6:e1000173.
35. Baker SG, Wax Y, Patterson BH. Regression analysis of grouped survival data: informative censoring and double sampling.
Biometrics . 1993;49:379-389.
36. Frangakis CE, Rubin DB. Addressing an idiosyncrasy in estimating survival curves using double sampling in the presence of self-selected right censoring.
Biometrics . 2001;57:333-342.
37. Geng EH, Emenyonu N, Bwana MB, et al. Sampling-based approach to determining outcomes of patients lost to follow-up in antiretroviral therapy scale-up programs in Africa.
JAMA . 2008;300:506-507.
38. Yiannoutsos CT, An MW, Frangakis CE, 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:e3843.
APPENDIX I: IeDEA-SA Participating Sites
Margaret Pascoe, Newlands Clinic, Harare, Zimbabwe; Robin Wood, Gugulethu ART Program, Cape Town, South Africa; Harry Moultrie, University of Witwatersrand Paediatric HIV Clinics (Harriet Shezi Clinic, Chris Hani Baragwanath Hospital), Johannesburg, South Africa; Karl Technau, University of Witwatersrand Paediatric HIV Clinics (Rahima Moosa Mother and Child Hospital), Johannesburg, South Africa; Gilles van Cutsem, Khayelitsha ART Program, Médecins sans Frontières and University of Cape Town, Cape Town, South Africa; Paula Vaz, Paediatric Day Hospital, Maputo, Mozambique; Janet Giddy, McCord Hospital, University of Kwazulu Natal, Durban, South Africa; Brian Eley, Red Cross Children's Hospital and the University of Cape Town, Cape Town, South Africa; Hans Prozesky, Tygerberg Academic Hospital, University of Stellenbosch, Cape Town, South Africa.