WHO guidelines for adult patient monitoring in 2010 recommended that patients be switched to second-line antiretroviral therapy (ART)  based on clinical criteria, immunologic criteria, or virologic criteria (viral load over 5000 copies/ml while on treatment). Viral load monitoring can be either routine or targeted; in the latter case, it is used to confirm virological failure in patients meeting either clinical or immunologic criteria, and in the former case it is performed routinely irrespective of other clinical or immunologic results. However, viral load tests are currently more expensive than CD4+ cell count tests, and the high cost of second-line therapy can drive up costs associated with viral load testing as a result of earlier and more frequent switching to second-line therapy. For this reason, the 2013 WHO guidelines overlay considerable debate about the benefits and frequency of routine viral load testing.
When allocating resources in low-resource settings, the benefits of more frequent viral load monitoring need to be balanced against its costs and disadvantages. The benefits of more frequent routine virological monitoring include timely switching to second-line therapy for those failing to suppress virus and therefore at higher risk of transmitting infection to others; avoiding unnecessary switches to second-line therapy for patients who do not need a new regimen; and limiting the development and onward transmission of virus resistant to first-line regimens. The costs and disadvantages of more frequent routine virological monitoring include those of tests, equipment, staff time, decrements in patient adherence associated with lengthy travel to sites where tests can logistically occur, and subsequent switches to more expensive second-line and third-line regimens.
It is particularly important to consider any additional costs of routine virological monitoring in light of its opportunity costs; that is, the gains that could be had by devoting resources alternatively to other simultaneously resource-constrained decisions, such as increasing the eligibility for first-line ART by raising the CD4+ cell initiation criteria, and/or by adopting more aggressive measures to pursue more complete ART coverage. These tradeoffs are complex. Mathematical modeling allows for systematic and detailed consideration of the costs and benefits of a broad potential repertoire of strategies over a range of timescales, which would be impossible to do experimentally with clinical trials.
Accordingly, our aim was to use a published validated computer simulation of the HIV epidemic in Kenya, now enhanced to incorporate HIV transmission as well as disease progression, to evaluate the cost–effectiveness of a wide spectrum of routine viral load monitoring strategies, seeking to identify the intensity of routine viral load monitoring at which its incremental cost–effectiveness ratio (ICER) became more favorable than the ICER of more expansive ART eligibility (in other words, when additional routine viral load monitoring bought more ‘health’ than using those resources, alternatively, to expand ART eligibility).
We used a computer simulation to explore alternative scenarios regarding cost, frequency, and switching threshold of alternative laboratory monitoring strategies, including routine viral load monitoring (every 6 months or every 12 months; switching thresholds of 1000 or 10 000 copies/ml), routine CD4+ cell monitoring, and routine CD4+ cell monitoring with confirmatory viral load. In addition, we explored alternative scenarios regarding ART initiation thresholds (CD4+ cell count of 200, 350, 500 cells/μl, and no CD4+ cell threshold). This model is different from our previous publications [2,3] because it considers HIV transmission as well as disease progression. For each ART initiation threshold, we sought to identify the viral load monitoring strategy at which the ICER of more frequent routine viral load testing became more favorable than the ICER of more expansive ART eligibility (in other words, when additional routine viral load monitoring bought more ‘health’ than using those resources, alternatively, to more fully comply with expanded ART eligibility). Cost inputs were based on previously reported expenditure data from a large healthcare system in western Kenya (Academic Model Providing Access to Healthcare, or AMPATH) , and disease progression and transmissibility inputs were based on prior published work. We used a discount rate of 3%. We used a payer perspective in base case analyses, exploring a societal perspective in sensitivity analyses. We used a time horizon of 20 years. Outcome measures included cost per quality-adjusted life year gained. For analyses in this article, we assumed two ART regimens were available; however, the simulation has the ability to consider alternative scenarios regarding numbers of ART regimens available.
We sought to identify ‘efficient frontiers’, defined as those strategies delivering the greatest health benefit given a particular budget scenario . Strategies within an efficient frontier confer the greatest benefit for a specified budget. Strategies outside this frontier are unable to deliver the greatest benefit regardless of budget, and therefore are not preferred choices regardless of the magnitude of available resources. We identified efficient frontiers by calculating the ICERs of HIV treatment initiation and monitoring strategies. ICERs measure the additive benefit of each strategy compared with its next best alternative, and interpret this benefit together with its additive cost.
This simulation is composed of a stochastic progression module (e.g., hypothetical patients are followed over time in a microsimulation based on a state transition model, and depending on ART adherence and other factors, may be more or less likely to die of AIDS versus other causes) that provides data to inform a deterministic transmission module (e.g., hypothetical compartments of persons interact with one another, and HIV-infected compartments may transmit the infection to non-HIV-infected groups).
Notable features of our simulation include that it can examine nearly infinite combinations of testing strategies and ART initiation scenarios, explicitly considers ART adherence, explicitly considers accumulation and transmission of resistance mutations, can perform evaluations over long time horizons (thereby capturing delayed harms that may not be evident in shorter term analyses), can partition effects to transmission versus progression, captures important heterogeneity in sexual mixing patterns, and has been well validated based on clinical and epidemiological data.
The simulation projects the course of the HIV epidemic over varying time horizons, and tracks the benefits of potential interventions using a variety of outcome measures, including life expectancy, measured in life years; quality-adjusted life expectancy, measured in quality-adjusted life years (QALYs); and HIV infections averted.
HIV progression module
Disease progression was modeled by evaluating mortality rates and trajectories of CD4+ cell counts and HIV-1 viral load within a previously described HIV progression simulation calibrated and validated on East African populations. This model explicitly represents the main causes of ART failure, nonadherence, and the accumulation of genotypic resistance, and has been well validated in multiple populations (see Table 1[2,5–29] for selected model input parameters). Details of this module are described elsewhere .
HIV transmission module
A dynamic compartmental model of HIV transmission was developed, specified by sets of differential equations. The model includes heterosexual transmission, but does not include homosexual transmission or transmission from needle-sharing during injection drug use.
At any particular time, hypothetical people in the module occupy one among a set of mutually exclusive and collectively exhaustive compartments (Fig. 1). These compartments describe health characteristics as well as behavioral risk characteristics. As time proceeds, hypothetical people may change the compartment that they occupy, for example, proceeding from being uninfected to having a primary HIV infection to having a chronic HIV infection to death. (Alternatively they may die without ever contracting HIV). In addition, groups of people in the model may alternatively be abstinent (in which case they will not contract HIV), monogamous, nonmonogamous, or may occupy an even higher risk group [community sex workers (CSWs) if women, and migrant workers if men], with correspondingly greater chances of contracting or spreading HIV. Women in any of the nonabstinent states may ‘mix’ (have sexual contact) with men in any of the nonabstinent states. Probability of transmission is a function of multiple factors, including rate of acquiring new partners, duration of partnership, frequency of sexual contact within a partnership, and likelihood of condom use. For example, CSWs are especially likely to transmit HIV to nonmonogamous men because the likelihood of initiating contact is high, condom use is low, and duration of the relationships may be high.
The HIV-infected population in the transmission simulation at baseline was divided into compartments based on CD4+ cell and viral load strata. Five CD4+ cell strata were represented (<50, 50–200, 200–350, 350–500, >500 cells/μl) and five log viral load logarithmic strata were represented (<2.5, 2.5–3.5, 3.5–4.5, 4.5–5.5, >5.5 log units/ml). Data regarding transition rates and trajectories by treatment and adherence level were captured from the progression model and interfaced with the transmission simulation in the form of rate tables. The spectrum of infection and care was modeled as a stepwise progression from HIV acquisition/primary infection to chronic infection, to HIV detection through testing, followed by linkage to care, and finally initiation of treatment with ART. Details of this module are described in Appendix A – Technical Appendix, http://links.lww.com/QAD/A428.
Other design features include the following:
- All sexually active persons may have untreated sexually transmitted infections (STIs), and this increases HIV transmission risk. Risk of having untreated STIs, like HIV itself, will increase with greater numbers of sexual contacts and with unhealthy alcohol use.
- Men may or may not be circumcised. Men who are circumcised have a decreased chance of being infected with HIV.
- Likelihood of transmission per sexual encounter will vary with viral load (greater with higher viral load). For this reason, individuals who are receiving antiretroviral treatment will have a lower risk of transmission than individuals who are not receiving antiretroviral treatment, and individuals with primary HIV infection (and a correspondingly high viral load) or AIDS will be especially likely to transmit the virus.
Whereas many transmission models assume ‘homogenous mixing’ (i.e., each hypothetical individual has an equal chance of transmitting the infection to each other hypothetical individual) or homogenous mixing stratified by age (i.e., each hypothetical individual within a particular age stratum has an equal chance of transmitting the infection to each other hypothetical individual within that age stratum), our transmission module assumes heterogeneous mixing patterns that are informed by sub-Saharan African data. Therefore, it represents the common phenomenon of assortative mixing, for example, individuals who engage in risky sex may be more likely to partner with other people who have risky sex than with people who do not engage in risky sex, even after controlling for the increase in partnering opportunities that may be available. Additionally, because older men may be more likely to have sexual contact with younger women than vice versa, the model includes a parameter that can be varied to include this age asymmetry of mixing.
Decomposition of risk per partnership
Transmission models often assume a single, aggregate risk that accumulates from all contacts that occur over the duration of a partnership. However, this approach has the disadvantage of implicitly assuming that the number of acts per partnership is static, and is independent of the number of concurrent partnerships. In other words, if a person has 10 times as many simultaneously sexual partners, conventional model approaches may implicitly assume that each person has 10 times as many sexual encounters in a particular time period. However, data suggest that people with many simultaneous partners do not have a proportionate increase in frequency of sexual encounters per unit time, and therefore this implicit assumption may exaggerate the impact of concurrency on transmission risk .
We chose not to represent the composite risk that accumulates from all events over the duration of a partnership, but instead represented the underlying determinates of this composite risk: frequency of sexual contacts, concurrency of sexual partnerships, and duration of sexual partnerships. Each characteristic is a separate ‘dial’ that can be increased or decreased without directly affecting the other ‘dials’; however, it is important to note that there may be important correlations between these inputs so they should only be varied independently with great caution, such as in exploratory analyses. For example, increasing the frequency of sexual encounters without changing the concurrency or duration of partnerships will increase the transmission risk per partnership without altering the partner change rate. As another example, increasing a person's concurrency of sexual relationships without increasing their overall sexual frequency or duration of those relationships will reduce the transmission risk per partnership because of fewer sexual acts per partnership. However, the transmission risk per unit time may increase because the number of simultaneous partnerships increases. Decomposing risk per partnership into its constituent characteristics is particularly important for modeling the sub-Saharan HIV pandemic because nonmonogamous contacts (e.g., CSWs) may have long durations . The model's calibration is described in Appendix B – Model Calibration, http://links.lww.com/QAD/A429.
Comparing results with and without transmission
As a validity test of our ability to discriminate how results varied with versus without transmission, we ran the transmission model in a mode that would approximate the fixed-sized, noninteractive cohort of the progression model by turning off births and transmission. We found that, by only incorporating results from known HIV-infected patients, the transmission model results converged with corresponding results from the progression module alone.
Other sensitivity analyses
We varied input assumptions that were likely to impact which strategies lied on the efficient frontier, based on our previous published work . In particular, we varied assumptions regarding the cost of second-line ART relative to first-line ART, the expense of the viral load test, and both of these together. We also varied assumptions regarding the overall transmissibility of HIV (likelihood of transmission per sexual encounter of specified risk).
Our model suggested different efficient frontiers regardless of whether our simulation considered the implications of different monitoring strategies on disease progression alone, or whether it also considered the implications of different monitoring strategies on transmissibility.
Considering progression alone
When our simulation evaluated different monitoring strategies considering only its impact on progression, efficient frontiers increased stepwise, beginning at strategies that did not involve any laboratory monitoring, increasing to strategies that involved CD4+ cell testing with follow-up virological testing upon positive results, and further increasing to strategies that involved progressively more frequent routine virological testing with progressively lower switching thresholds (Fig. 2a). The ICER of earlier ART (CD4+ of 350 cells/μl rather than 200 cells/μl, $857–2388 per QALY, considering progression only; Appendix C – Additional Results, http://links.lww.com/QAD/A430) was more favorable than any of the routine virological testing strategies, which started at $6557 per QALY, and rose to as high as $25 370 per QALY (Table 2a–d). In other words, it would only maximize health to obtain routine virological testing in situations in which coverage with first-line ART is already at goal, and in situations in which ART can be started widely at CD4+ cell count of 350 cells/μl rather than 200 cells/μl.
Considering transmission as well as progression alone
When our simulation evaluated different monitoring strategies considering transmission as well as progression, the value of routine virological monitoring increased greatly, with ICERs cut by approximately one-third (decreasing from $6557/QALY to $4723/QALY; Table 2a–d). However, the value of routine virological monitoring strategies remained less favorable than that of earlier ART initiation (CD4+ cell count of 350 versus 200 cells/μl, $919 per QALY; CD4+ cell count of 500 versus 350 cells/μl, $1522 per QALY; considering transmission as well as progression; Table 3). Considering transmission in addition to progression increased the value of routine virological testing by taking into account the reduced infectiousness of patients who are effectively managed (lowering the average number of infections per infected person per year assuming immediate ART from 0.0679 to 0.0631), and decreased the number of new infections over 20 years by 9.9% from 1 268 591 to 1 142 956 (Table 3). Efficient frontiers (Fig. 2b) continued to increase stepwise, beginning at strategies that did not involve any laboratory monitoring, increasing to strategies that involved CD4+ cell testing with follow-up virological testing upon positive results, and further increasing to strategies that involved progressively more frequent routine virological testing with progressively lower switching thresholds.
As a result, even when considering transmission as well as progression, and considering the resulting decrease in infectiousness that more aggressive monitoring could yield, frequent virological routine monitoring would only be preferred in those settings in which coverage of first-line regimens is already maximized, and in settings in which coverage is already initiated at earlier CD4+ cell counts (Table 3).
Considering transmission along with progression, when we assumed that second-line ART was no more expensive than first-line ART (and, therefore, reducing the penalization of testing regimens that had frequent false positives and would consequently result in far more expensive ART), we found that routine viral load testing became far more attractive. The strategies on the frontier not involving routine viral load testing became sparse, and routine viral load testing offered favorable value (as low as $1331/QALY), now in the range of earlier ART initiation (data not shown).
When we assumed that viral load testing was extremely cheap, as might happen with advances in technology and point-of-care testing, the ICER of routine virological testing decreased but not nearly as much as when we assumed that second-line ART was no more expensive than first-line ART. Routine viral load testing every 6 months now had an ICER of $2386/QALY if the switching threshold was 10 000 copies/ml, rising to $5917/QALY if the switching threshold was 1000 copies/ml (data not shown). However, these values remained less favorable than earlier ART initiation.
When we assumed that second-line ART was no more expensive than first-line ART, and when we also assumed that viral load testing was extremely cheap, we found that routine virologic testing was clearly the preferred choice. Indeed, there was only one strategy on the efficient frontier, and this strategy specified the most frequent viral load testing (every 6 months) and the least conservative viral load switching threshold (1000 copies/ml).
Our analyses underscore the difficult tradeoffs that are necessary to maximize health when there are simultaneous resource-constrained decisions. Notably, routine virological testing only achieved favorable expectancy after coverage with first-line ART was already at goal levels, and after initiation was achieved at a CD4+ cell count of 500 cells/μl. This remained true regardless of whether we considered progression alone, or whether we considered transmission as well as progression.
However, considering transmission in addition to progression substantially improved the value of routine viral load testing. The approximately one-third drop in ICER of routine viral load testing once transmission is considered is not entirely surprising, as more vigilant monitoring and switching would be expected to reduce community viral load and to reduce transmission. Indeed, our results underscore the important relationship between HIV treatment and progression. Our results suggest that HIV monitoring strategies that treat HIV more aggressively have the potential to prevent 9.9% of new HIV infections if ART is initiated instantly. However, this health improvement arose at a sufficiently high cost per benefit that our results suggest that it would only maximize health benefit to use routine virological testing in situations in which ART is already initiated at 500 cells/μl and goal levels of coverage are already met. The only exception to this inference occurred in sensitivity analyses when we assumed that first-line and second-line ART regimens had the same price, because most of the costs of viral load testing originate from downstream costs from more expensive ART regimens, rather than from the costs of the tests themselves.
It is likely that our results are generalizable to East African regions with similar HIV epidemics and similar cost structures, particularly as calibration of the model required varying only two parameters, lessening the danger that it was ‘overfit’ to one particular country. Additionally, most input data were from Kenya, whereas most calibration data were from Uganda, further increasing the likelihood that the model is generalizable across East African regions. However, we make this observation with considerable caution, as there is extensive variability between pandemic epidemiology and cultural factors among different regions in East African countries.
Our results are consistent with recent analyses of three mathematical models by Keebler et al.  examining the cost–effectiveness of different laboratory monitoring strategies, one of which is the progression module of the simulation described in the current report except applied to different sub-Saharan countries and using country-specific costs estimated by WHO rather than AMPATH. Although these three models were developed independently with different assumptions and structures, they consistently found that routine viral load testing delivered lesser incremental health benefits than using those same resources to increase ART coverage and initiating ART at earlier thresholds. We found that our current estimates for the ICER of earlier ART were consistent with estimates by Eaton et al.  of the ICER of earlier ART (ranging from approximately $200/DALY averted to $3800/DALY averted). Our progression-only estimate of the ICER of earlier ART (between $857/QALY and $2388/QALY) was more favorable than our own prior estimate of the ICER of earlier ART initiation ($2600/QALY)  because the prior analysis made different assumptions regarding available ART regimens and used a different initial CD4+ cell distribution for the cohort.
Our results have substantial limitations. All results reflect the assumption that two ART regimens were available, and preferred monitoring strategies clearly vary with the number of available regimens. On the one hand with only one regimen available, no virological testing is favored, as there are no decisions to be made as a result of that testing. On the other hand, when more regimens are available, more virological testing may be preferred because there is less of a danger of ‘burning through’ available regimens. Our model does not consider transmission through MSM or IDU mechanisms. While these mechanisms are likely to account for a fairly small proportion of infections in East Africa overall, same sex is very stigmatized and therefore difficult to estimate its true importance, and IDU is concentrated in the eastern Coastal region, where our results may be less generalizable.
In summary, our results suggest an important role of routine virological testing in improving life expectancy and quality-adjusted life expectancy of individuals with HIV, as well as a potential important role of routine virological testing in preventing new HIV infections in East Africa if the cost of second-line ART regimens were to decrease. However, if all HIV patients are not yet treated with ART starting at 500 cells/μl, resources would buy more health for the population if they are spent on heading toward the ART goal rather than being spent on increasing the frequency of routine virological testing.
This work is supported by National Institute of Allergy and Infectious Disease Award UO1AI069911–01 (IeDEA East Africa) and by the National Institute of Alcohol Abuse and Alcoholism award R01 AA017385, US National Institutes of Health. This work has utilized computing resources at the High Performance Computing Facility of the Center for Health Informatics and Bioinformatics at New York University Langone Medical Center.
Conflicts of interest
1. 2010; World Health OrganizationAntiretroviral therapy for HIV infection in adults and adolescents: recommendations for a public health approach
2. Braithwaite RS, Nucifora KA, Yiannoutsos CT, Musick B, Kimaiyo S, Diero L, et al. Alternative antiretroviral monitoring strategies for HIV-infected patients in East Africa: opportunities to save more lives?
. J Int AIDS Soc
3. Keebler D, Revill P, Braithwaite S, Phillips A, Blaser N, Borquez A, et al. Cost-effectiveness of different strategies to monitor adults on antiretroviral treatment: a combined analysis of three mathematical models
. Lancet Glob Health
4. Markowitz H. Portfolio selection
. J Finance
5. Bajunirwe F, Arts EJ, Tisch DJ, King CH, Debanne SM, Sethi AK. Adherence and treatment response among HIV-1-infected adults receiving antiretroviral therapy in a rural government hospital in southwestern Uganda
. J Int Assoc Physicians AIDS Care (Chic)
6. Johnson VA, Brun-Vezinet F, Clotet B, Gunthard HF, Kuritzkes DR, Pillay D, et al. Update of the drug resistance mutations in HIV-1: December 2009
. Top HIV Med
7. Braithwaite RS, Shechter S, Chang CC, Schaefer A, Roberts MS. Estimating the rate of accumulating drug resistance mutations in the HIV genome
. Value Health
8. Braithwaite RS, Kozal MJ, Chang CC, Roberts MS, Fultz SL, Goetz MB, et al. Adherence, virological and immunological outcomes for HIV-infected veterans starting combination antiretroviral therapies
9. Freedberg KA, Scharfstein JA, Seage GR 3rd, Losina E, Weinstein MC, Craven DE, et al. The cost-effectiveness of preventing AIDS-related opportunistic infections
10. United Nations, Department of Economic and Social Affairs, Population Division. World Population Prospects: the 2008 Revision. 2010.
11. WHO. Epidemiological Fact Sheet on HIV
and AIDS, Uganda. The UNAIDS/WHO Working Group on Global HIV
/AIDS and STI Surveillance. 2008.
12. Boily MC, Baggaley RF, Wang L, Masse B, White RG, Hayes RJ, et al. Heterosexual risk of HIV-1 infection per sexual act: systematic review and meta-analysis of observational studies
. Lancet Infect Dis
13. Attia S, Egger M, Muller M, Zwahlen M, Low N. Sexual transmission of HIV according to viral load and antiretroviral therapy: systematic review and meta-analysis
14. WHO. World Health Organization. Global prevalence and incidence of selected curable sexually transmitted infections. 2001.
15. The Voluntary HIV
-1 Counseling and Testing Efficacy Study GroupEfficacy of voluntary HIV-1 counselling and testing in individuals and couples in Kenya, Tanzania, and Trinidad: a randomised trial
16. Irungu TK, Varkey P, Cha S, Patterson JM. HIV voluntary counselling and testing in Nakuru, Kenya: findings from a community survey
. HIV Med
17. Mossdorf E, Stoeckle M, Vincenz A, Mwaigomole EG, Chiweka E, Kibatala P, et al. Impact of a national HIV voluntary counselling and testing (VCT) campaign on VCT in a rural hospital in Tanzania
. Trop Med Int Health
18. Rosen S, Fox MP. Retention in HIV care between testing and treatment in sub-Saharan Africa: a systematic review
. PLoS Med
19. Kapiga SH, Sam NE, Mlay J, Aboud S, Ballard RC, Shao JF, et al. The epidemiology of HIV-1 infection in northern Tanzania: results from a community-based study
. AIDS Care
20. Kapiga SH, Sam NE, Shao JF, Renjifo B, Masenga EJ, Kiwelu IE, et al. HIV-1 epidemic among female bar and hotel workers in northern Tanzania: risk factors and opportunities for prevention
. J Acquir Immune Defic Syndr
21. Ao TT, Sam NE, Masenga EJ, Seage GR 3rd, Kapiga SH. Human immunodeficiency virus type 1 among bar and hotel workers in northern Tanzania: the role of alcohol, sexual behavior, and herpes simplex virus type 2
. Sex Transm Dis
22. Mbizvo MT, Machekano R, McFarland W, Ray S, Bassett M, Latif A, et al. HIV seroincidence and correlates of seroconversion in a cohort of male factory workers in Harare, Zimbabwe
23. Quigley M, Munguti K, Grosskurth H, Todd J, Mosha F, Senkoro K, et al. Sexual behaviour patterns and other risk factors for HIV infection in rural Tanzania: a case-control study
24. Mishra V. Concurrent sexual partnerships and HIV
infection: evidence from national population-based surveys. DHS Working Papers; 2009; Calverton, MD.
25. Volle J, Letsatsi T, Tan A. A baseline survey of multiple and concurrent sexual partnerships among Basotho men in Lesotho
. Washington DC:C-Change/AED; 2009.
26. Mmbaga EJ, Leyna GH, Mnyika KS, Klepp KI. Sexually transmitted infections knowledge and its impact in the practice of risky sexual behaviours and HIV serostatus: results from rural Kilimanjaro, Tanzania
. Sex Transm Infect
27. Vandepitte J, Lyerla R, Dallabetta G, Crabbe F, Alary M, Buve A. Estimates of the number of female sex workers in different regions of the world
. Sex Transm Infect
2006; 82 (Suppl 3):iii18–iii25.
28. Westercamp N, Mattson CL, Madonia M, Moses S, Agot K, Ndinya-Achola JO, et al. Determinants of consistent condom use vary by partner type among young men in Kisumu, Kenya: a multilevel data analysis
. AIDS Behav
29. 2005; Marks G, Crepaz N, Senterfitt JW, Janssen RS. Meta-analysis of high-risk sexual behavior in persons aware and unaware they are infected with HIV in the United States: implications for HIV prevention programs
30. Sawers L, Isaac AG, Stillwaggon E. HIV and concurrent sexual partnerships: modelling the role of coital dilution
. J Int AIDS Soc
31. Voeten HA, Egesah OB, Ondiege MY, Varkevisser CM, Habbema JD. Clients of female sex workers in Nyanza province, Kenya: a core group in STD/HIV transmission
. Sex Transm Dis
32. Eaton JW, Menzies NA, Stover J, Cambiano V, Chindelevitch L, Cori A, et al. Health benefits, costs, and cost-effectiveness of earlier eligibility for adult antiretroviral therapy and expanded treatment coverage: a combined analysis of 12 mathematical models
. Lancet Glob Health