KwaZulu-Natal in South Africa has one of the highest HIV prevalences in the world. Following the roll-out of HIV care and antiretroviral treatment (ART) in 2004, HIV-associated mortality decreased rapidly , leading to a substantial rise in life expectancy from 49.1 years in 2003 to 60.5 years in 2011 in our demographic surveillance site (DSS) in rural KwaZulu-Natal . This rise in life expectancy is probably largely driven by the fact that HIV-positive adults in South Africa can have a near-normal life expectancy, provided that they start ART early enough and adhere to treatment . Whereas in the initial phases of the ART programme, we observed little change in HIV prevalence and incidence locally , more recently, we observed both an increase of HIV prevalence in older age groups  and a decrease in risk of HIV acquisition with high ART coverage .
In the most recent National Strategic Plan, the South African National AIDS Council and the South African government confirm their strong commitment to reduce HIV incidence (the rate of new HIV infections) by 50% between 2012 and 2016 . However, unlike HIV prevalence, which can be measured from nationally representative cross-sectional serological surveys, measurement of HIV incidence with any degree of accuracy and reliability is much more challenging . Direct estimation of HIV incidence requires a cohort of individuals to be followed up prospectively and tested multiple times. Not only are such multiple measurements on the same individuals logistically difficult to achieve at a regional or national scale, but such data are also potentially subject to biases due to low consent rates for repeat testing  and loss to follow up due to migration. Although much effort has gone into trying to develop and validate novel biological assays that can reliably predict recent HIV infection [8,10–13], such assays were often found to require local calibration. An alternative approach to direct incidence estimation is to use a ‘catalytic’ modelling technique pioneered in the 1950s , which can potentially be applied to any age-stratified seroprevalence data sets, as long as sample sizes are large enough over wide age ranges.
Although a modelling approach has been previously applied to data from other sub-Saharan countries [15–23], we apply it here for the first time to multiple rounds of cross-sectional HIV surveys from rural KwaZulu-Natal, and for the first time in the context of a large-scale ART roll-out. Moreover, whereas previous work often relied on using external sources for HIV-associated mortality , we estimate differential HIV survival in parallel to HIV incidence using seroprevalence data only.
The aim of our study was thus to infer trends of HIV incidence and survival in rural KwaZulu-Natal from HIV prevalence data collected on a yearly basis between 2003 and 2011. This was achieved by applying a catalytic modelling approach on age-specific and sex-specific HIV seroprevalence patterns during three calendar periods (2003–2005, 2006–2008, 2009–2011), which coincide with increasing ART coverage. We used the model to simultaneously infer HIV incidence and survival during the three calendar periods and validate model-based incidence estimates with directly observed cohort incidence estimates obtained in the same setting between 2003 and 2011.
Materials and methods
Since 2003, consecutive annual HIV sero-surveys have been conducted in the DSS in a rural subdistrict of uMkhanyakude in northern KwaZulu-Natal with a resident population of approximately 60 000 Zulu-speaking people. Although logistic details of the HIV sero-surveys have been described elsewhere [24,25], briefly, all adult persons eligible to participate on the basis of age and residency status are attempted to be contacted by a field team for participation in an HIV survey. During the interview, written informed consent is obtained before collecting information on health, sexual behaviour and biomedical measures including a finger prick sample for HIV assessment. In order to ensure comparability over time and across sexes, the data used in this analysis are restricted to persons aged 15–49 years who participated in the HIV surveillance in 2003–2011, as women aged 50 years or older and men aged 55 years or older were only included in HIV surveillance from 2007 onwards. Ethical approval for the HIV surveys and the Africa Centre Demographic Surveillance System was obtained from the Research Ethics Committee at the College of Health Sciences, University of KwaZulu-Natal.
Mathematical model-based estimates of incidence and survival
We apply a modified cumulative incidence and survival model [16,18] to estimate simultaneously age-dependent HIV incidence and differential mortality in HIV-infected persons. We evaluate trends by estimating model parameters (see below) separately for women and men during three calendar periods (2003–2005, 2006–2008, 2009–2011), which coincide with increasing local ART coverage by the local HIV treatment and care programme, which was a partnership between the Department of Health and the Africa Centre. The total number of patients initiating ART in the six clinics serving the DSS increased from 372 in 2003–2005 to 3528 in 2006–2008 and attained 5573 in 2009–2011 (C. Newell, personal communication). During 2011, it is estimated that 30.7% of all HIV-infected adults aged 15–49 years in the DSS were receiving ART .
We adapt the catalytic model of Gregson et al., first, by using a lognormal functional form of incidence as proposed in Williams et al. and, second, by fitting one survival parameter and three incidence parameters to yearly age-specific HIV-seroprevalence profiles of persons aged 15–49 years assuming homogeneity over time, grouped within three calendar periods.
The incidence function f(a) (i.e. the age-dependent risk of acquiring HIV infection for a susceptible person of age a years) is assumed to follow a scaled offset lognormal distribution function  with parameters mean μ, standard deviation σ and scale θ and an offset corresponding to sexual debut fixed at 12 years:
The probability of survival s(x) after a duration of HIV infection of x years was assumed to follow the Weibull distribution,
where α and β are the shape and scale parameters, respectively. Unlike previous authors who fixed shape and scale parameters on the basis of external sources [26–28], we attempt to estimate the scale parameter β simultaneously from the seroprevalence data. We further conduct a model sensitivity analysis by varying the shape parameter α in a plausible range from 2 to 5 to see the effect of the shape parameter on incidence and survival summary statistics. Model comparisons were performed using Akaike information criterion (AIC). The results presented in the abstract and figures are based on models with the lowest AIC in the plausible range.
We obtain maximum likelihood estimates of the four parameters β, μ, σ, θ separately for men and women during three calendar periods C1 (2003–2005), C2 (2006–2008) and C3 (2009–2011) by minimizing the negative binomial log likelihood :
where Ra(t) and Na(t) are the observed number of HIV-positive individuals and total number of individuals who participated in HIV surveillance, respectively, having age a (years) during calendar year t, Cj is the calendar period [j = 1,2,3 and C1=(2003,2004,2005), C2=(2006,2007,2008), C3=(2009,2010,2011)] and is the estimated HIV prevalence at age a given as in ref. ,
The functions f and s in Eq. (4) correspond to the functions defined in Eqs. (1) and (2). The numerator of Eq. (4) represents individuals infected with HIV at age a, which is the integral of those initially infected at an earlier age who have survived to age a. The denominator represents the sum of the expression in the numerator and the individuals who have not yet been infected by age a, so the denominator encompasses all individuals who have survived to age a.
The HIV prevalence data used for fitting the incidence and survival model thus consist of observed number of HIV-positive and number of participants in HIV surveillance stratified by year of age, by year of surveillance and by sex (see Supplemental Digital Content 1, https://links.lww.com/QAD/A372).
Optimization was performed using the default Nelder-Mead method of the mle function of the package stats4 in R . Differences in parameter estimates between calendar periods were assessed assuming that parameter estimates are (asymptotically) normally distributed  and that the standard error of the difference of two parameter estimates is the sum of the squares of the standard errors .
Observed cohort incidence
Observed cohort incidence was estimated from resident persons in the demographic surveillance area who participated at least twice in the annual HIV sero-surveys, whose earliest HIV result was negative and who were aged between 15 and 49 years at the time of the survey . A uniformly distributed random seroconversion date was assigned between the last negative test and the first positive test. Incidence was calculated as the number of seroconversions in a given age and sex group divided by the total number of person-years of observations in that age and sex group between 2003 and 2011. Estimates and 95% confidence intervals (CIs) were obtained with Stata 11 (StataCorp, College Station, Texas, USA) using stptime.
The age-specific observed and fitted model HIV prevalence estimates in women and men during the three calendar periods 2003–2005, 2006–2008 and 2009–2011 in Fig. 1 indicate a good visual model fit with observed data. The higher dispersion of the observed prevalence in older age groups (most pronounced among men) is indicative of lower numbers of HIV surveillance participants in these age groups (see Supplemental Digital Content 1, https://links.lww.com/QAD/A372). This lower sample size is both a reflection of smaller numbers of older men living in the DSS (see Fig. 2 in Tanser et al. ) and lower participation rates of older men in HIV surveillance .
When the HIV prevalence model fits of the three calendar periods are compared (Fig. 2), we observe virtually no change of HIV prevalence in younger women aged 15–24 years, but a substantial increase of HIV prevalence in women aged 30–49 years in 2009–2011 compared with the two earlier periods. In men, HIV prevalence in 2009–2011 is also higher than earlier periods, but the increase occurs at a later age than in women. During all three calendar periods, HIV prevalence in women is generally higher than in men throughout all age groups, with the exception of 45–49 year olds.
A comparison between fitted model-based and observed cohort-based incidence estimates for the full calendar period 2003–2011 is shown in Fig. 3. For women, the model incidence estimates are in broad agreement and overlap with the CIs of the observed cohort estimates, particularly for the two later calendar periods 2006–2008 and 2009–2011. In men, except for teenagers, the model-based incidence estimates tend to be considerably higher than the observed cohort estimates during the periods 2003–2005 and 2006–2008. For the latest period 2009–2011, the agreement between the two methods is somewhat better with an overlap of CIs occurring for four of the seven observed cohort age group estimates.
In women, incidence was marginally lower in the two later calendar periods over all age ranges, although this was not statistically different (P = 0.932 for parameter θ in Table 1). The age of peak incidence in women varied very little between calendar periods ranging between 21.4 and 22.1 years, with peak incidence ranging from 0.064 to 0.075 per year. Similarly in men, the age of peak incidence varied little ranging from 27.3 to 29.4 years, with peak incidence ranging from 0.043 to 0.057 per year (see Table 1). The age difference between men and women at peak incidence ranged between 5.5 and 7.3 years during the three calendar periods.
The cumulative incidence by age 50 years, which is equivalent to the probability of having acquired HIV by the age of 50 (assuming current age-specific incidence and in the absence of mortality), varied very little by calendar period, ranging from 0.713 to 0.768 in women and from 0.663 to 0.717 in men (Table 1 and Fig. 4). Although incidence in younger men was lower than in women, a higher incidence during middle ages translates into a similar cumulative incidence by the time men turn 50 years old.
For women, median survival after HIV infection was estimated at 10.0 (95% CI 8.8–11.2), 12.6 (95% CI 10.9–14.3) and 14.2 (95% CI 12.6–15.8) years during the calendar periods 2003–2005, 2006–2008 and 2009–11, respectively (Table 1). The difference of the survival scale parameter β between the first and last calendar period was highly significant (P < 0.001). For men, the respective median survival after HIV infection was estimated at 10.0 (95% CI 9.2–10.8), 9.8 (95% CI 6.4–13.2) and 14.0 (95% CI 10.6–17.4). The difference of the survival scale parameter β between the first and last calendar was borderline significant (P = 0.03).
Table 2 provides details of a sensitivity analysis of model estimates when varying the shape parameter α from 2 to 5. Higher values of α make the survival distribution more rectangular-shaped by lowering the variance. Although previous authors have tended to use values ranging from 2 to 2.5 on the basis of observed mortality in sero-converters prior to ART , our results suggest that a value of α equal to 5 provided the best fit for our data. Although model fit according to AIC could be improved by increasing the value of α beyond 5 (data not shown), higher values of α had a little impact on estimates of median survival, peak incidence and age at peak incidence. In particular, the increase of median survival after HIV infection between calendar periods 2003–2005 and 2009–2011 remains robust throughout a realistic parameter range.
In this study, we use a catalytic model to provide useful insights on the evolution of the South African HIV epidemic following the large-scale roll-out of antiviral treatment. Our modelling approach is appealing because it is solely based on age-specific seroprevalence data. Such information is much more commonly available than HIV incidence derived from cohorts followed up over time. Although other authors have advocated synthetic cohort approaches on serial seroprevalence surveys [15,33], our substantially simpler approach provides reasonable HIV incidence estimates when validated against the ‘gold standard’ of the cohort methodology. However, even if the observed cohort and model-based incidence estimates are similar, one should be careful not to overinterpret results because there are potential selection biases in the data used for the two incidence estimation approaches and the biases could differ across the two approaches, for example due to lower consent rates for repeat testing . However, it is important to note that recent work [5,34] has confirmed that survey nonparticipation does not lead to large biases in cross-sectional HIV prevalence estimation in this community. It would be interesting to apply our model to other age-specific HIV prevalence data sets in high prevalence settings to obtain incidence estimates that are impossible to obtain otherwise.
Our model further provides additional evidence of a survival benefit of HIV-infected persons following ART roll-out between 2003 and 2011 in a high HIV prevalence setting in rural South Africa. Our findings, which are solely derived from cross-sectional seroprevalence data, are in agreement with previous epidemiological findings from longitudinal demographic surveillance in our setting [1,35,36] showing a significant reduction of HIV-associated mortality and increased life expectancy . The increase of survival occurred at the same time as ART coverage in the area was scaled up by the South African government with the help of PEPFAR funding . One of the side effects of this increased survival is a noticeable increase of HIV prevalence at the population level, particularly in women aged 30–49 years and in men aged 45–49 years. This in agreement with recent findings from our DSS showing that the increase in HIV prevalence was largely driven by the older age groups in which the largest proportions of HIV-infected people received ART .
Our study confirms previous findings of age-mixing patterns in this [6,9,37] and other sub-Saharan settings , particularly in that the peak of HIV incidence occurs approximately 5–7 years earlier in women than in men. In this rural setting in KwaZulu-Natal, recent work has suggested that age-disparity within a sexual partnership (i.e. ‘sugar-daddies’) did not increase the risk of HIV acquisition among young women (Harling et al., unpublished observation). Although our model allowed incidence functions to vary between sexes, for future work, it would be interesting to extend our modelling framework to directly estimate sexual mixing patterns (‘who acquires infection from whom’), which was a popular approach for respiratory infections [39,40] prior to the availability of empirical social contact data .
From a public health perspective, our modelled incidence estimates were not suggestive of a major decrease of HIV incidence, which is one of the major goals set by the South African government in the years to come . One possible reason for this is that our model could have been underpowered to detect a small reduction of incidence . A review of several transmission models has predicted that an incidence decrease in South Africa is likely to fall in the range of 15–25% by 2010 following the roll-out of ART . Moreover, a recent study  in our setting has shown that treatment levels have to reach a relatively high coverage (between 20 and 30% of all HIV-positive individuals on treatment) to have a significant effect on reducing the risk of acquiring HIV infection at the population level. Coverage levels of 30% in adults aged between 15 and 49 years were only obtained in 2011 . It should be noted that the study by Tanser et al., which also included individuals older than 50 years, found that the decline in the risk of acquisition of HIV infection with increasing ART coverage was more pronounced in the older age groups (>35 years of age). Taken together, this might suggest that incidence reductions are more likely to have occurred in older age groups due to higher access and uptake of ART.
One of the innovative aspects of our modelling approach is that we were able to estimate the increase of survival after HIV acquisition without having to rely on external mortality sources. Although previous authors have fixed both shape and scale parameters on the basis of observed survival in HIV-infected cohorts, we were able to estimate one of these, the scale parameter. Although it is possible in principle to simultaneously fit the shape and scale parameters, in practice, we found that the maximum likelihood estimates of the shape parameter tended to infinity whereas the likelihood converged to a finite value. This could be due to insufficient sample sizes at older ages, so that our data did not contain enough information to successfully estimate all parameters. The problem is likely to be amplified by correlation among parameters, which is why we decided to keep the shape parameter constant and estimate the remaining parameters as other researchers have done.
One limitation of our study is that our catalytic model does not take advantage of the longitudinal set-up of our HIV-surveillance, but rather treats yearly cross-sectional sero-surveys as independent data points. Obviously, this is a crude simplification and other authors have proposed more complicated models with a demographically correlated period/cohort or serial approaches [15,44]. However, we do take advantage of the longitudinal set-up in our model validation exercise. For future work, it would be of great benefit to refine model estimates within a Bayesian framework by integrating other existing data sources such as age-specific HIV-associated mortality that can provide additional information on the underlying relevant age-dependent epidemiological processes.
Our modelling study thus provides further evidence of the substantial demographic impact of ART scale-up between 2004 and 2011 in our rural community in South Africa in terms of improving survival of HIV-infected persons. If current trends of local ART coverage are maintained or increased, we will expect survival of HIV-infected persons to increase further, as predicted by mathematical modelling . Because our modelling approach is based on cross-sectional seroprevalence data rather than longitudinal cohort data, it could prove useful to monitor trends of HIV survival and incidence in other African settings with increasing ART coverage.
J.M., F.T., T.B. and M.L. designed the study; J.M. and E.G. coded the catalytic model in R and conducted the analysis; E.G. and F.T. provided the cohort-based incidence estimates; all authors contributed to the drafting of the manuscript.
This work was supported by the Wellcome Trust (core grant #050534), by the National Institute of Child Health and Human Development (1R01-HD058482-01) and by the German government through CIM (Centre for International Migration and Development). The authors are grateful for the help of Colin Newell in providing data. We also thank all community participants, data collectors and processors, Africa Centre virology laboratory and the support staff of the Africa Centre.
Conflicts of interest
There are no conflicts of interest.
1. Herbst AJ, Mafojane T, Newell ML. Verbal autopsy-based cause-specific mortality trends in rural KwaZulu-Natal, South Africa, 2000–2009
. Popul Health Metr
2. Bor J, Herbst AJ, Newell ML, Bärnighausen T. Increases in adult life expectancy in rural South Africa: valuing the scale-up of HIV treatment
3. Johnson LF, Mossong J, Dorrington RE, Schomaker M, Hoffmann CJ, Keiser O, et al. Life expectancies of South african adults starting antiretroviral treatment: collaborative analysis of cohort studies
. PLoS Med
4. Bärnighausen T, Tanser F, Newell ML. Lack of a decline in HIV incidence in a rural community with high HIV prevalence in South Africa, 2003–2007
. AIDS Res Hum Retroviruses
5. Zaidi J, Grapsa E, Tanser F, Newell ML, Bärnighausen T. Dramatic increases in HIV prevalence after scale-up of antiretroviral treatment
2013; [Epub ahead of print].
6. Tanser F, Bärnighausen T, Grapsa E, Zaidi J, Newell ML. High coverage of ART associated with decline in risk of HIV acquisition in rural KwaZulu-Natal, South Africa
7. National Strategic Plan for HIV, STIs and TB (2012–2016). Pretoria, South Africa: South African National AIDS Council; 2011; http://www.info.gov.za/view/DownloadFileAction?id=155622
(Accessed 11 July 2013).
8. Brookmeyer R. Measuring the HIV/AIDS epidemic: approaches and challenges
. Epidemiol Rev
9. Bärnighausen T, Tanser F, Gqwede Z, Mbizana C, Herbst K, Newell ML. High HIV incidence in a community with high HIV prevalence in rural South Africa: findings from a prospective population-based study
10. Bärnighausen T, McWalter TA, Rosner Z, Newell ML, Welte A. HIV incidence estimation using the BED capture enzyme immunoassay: systematic review and sensitivity analysis
11. Bärnighausen T, Wallrauch C, Welte A, McWalter TA, Mbizana N, Viljoen J, et al. HIV incidence in rural South Africa: comparison of estimates from longitudinal surveillance and cross-sectional cBED assay testing
. PLoS One
12. Hallett TB, Ghys P, Bärnighausen T, Yan P, Garnett GP. Errors in ‘BED’-derived estimates of HIV incidence will vary by place, time and age
. PLoS One
13. Hallett TB. Estimating the HIV incidence rate: recent and future developments
. Curr Opin HIV AIDS
14. Muench H. Catalytic models in epidemiology
. Cambridge, MA:Harvard University Press; 1959.
15. Hallett TB, Zaba B, Todd J, Lopman B, Mwita W, Biraro S, et al. Estimating incidence from prevalence in generalised HIV epidemics: methods and validation
. PLoS Med
16. Gregson S, Donnelly CA, Parker CG, Anderson RM. Demographic approaches to the estimation of incidence of HIV-1 infection among adults from age-specific prevalence data in stable endemic conditions
17. Sakarovitch C, Alioum A, Ekouevi DK, Msellati P, Leroy V, Dabis F. Estimating incidence of HIV infection in childbearing age African women using serial prevalence data from antenatal clinics
. Stat Med
18. Williams B, Gouws E, Wilkinson D, Karim SA. Estimating HIV incidence rates from age prevalence data in epidemic situations
. Stat Med
19. White RG, Vynnycky E, Glynn JR, Crampin AC, Jahn A, Mwaungulu F, et al. HIV epidemic trend and antiretroviral treatment need in Karonga District, Malawi
. Epidemiol Infect
20. Saidel T, Sokal D, Rice J, Buzingo T, Hassig S. Validation of a method to estimate age-specific human immunodeficiency virus (HIV) incidence rates in developing countries using population-based seroprevalence data
. Am J Epidemiol
21. Fontanet AL, Messele T, Dejene A, Enquselassie F, Abebe A, Cutts FT, et al. Age- and sex-specific HIV-1 prevalence in the urban community setting of Addis Ababa, Ethiopia
22. Kim AA, Hallett T, Stover J, Gouws E, Musinguzi J, Mureithi PK, et al. Estimating HIV incidence among adults in Kenya and Uganda: a systematic comparison of multiple methods
. PLoS One
23. Hallett TB, Stover J, Mishra V, Ghys PD, Gregson S, Boerma T. Estimates of HIV incidence from household-based prevalence surveys
24. Welz T, Hosegood V, Jaffar S, Batzing-Feigenbaum J, Herbst K, Newell ML. Continued very high prevalence of HIV infection in rural KwaZulu-Natal, South Africa: a population-based longitudinal study
25. Tanser F, Hosegood V, Bärnighausen T, Herbst K, Nyirenda M, Muhwava W, et al. Cohort profile: Africa Centre Demographic Information System (ACDIS) and population-based HIV survey
. Int J Epidemiol
26. Todd J, Glynn JR, Marston M, Lutalo T, Biraro S, Mwita W, et al. Time from HIV seroconversion to death: a collaborative analysis of eight studies in six low and middle-income countries before highly active antiretroviral therapy
2007; 21 (Suppl 6):S55–S63.
27. Lopman B, Gregson S. When did HIV incidence peak in Harare, Zimbabwe? Back-calculation from mortality statistics
. PLoS One
28. Zaba B, Marston M, Crampin AC, Isingo R, Biraro S, Bärnighausen T, et al. Age-specific mortality patterns in HIV-infected individuals: a comparative analysis of African community study data
2007; 21 (Suppl 6):S87–S96.
29. Grenfell BT, Anderson RM. The estimation of age-related rates of infection from case notifications and serological data
. J Hyg (Lond)
30. R Development Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2011.
31. Podgor MJ, Leske MC. Estimating incidence from age-specific prevalence for irreversible diseases with differential mortality
. Stat Med
32. Altman DG, Bland JM. Interaction revisited: the difference between two estimates
33. Rehle TM, Hallett TB, Shisana O, Pillay-van Wyk V, Zuma K, Carrara H, et al. A decline in new HIV infections in South Africa: estimating HIV incidence from three national HIV surveys in 2002, 2005 and 2008
. PLoS One
34. McGovern M, Newell ML, Tanser F, Canning D, Bärnighausen T. Accounting for selection effect in HIV prevalence estimation in the era of expanding ART access. Working Paper
. Mtubatuba:South Africa: Africa Centre for Health and Population Studies; 2013.
35. Herbst AJ, Cooke GS, Bärnighausen T, KanyKany A, Tanser F, Newell ML. Adult mortality and antiretroviral treatment roll-out in rural KwaZulu-Natal, South Africa
. Bull World Health Organ
36. Ndirangu J, Newell ML, Tanser F, Herbst AJ, Bland R. Decline in early life mortality in a high HIV prevalence rural area of South Africa: evidence of HIV prevention or treatment impact?
37. Ott MQ, Bärnighausen T, Tanser F, Lurie MN, Newell ML. Age-gaps in sexual partnerships: seeing beyond ‘sugar daddies’
38. Gregson S, Nyamukapa CA, Garnett GP, Mason PR, Zhuwau T, Carael M, et al. Sexual mixing patterns and sex-differentials in teenage exposure to HIV infection in rural Zimbabwe
39. Farrington CP, Whitaker HJ. Contact surface models for infectious diseases: estimation from serologic survey data
. J Am Stat Assoc
40. Anderson RM, May RM. Infectious diseases of humans. Dynamics and control
. Oxford:Oxford University Press; 1991.
41. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases
. PLoS Med
42. Brookmeyer R, Konikoff J. Statistical considerations in determining HIV incidence from changes in HIV prevalence
. Stat Commun Infect Dis
43. Eaton JW, Johnson LF, Salomon JA, Bärnighausen T, Bendavid E, Bershteyn A, et al. HIV treatment as prevention: systematic comparison of mathematical models of the potential impact of antiretroviral therapy on HIV incidence in South Africa
. PLoS Med
44. Keiding N. Age-specific incidence and prevalence: a statistical perspective
. J Roy Stat Soc Series A (Statistics in Society)