Objectives This study assesses whether the observed declines in HIV prevalence since the beginning of the ‘Avahan’ India HIV/AIDS prevention initiative are consistent with self-reported increases in condom use by female sex workers (FSWs) in two districts of southern India, and provides estimates of the fraction of new infections averted among FSWs and clients due to increases in condom use in commercial sex after 2004.
Methods A deterministic compartmental model of HIV/sexually transmitted infection (STI) transmission incorporating heterogeneous sexual behaviour was developed, parameterised and fitted using data from two districts in Karnataka, India. Three hypotheses of condom use among FSWs were tested: (H0), that condom use increased in line with reported FSW survey data prior to the Avahan initiative but remained constant afterwards; (H1) that condom use increased following the Avahan initiative, in accordance with survey data; (H2) that condom use increased according to estimates derived from condom distribution data. The proportion of fits to HIV/STI prevalence data was examined to determine which hypothesis was most consistent.
Results For Mysore 0/36/82.7 fits were identified per million parameter sets explored under hypothesis H0/H1/H2, respectively, while for Belgaum 9.7/8.3/0 fits were identified. The HIV epidemics in Belgaum and Mysore are both declining. In Mysore, increases in condom use during commercial sex between 2004 and 2009 may have averted 31.2% to 47.4% of new HIV infections in FSWs, while in Belgaum it may have averted 24.8% to 43.2%, if there was an increase in condom use.
Discussion Increased condom use following the Avahan intervention is likely to have played a role in curbing the HIV epidemic in Mysore. In Belgaum, given the limitations in available data, this method cannot be used alone to decide if there has been an increase in condom use.
- mathematical model
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
Karnataka, in southern India, has one of the highest prevalences of HIV among Indian states,1 and is among the states included in the ‘Avahan’ India HIV/AIDS initiative established by the Bill & Melinda Gates Foundation in 2004. Working through state-level partners and local non-governmental organisations (NGOs), the initiative is a multifaceted HIV preventive intervention promoting prevention strategies among high-risk groups, including distribution of condoms to and promotion of condom use among female sex workers (FSWs).2–4
Assessing the impact of the Avahan intervention is crucial, not only to ensure that Avahan achieves its goals of reducing HIV transmission in India, but also to inform future large-scale interventions.5 However, the dynamics of an HIV epidemic are complicated: declines in HIV prevalence can occur without any intervention effects, due to natural infection dynamics and, conversely, increases can occur even in the presence of an effective intervention.6 7 In order to separate out the impact of the Avahan initiative from natural HIV transmission dynamics, and because of the difficulty and huge expense of implementing community randomised controlled trials, a tailor-made transmission dynamics model will be used as one important component of the impact assessment, the framework for which is described in Boily et al.8
The final impact and cost effectiveness analysis, planned for 2011, will estimate the impact of the Avahan initiative on the HIV and sexually transmitted infection (STI) epidemics among high-risk groups and the broader population across more than 30 intervention sites.5 8 This paper presents the results of an interim analysis focusing on FSWs and their clients in two districts where sufficient data are currently available, and which represent different HIV/STI epidemic trends and behavioural profiles: Belgaum Urban, a higher prevalence district in which Avahan was not the first but is now the only intervention; and Mysore Urban, a lower prevalence district in which Avahan was the first and remains the only intervention.
The objective of this analysis is to determine whether the observed changes in HIV prevalence among FSWs could be solely due to natural disease dynamics rather than to a change in condom use following the intervention. We address this question by testing whether the null hypothesis (H0) of stable (ie, no increase in) condom use following the initiation of the Avahan intervention in 2004 is less likely than the alternative hypotheses (H1) that condom use actually increased by as much as reported in FSW survey data, or (H2) that there was a lower level of preintervention condom use and a sharper postintervention increase than indicated by FSW survey data, as suggested by records of condom availability.9
In order to model the HIV/STI epidemics in urban Belgaum and Mysore districts, an in-depth mathematical modelling analysis was undertaken with the key elements described below.
Methodological framework of evaluation
The structure of the evaluation follows the methodological framework designed and described in Boily et al.8 In short, a mathematical model was used within a Bayesian framework, where available data were first used to specify what was known about each model parameter by defining a plausible range of values (the ‘prior distribution’). This distribution was then randomly sampled repeatedly, to test the model with a large number of parameter combinations and identify those parameter sets (the ‘posterior distribution’) that agreed with the empirical HIV/STI prevalence data shown in table 1.
One of the advantages of this approach is that the percentage of parameter sets producing a model fit to the prevalence data gives a measure of how compatible the prior parameter sets were with the observed prevalence data for the specified model. This approach was used here to test the three hypotheses described above concerning possible trends in condom use before and after the Avahan intervention. Finally, the posterior distribution was used to simulate a control group without the intervention, which was then compared to the original scenario to produce estimates of intervention impact.
Model structure and parameterisation
A purpose-built deterministic compartmental model of HIV, herpes simplex virus (HSV)-2 and syphilis transmission was constructed and parameterised to incorporate key aspects and heterogeneities of the complex FSW and client structure specific to each setting.10 In summary, the model consisted of an open population, growing at the rate described by census data,11 stratified into ‘high-risk’ (FSWs and their clients) and ‘low-risk’ (the remaining population) groups. Based on work by Vickerman et al4 (see page 33), it was assumed in this model that there is no transmission between low-risk individuals. The high-risk groups (FSWs and their clients) were further stratified by variables showing strong associations with HIV prevalence and high-risk behaviour in analyses of data from serial cross-sectional surveys termed integrated behavioural and biological assessments (IBBAs).10 12
The model simulates the transmission of HIV/STIs between FSWs and their clients through commercial and longer-term non-commercial partnerships, and the bridging infections from clients to low-risk women through their non-commercial partnerships, as shown in figure 1. The force of infection depends on the disease-specific infectivity, type and duration of partnership, frequency of sex acts for the type of partnership, condom use and STI coinfection status of both partners. The three parameters influencing HIV prevalence the most in a sensitivity analysis of a representative sample of the unrestricted model runs were the number of sex acts per partnership between FSWs and occasional clients, the male to female HIV transmission per sex act and the RR for increased transmissibility of primary HIV. Further details of the model and equations are provided in the supplementary material.
HIV was modelled with an initial short acute phase of high infectivity, followed by a long low-infectivity phase and a pre-AIDS phase of increased infectivity.13 It was assumed that those with AIDS are chronically ill, and cease being sexually active. HSV-2 coinfection14 and syphilis15 were also modelled dynamically, with cofactors representing facilitation of HIV and HSV-2 acquisition and transmission.
The setting-specific behaviour-related, demographic-related and intervention-related prior model parameter distributions were derived from detailed serial cross-sectional behavioural and biological data (from IBBAs) collected from FSWs and their clients, and general population surveys (GPS) carried out as part of the Avahan monitoring programme, as well as from complementary data from other sources. The non-setting specific biological parameter ranges were based on reviews of relevant literature (tables 2 and 3). Prior distributions were assumed to be uniform across each parameter range; this was a ‘conservative’ choice to reflect the lack of precise knowledge and possibility of bias in the data used to estimate the parameters.
The parameterisation of condom use before and after the intervention was data driven and specific to each hypothesis tested, resulting in three different prior condom use distributions. The prior distribution of condom use for hypothesis H1 was based on estimates (with CIs) of site-specific condom use from 2001 onwards given by the analysis of Lowndes et al.75 That analysis used data from the IBBA round 2 surveys concerning time since starting to use condoms consistently with occasional clients to retrospectively ‘reconstruct’ the fraction of FSWs consistently using condoms each year from 2001 to the year of the survey. The null hypothesis H0 assumed the same increase in condom use prior to 2004 as hypothesis H1 but no further increase after the start of the Avahan intervention in 2004. A further hypothesis (H2) was based on HIV prevention programme data and other sources regarding the availability of condoms for FSWs9 from 2004 to 2008 across Karnataka. This study suggested that the proportion of consistent condom users among FSWs in 2004 was much lower than suggested in hypothesis H1 due to low availability, but that it rose steeply during 2006–2008 due to the Avahan intervention. The three different condom prior distributions are shown in figure 2, with the round 1 and 2 FSW IBBA point estimates shown for comparison.
Parametric uncertainty analysis, model fitting and validation
This step involved identifying combinations of model parameters that agreed with (‘fit’) the HIV/STI prevalence data from each site for each condom use hypothesis tested. Latin hypercube sampling (LHS) was used to sample several million parameter sets from all ‘prior’ parameter ranges for each hypothesis and district, with the sampling done so that the same set of parameters, excluding the intervention parameters being tested, were used for the different hypotheses in each district. The model was run with each sampled parameter set to predict HIV/STI prevalences. A parameter set was accepted as a fit if the model predictions simultaneously lay within the 95% CIs of the HIV, HSV-2 and syphilis prevalence data of the relevant FSW and client IBBAs, as described below.
The HIV/STI prevalence data from FSWs and clients for each district were used in two stages: (1) as validation and (2) for hypothesis testing. In stage 1, the model was first fitted to just the round 1 FSW and client IBBA HIV, HSV-2 and syphilis prevalence data (see table 1), and the predictive value of these model fits was assessed against the HIV prevalence data not used directly at this stage of the fitting process, namely FSW HIV prevalence from the second round FSW IBBA in Mysore and Belgaum (table 1) and third round in Mysore (table 1), and FSW HIV prevalence by typology (street, brothel or home based) in rounds 1 and 2 (see supplementary material). In stage 2, the second round of IBBA HIV prevalence data were included in the fitting procedure along with the first round of FSW and client IBBA HIV/STI prevalence data (table 1) to produce a reduced set of model fits for hypothesis testing.
Testing hypotheses of condom use
The most likely condom use hypothesis was determined based on the relative proportions of fits to the first round FSW and client IBBA HIV/STI prevalence data and the second round FSW HIV prevalence data. For each hypothesis, the fraction of model fits from all parameter sets was compared to determine which hypothesis was more likely. This method, which is related to the version of Approximate Bayesian Computing as described by Weiss and von Haeseler,76 is also similar to the Bayesian methods used by Vickerman et al for determining the likelihood of different injecting drug user risk behaviour scenarios in Pakistan,77 and Hallett et al for examining behaviour change in Zimbabwe.78
The posterior distribution of the most likely condom use hypothesis was used to simulate HIV prevalence and incidence over time following the start of the intervention, and in a matched control population. In the latter, each model fit was run again with condom use kept constant at 2004 levels (ie, at the start of the Avahan initiative). Comparing the impact and control scenarios, estimates of the impact of the increase in condom use following the start of the intervention among high-risk groups were derived for a 5-year period (2004–2009). These estimates derived from the posterior distribution include credibility intervals (CrI) generated by using all the parameter combinations that fit the prevalence data to reflect the uncertainty in parameter assumptions, since more than one set of parameters may produce an equally good fit to epidemiological trends.
The proportion of fits to each hypothesis is shown in table 4. In Mysore, the hypothesis H2, derived from analysing availability of condoms, gives the largest proportion of simulations which are consistent with prevalence data, and scenarios H2 and H1 of increasing condom use following the introduction of Avahan are more likely than the null hypothesis H0 of no condom increase after the start of Avahan, for which no model fits were obtained. In Belgaum, however, the null hypothesis H0 of no increase use in condom use following Avahan and the alternative hypothesis H1, based on the historical condom use reconstruction, cannot be distinguished by this method, although both are much more likely than the hypothesis H2. An alternative method to comparing the proportion of fits would be a likelihood-based approach, calculating the likelihood for each run; however this method gave comparable results (data not shown).
Crossvalidation of results
Table 4 shows that fits to only the round 1 FSW and client IBBA prevalence data (table 4) also generally fit the round 2 FSW IBBA HIV prevalence data well for all hypotheses (90% or more of posterior runs also fit FSW HIV round 2 prevalence). This then validates the model with respect to time trends in HIV (see also supplementary material). A third round of IBBA data was available for Mysore only, and the model HIV prevalence projections were slightly higher (12.1% to 20.6% for H2; 12.4% to 19.6% for H1) than the third round IBBA data (8.1% to 14.1%). This is consistent with the fact that any further increases subsequent to the IBBA round 2 in condom use were not incorporated in the model.
Lastly, validation of the model by comparing its projections with FSW HIV prevalence by typology generally showed good agreement, although there was some residual heterogeneity, as might be expected (see supplementary material).
Trends in HIV/STI prevalence and incidence
The predicted prevalences over time for the model fits from the favoured hypotheses (H0 and H1 for Belgaum and H2 for Mysore) suggest that the HIV epidemics among FSWs and clients are declining in both districts (figure 3A–F). In Belgaum, it is likely that this decline began before Avahan started, around the year 2000, while in Mysore the decline probably started shortly after Avahan, in 2004–2005. Modelled incidence and the prevalence for hypothesis H1 for Mysore show similar trends (see supplementary material). More details of the general characteristics of the model runs and fits are included in the supplementary material.
In Mysore, the median relative decline in HIV prevalence in 2009 compared to 2004 (‘time decline’) under H2 was 36.2% (95% CI 21% to 51%) in FSWs (median 25.6%, 95% CI 7% to 42% in clients), whereas the median relative decline in 2009 prevalence compared to a simulated scenario without any increase in condom use since 2004 (‘intervention decline’) was 44.7% (95% CI 35% to 52%) in FSWs (37.9%, 95% CI 30% to 45% in clients). In Belgaum, under H1, the median time decline was 47.6% (95% CI 39% to 58%) in FSWs (median 55.7%, 95% CI 48% to 66% in clients) while the median intervention decline was actually lower at 23.1% (95% CI 17% to 30%) in FSWs (28.4%, 95% CI 22% to 36% in clients). Thus, the time decline in HIV prevalence slightly underestimates the intervention decline in Mysore for hypothesis H2 and overestimates it in Belgaum for hypothesis H1.
Intervention impact among high-risk groups
We estimated the impact of the increase in condom use following the start of Avahan using the most likely hypothesis as the simulated intervention group/population, and using the null hypothesis H0 as the simulated matched control group/population. In Mysore, the increase in condom use between 2004 and 2009 may have prevented 31.3% to 47.4% of new HIV infections in FSWs and 32.7% to 47.2% in clients under hypothesis H2 (figure 4 shows this in comparison to the impact projections for hypothesis H1). Impact increases over time, with 58.7% to 82.2% of new HIV infections prevented among FSWs in 2008 and 64.1% to 85.1% among clients under hypothesis H2 (41.2% to 72.2% and 47.1% to 77.6% under H1).
In Belgaum, hypothesis H0, where condom use remains constant from 2004 onwards, suggests no impact. In comparison, the H1 scenario projects that 24.8% to 43.2% of new FSW HIV infections and 30.9% to 53.1% of new client infections were prevented in Belgaum by the increase in condom use between 2004 and 2009 (figure 4). As for Mysore, impact increases over time in the H1 scenario, with 32.7% to 59.7% of new infections in FSWs and 42.0% to 72.6% in clients prevented in 2008.
With sufficient data to inform model parameters and validate model predictions, the use of mathematical models within a Bayesian framework enables the testing of hypotheses about the value or range of more uncertain parameters,77 78 such as the evolution of condom use over time. In this analysis we have tested whether condom use is likely to have increased following the implementation of the Avahan intervention. In theory this could have been done more directly if unbiased baseline data on condom use had been collected; however, in practice it was not possible because such data could only be collected once the intervention was sufficiently well established to map and sample FSWs.
The success rates obtained in the hypothesis-testing analysis support the conclusion that in Mysore condom use during commercial sex increased substantially after the beginning of the intervention. This increase most likely occurred shortly after the beginning of the intervention, as suggested by data on availability of condoms. Although there was no significant decrease in the FSW IBBA prevalence data between rounds 1 and 2, the model suggests that, in the absence of an increase in condom use after 2004, prevalence would have increased, and so to obtain the required trend a significant increase in condom use was required. In Belgaum, the model results show that decreases in prevalence alone from the IBBA surveys are insufficient to distinguish between the hypotheses of increased condom use as reported by FSWs and constant condom use following the beginning of the intervention. Without information from other sources, it is not possible to definitely conclude that there has been an increase in condom use in commercial sex in Belgaum. The results also suggest that the level of condom use was higher in Belgaum than in Mysore prior to 2004. In Belgaum, a smaller-scale intervention was present for some years prior to 2004,79 before it was scaled up by Avahan. Condom use by FSWs with clients may therefore have been higher at the start of the intervention than in Mysore, where there was no prior intervention.
The model results indicate that HIV incidence and prevalence have declined over time. The decline occurred earlier in Belgaum than in Mysore because condom use became widespread earlier, but perhaps also because the epidemic matured earlier. In Mysore, the decline roughly coincides with the beginning of the intervention. This pattern of decrease in prevalence agrees with that seen by Vickerman et al,4 and by Boily et al80 examining ANC data. While this study only considers high-risk individuals, the decline observed among FSWs and clients should translate to a slower decline in the general population in both districts, as shown in the analysis by Vickerman et al.4
The relative decline in HIV prevalence in Mysore in 2009 compared to the simulated control group (‘intervention decline’) is higher than the decline over time since 2004 (‘time decline’). This is because the epidemic was still growing in 2004. In Belgaum, under the hypothesis of increasing condom use after Avahan, the intervention decline is smaller than the time decline, as the epidemic there was already declining prior to Avahan. This highlights that time trends alone should not be used to evaluate an intervention, and that modelling results are important to complement prevalence data in carrying out an evaluation.6 8 80
Strengths and limitations
The model described in this paper was specifically tailored to reflect the behavioural and biological heterogeneities present in these settings, with dynamical modelling of HSV-2 and syphilis in addition to HIV in each site. It was parameterised using data from multiple surveys among FSWs, clients and general population, and captured key epidemiological trends in HIV, HSV-2 and syphilis by risk groups. This model also compared well to an independent model by Vickerman et al,4 and it was also validated with respect to time trends in HIV and FSW HIV prevalences by typology.
The framework described in Boily et al8 was used to produce point estimates and credibility intervals of intervention impact by the creation of control groups that assumed no increase in condom use. This control group gives the impact due to the increase in condom use since the beginning of the intervention, and, in situations where condom use would have increased naturally even without an intervention to drive it, will overestimate intervention impact. Conversely, in the absence of a likely increase in condom use following Avahan, as may have been the case in Belgaum, this choice of control group results in the model projecting no HIV infections averted, while it is possible that there would still have been a decline in condom use in this situation without an intervention to sustain it. Indeed, analysis of a slow, modest decrease in consistent condom use of 10% over 5 years suggests that even such a small decline in condom use leads to a large impact (17.3% to 29.7% additional new FSW HIV infections between 2004–2009 due to the decrease in condom use), and thus the contribution of the intervention in sustaining condom use is important, although it is difficult to attribute directly.
The precision of the estimates was limited by data availability. In both districts, the fitting success rate was modest because conservative (wide) prior ranges were used for most parameters to capture their uncertainty, and because the model had to fit seven separate prevalence outcomes. If more data were available, then it would be possible to use this additional information to reduce the width of these priors and so possibly identify a higher proportion of model fits. This would allow a stronger comparison between the condom use scenarios, and would better reflect uncertainty in the impact projections. HIV prevalence data for earlier in the epidemic would also enable us to better constrain the trajectory of the early epidemic and so rule out unlikely epidemics with high early prevalence. However, using a constraint such as FSW prevalence being <10% in 1985 did not change the impact estimates significantly (data not shown).
The picture is further complicated in Belgaum because the first IBBA survey was carried out 16 months after the start of the Avahan initiative, and so behavioural and epidemiological data reflected a situation in which Avahan was already present and where consistent condom use rates reported by FSW in the 2005 IBBA were already very high, at 91%. In contrast, the survey in Mysore is likely to be more comparable to the pre-Avahan baseline as it took place 8 months after the start of Avahan. It is difficult to quantify the impact of a change in behaviour in situations where there was a pre-existing intervention and a long delay between the start of the intervention and the first IBBA survey, such as in Belgaum. Impact estimation is more uncertain in such settings because it is hard to define and simulate a control group reflecting how condom use would have decreased without the intervention, for which no data are available. An additional round of behavioural and HIV prevalence data for FSWs in Belgaum could help to further inform impact estimates by giving extra information about how the epidemic there is changing.
In Belgaum, the low number of fits from the model to the prevalence data may be largely explained by the difficulty of fitting the relatively high HIV prevalence in clients. This analysis modelled the transmission of HIV/STIs among urban clients, since the IBBA was carried out in urban areas where the intervention started. However the size of the FSW population is higher in rural areas, and as city clients may also go to smaller towns and villages to visit FSWs, the urban client HIV prevalence may be influenced by the large rural FSW population. Future data collection strategies for programme evaluation should be adapted to reflect the reality of how the epidemic is driven and the nature of the intervention.
Dynamical modelling, used within a Bayesian framework, is useful for testing hypotheses and providing informative and less subjective impact estimates on which to base decisions. This method of comparing proportions of fits of different scenarios is relatively straightforward to implement, as shown here and in another recent analysis,77 and helps to further augment survey findings of behaviours that are prone to reporting biases, such as condom use.
This analysis suggests that the HIV epidemics in Belgaum and Mysore are declining. In Belgaum, it is probable that the decline in HIV prevalence started before the beginning of the Avahan intervention, since condom use by FSWs with their clients was already high in this district due to pre-existing interventions. The model, using current prevalence data alone, is unable to distinguish between scenarios where condom use in commercial sex was sustained at existing levels and where it increased after Avahan. If there was a rise in condom use, as suggested by other data9 and equally likely compared to the null hypothesis from this method, this increase would have averted 24.8% to 43.2% of HIV infections in FSWs and 30.9% to 53.1% of infections in clients between 2004 and 2009. If, however, Avahan only sustained the existing level of condom use attained by FSWs prior to Avahan, this is still important for controlling the epidemic, as even a small, gradual decline in condom use can greatly increase the number of new infections. This conclusion on impact, of including both scenarios, is conservative (which is deemed preferable), since the difference in proportion of fits between the hypotheses where condom use does or does not increase after 2004 is small. However, combined with the additional data presented in Bradley et al9 regarding the availability of condoms in Belgaum, it can be argued that it is likely that there has been an increase in condom use since the start of Avahan.
In Mysore, the findings strongly suggest that condom use during commercial sex has increased since the introduction of the Avahan programme, and that this increase in condom use has averted 31.3% to 47.4% of HIV infections in FSWs and 32.7% to 47.2% in clients in Mysore over the first 5 years.
Transmission dynamics models used within a Bayesian framework are useful tools to help evaluate large-scale intervention in absence of a control group. Such a framework enables testing of different hypotheses and provides credibility intervals for estimates of infections averted.
Our modelling analysis suggests there has been a decrease in HIV prevalence in high-risk groups due to an increase in condom use since the beginning of the Avahan intervention in Mysore.
Due to data limitations, our analysis could not say conclusively whether the change in prevalence in Belgaum was due to an increase in condom use following the start of the intervention, or simply the result of the natural dynamics and increase in condom use before 2004.
This research was funded by the Bill & Melinda Gates Foundation. The views expressed herein are those of the authors and do not necessarily reflect the official policy or position of the Bill & Melinda Gates Foundation, the London School of Hygiene & Tropical Medicine, or Imperial College London. Many thanks to Jan Bradley of the CHARME-India project, Bangalore for useful discussions and to the London School of Hygiene and Tropical Medicine for use of the HPC.
Supplementary figures, tables and appendix are published online only at http://sti.bmj.com/content/vol86/issue1
Funding The Bill and Melinda Gates Foundation, Seattle, Washington, USA.
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.