Modeling secondary level of HIV contact tracing: its impact on HIV intervention in Cuba
 YingHen Hsieh^{1},
 YunShih Wang^{2},
 Hector de Arazoza^{3} and
 Rachid Lounes^{4}Email author
DOI: 10.1186/1471233410194
© Hsieh et al; licensee BioMed Central Ltd. 2010
Received: 8 September 2009
Accepted: 1 July 2010
Published: 1 July 2010
Abstract
Background
Universal HIV testing/treatment program has currently been suggested and debated as a useful strategy for elimination of HIV epidemic in Africa, although not without practical issues regarding the costs and feasibility of a fully implemented program.
Methods
A mathematical model is proposed which considers two levels of detection of HIVinfectives through contact tracing of known infectives in addition to detections through other means such as random screening. Simulations based on Cuban contact tracing data were performed to ascertain the potential impact of the different levels of contact tracing.
Results
Simulation studies illustrate that: (1) contact tracing is an important intervention measure which, while less effective than random screening, is perhaps less costly and hence ideal for largescale intervention programs in developing countries with less resources; (2) the secondary level of contact tracing could significantly change the basic disease transmission dynamics, depending on the parameter values; (3) the prevalence of the epidemic at the time of implementation of contact tracing program might be a crucial factor in determining whether the measure will be effective in preventing disease infections and its eventual eradication.
Conclusions
Our results indicate that contact tracing for detection of HIV infectives could be suitably used to remedy inadequacies in a universal HIV testing program when designing timely and effective intervention measures.
Background
A recent modeling study [1] on universal HIV screening followed by immediate antiretroviral treatment (ART) for those tested positive concludes that, assuming that the infectiousness of those treated fell to 1% of their infectiousness before treatment, this strategy could have a major effect on HIV/AIDS epidemic. However, as noted in that article, screening everyone periodically, whether voluntary or not, is not feasible in practice. Moreover, the costs of periodic universal testing program might be prohibitively expensive, especially in the developing countries.
The HIV/AIDS epidemic in Cuba has had significantly low prevalence compared to its neighboring countries in the Caribbean Basin, which has the second highest rate of HIV/AIDS in the world after subSaharan Africa [2]. The Joint United Nations Program on HIV/AIDS (UNAIDS) report on the global HIV/AIDS epidemic for 2007 indicates that the region's HIV epidemics are in general highly prevalent and driven primarily by heterosexual intercourse. In contrast, estimates for Cuba indicate an adult HIV prevalence of 0.1%, out of 6.1 million persons in the age group 1549 [3]. Moreover, Men having sex with men (MSM) continue to be the main high prevalence group in Cuba, with most of the HIVpositive adults being men (77%) and most of the detected HIVpositive men being reported as MSM (85.1%), while most of the HIVpositive women reported having had sex with MSM [4].
The continued low HIV prevalence in Cuba has its root in its distinct social factors as well as in its unique intervention program. The Cuban HIV/AIDS programme included a system to detect HIV cases from several sources. Some of these sources were used at the beginning of the programme, others were introduced later and some have been discontinued. Since 1993, this detection system has focused on 6 major sources: blood donors; persons treated for other sexually transmitted infections; persons admitted to hospital with suspected HIV infection or subject to specific procedures like dialysis; persons volunteering to be tested; persons whose general practitioner has recommended HIV testing; and sexual partner tracing. Other less important sources include testing of all pregnant women and prison inmates. One of the major reasons for the low HIV prevalence in Cuba compared to its neighbors in the Caribbean region is its sexual contact tracing program, also known as the Partner Notification Program. Since 1986, a person testing HIVpositive in Cuba is interviewed by health workers using a nonanonymous structured questionnaire. They are invited to give names and contact details of their sexual partners of the past two years. These partners are then traced and a recommendation for voluntary HIV testing is made with followup testing for up to one year [5].
Modeling of Cuba HIV has been carried out in recent years, some focusing mainly on underreporting [6–8] and others on estimation of HIVinfected population sizes using Generalized Removal Model for Open Population or GERMO [5, 9]. Several studies also focused on modeling the contact tracing program in Cuba includes [8, 10–13], mainly with the Cuba HIV/AIDS epidemic as the background due to the availability of detailed HIV case and contact tracing data.
Sexual transmission is the primary mode of infection in Cuba, in particular that of MSM group [4]. Hence seeking out sexual contact networks through contact tracing is a vital aspect of intervention and control of the HIV epidemic in Cuba. In the previous modeling work on contact tracing, a detection rate by contact tracing is assigned for all detected individuals, thus ignoring the difference between the first level of traced contacts (i.e., those who were traced earliest as members of a newly discovered network of direct sexual contacts of an infective detected through means other than contact tracing) and the succeeding levels of traced contacts (i.e., those in the network who were traced consequently as subsequent contacts of traced contacts in a network). For example, those infectives traced and diagnosed as reported contacts of the first level of traced contacts belong to the second level, and those traced as reported contacts of the second level of traced contacts belong to the third level, and so on. In this work, we will extend our previous model to include a second level of contact tracing, that is, those who were detected as sexual contacts of those who were the earliest detections through contact tracing as contacts of those who were detected through means other than contact tracing.
Fig. 2 gives evidence that the effectiveness of contact tracing decreases sharply as tracing into the tertiary levels is being conducted, most likely due to the gradual exhaustion of the undetected members in the sexual contact network to which the detected infectives belong, and also due to the overlapping between infection chains or sexual networks. Moreover, the information it provides becomes redundant and obsolete as the contact tracing program continues, since largescale contact tracing is a timeconsuming process and requires personnel that must be exclusively dedicates to this task, especially in large population centers. Therefore, we wish to ascertain the true effectiveness of the contact tracing program, by gauging the impact of secondary and perhaps tertiary levels of contact tracing.
Early mathematical models of screening and treatment for HIV include [14–17]. Ample work on modeling sexual partnership network had been carried out in recent years, e.g., [18–25]. However, these models were mostly proposed to construct sexual partnership (concurrent or otherwise) networks.
In the context of the HIV epidemic in Cuba, contact tracing exhibit only partial results of the partnership networks, filtered through nonrandom detection of infected individuals, recall bias of these infected individuals, tracing and testing of the contacts, timing of the tracing and testing, etc. Moreover, the Cuban contact tracing data contains only those tested positive. Therefore, to use the approaches proposed in the above models would lead to an even more complicated (compartmental or individualbased) model which requires additional data on sexual network which is rarely available, with additional assumptions and further uncertainties in the results obtained. Hence, in this work, we will consider a simple compartmental model which enables us to focus on the secondary level of contact tracing and its relative impact.
Methods
The Model
S(t) is the susceptible population.
X(t) is the number of HIVinfected persons that do not know they are infected at time t.
Y_{1}(t) is the number of HIVinfected persons that are known by the health authority to be infected at time t and were detected in a random type search.
Y_{2}(t) is the number of HIVinfected persons that are known by the health authority to be infected at time t and were detected at the first level of contact tracing as past contacts (or sexual partners) of Y_{1}.
Y_{3}(t) is the number of HIVinfected persons that are known by the health authority to be infected at time t and were detected at the second level of contact tracing as past contacts of Y_{2}.
Z(t) is the number of persons with AIDS at time t.
The model parameters are listed below:
Λ: the constant recruitment rate into the susceptible population.
λ: the rate of the new HIVinfected persons are infected by X(t).
k _{1}: the rate at which the unknown HIVinfected persons are detected by the system, independently of other HIVpositive persons (through "random" screening).
k _{2}: the rate at which the unknown HIVinfected persons are detected by the system through contact of Y_{1}.
k _{3}: the rate at which the unknown HIVinfected persons are detected by the system through contact of Y_{2}.
β: the rate at which the undetected HIVpositive persons develop AIDS.
β': the rate at which the detected HIVpositive persons develop AIDS.
μ: mortality rate of the sexually active population.
μ': mortality rate of the population with AIDS.
Following previous modeling work (e.g. [7, 8, 12, 13]), we assume a constant recruitment rate λ of the new infectives, infected by the infectives who did not know that they were infected. The assumption of linear rate of recruitment of new HIVinfected persons infected by X(t), λ X, is based on our background setting of Cuba, where the prevalence is low. We also assume linear detection rates by contact tracing, k_{2} and k_{3}, where a certain fraction of previous contacts of a known infective is successively traced and tested to be HIVpositive. A class of linear and nonlinear functions for detection via contact tracing was considered in [13], in which the nonlinear detection by contact tracing term used in this model is the one found to give best fit for Cuban HIV data. Moreover, estimates for k_{2} were also obtained in [7].
We also assume that the infection by detected infectives is negligible when compared with the infection rate by the unknown infectives, since it was estimated in [12], using Cuban HIV/AIDS data between 19862000, that the infection rate by the known infectives in Cuba is only around 5.79% (SD = 3.55) of that of unknown infectives. We also ignore the contact tracing at the third level and after, since the tertiary levels of detections are small when compared with that of the second level.
Note that System (2.1) is welldefined in the 4dimensional positive orthant of the XY_{1}Y_{2}Y_{3}space. For a theoretical analysis for models of this type near the origin, see [26].
Results
Steady states and threshold conditions
System (3.1) has several dynamic steady states (also known as equilibrium points), each with its corresponding threshold condition for stability.
1. Trivial steady state: E _{0} = (0, 0, 0)
2. Boundary steady states: E _{ B1 }= (0, r _{2},0), E _{ B2 }= (0, 0, r _{3}) and where , and are nonzero numbers between 0 and 1, and their values depend on the initial fractions of the solution trajectories at time t = 0, i.e.,(x(0), y _{2}(0), y _{3}(0)).
(i) E _{ B1 }= (0, r _{2},0), where 0 < r _{2} ≤ 1.
(ii) E _{ B2 }= (0, 0, r _{3}), where 0 < r _{3} ≤ 1.
(iii) , where and .
3. Positive endemic steady state: , .
Note that since we have assumed that the infections by the detected infectives are negligible, the endemic steady state E* is the only endemic steady state where some infectives have not been detected and continue to cause new infections in the community. The threshold condition for each of the steady states is as follows:
(a) If R _{0} < 1, then E _{0} = (0, 0, 0) is locally asymptotically stable for (2.1) where R _{0} = (λ  β + β')/(k _{1} + k _{2}).
(b) If R _{ B1 }< 1, then E _{ B1 }= (0, r _{2},0) is locally asymptotically stable for (2.1) where R _{ B1 }= (λ  β + β')/[k _{1} + k _{2}  (k _{2}  k _{3})r _{2}] and 0 < r _{2} ≤ 1.
(c) If R _{ B2 }< 1, then E _{ B2 }= (0, 0, r _{3}) is locally asymptotically stable for (2.1) where R _{ B2 }= (λ  β + β')/(k _{1} + k _{2}  k _{2} r _{3}) and 0 < r _{3} ≤ 1.
(d) If R _{ B3 }< 1, then is locally asymptotically stable for (2.1) where , and .
The above results can be easily derived by making use of Jacobian matrix and hence is omitted.
Simulations Studies
We will now give some numerical examples to illustrate our results. Recall that all equilibria except the endemic steady state E* listed in the previous section are essentially disease free, since the fraction of undetected HIVinfectives, where individuals in all other infective classes must come from, is 0. We also note that the sizes of threshold parameters, R _{ B1}, R _{ B2}, R _{ B3}, and R*, depend on , and , and subsequently on the initial fractions used for the simulations.
Simulations for model with first and second levels of contract tracing
Parameters values  Reproduction numbers  Limiting steady state  

k _{1}  k _{2}  k _{3}  R _{0}  R _{ B1}  R _{ B2}  R _{ B3}  R* 

0.25  0.25  0.25  1.19  1.19  1.24  1.24  0.67  endemic steady state E* 
0.375  0.25  0.25  0.95  0.95  0.996  0.996  1.004  DFE E _{ B3} 
0.25  0.5  0.25  0.79  0.90  0.88  1.01  0.92  endemic steady state E* 
0.25  0.5  0.25  0.79  0.87  0.81  0.89  0.79  DFE E _{ B3} 
If we increase detection by means other than contact tracing (and in the process increase the number of infectives whose contacts are to be traced) by 50%, all other reproduction numbers are less than 1 with the exception R* which is greater than 1, thus guaranteeing that solutions will become diseasefree (second row of Table 1). In this instance, the solution approaches E _{ B3}, the steady state with nonzero fractions of infectives traced by first and second level contact tracing. However, if instead we double the first level contact tracing rate, only R _{ B3 }is greater than 1 and our solution goes to the endemic equilibrium (third row, Table 1). Hence increasing detection through means other than contact tracing might be more effective, although perhaps also more costly.
Simulations for model with first level of contact tracing only after initial time t = 0
Parameters values  Reproduction numbers  Limiting steady state  

k _{1}  k _{2}  k _{3}  R _{0}  R _{ B1}  R _{ B2}  R _{ B3}  R* 

0.25  0.25  0  1.19  1.31  NA  NA  0.60  endemic steady state E* 
0.25  0.85  0  0.54  0.99  NA  NA  1.02  DFE E _{ B1} 
Simulations for model without contract tracing
Parameters values  Reproduction numbers  Limiting steady state  

k _{1}  k _{2}  k _{3}  R _{0}  R _{ B1}  R _{ B2}  R _{ B3}  R* 

0.25  0  0  2.38  NA  NA  NA  NA  endemic steady state 
0.6  0  0  0.99  NA  NA  NA  NA  DFE E _{ B3} 
Discussion and Conclusions
Our results indicate the following
1. Increasing the efficiency of contact tracing (higher contact tracing rates k_{2} and k_{3}) is less effective than other means of detection (e.g., random screening), but is perhaps less costly and more costeffective, and hence is ideal for largescale intervention programs in developing countries with limited resources.
2. Without contact tracing (k_{2} = k_{3} = 0), detection by other means must be increased substantially (more than twofold in Table 3) to lower R _{0} to less than 1.
3. If the second level contact tracing is not carried out (k_{3} = 0), the proportion of undetected infectives x* will increase by about 22% (from 0.331 to 0.404) in our simulation. Moreover, drastic increase (by more than threefold) in the first level contact tracing rate is need for the system to become disease free.
4. The results in Fig. 4 indicate that if random screening had been highly effective (90% in our simulation), implementation of contact tracing program could help to eradicate the disease in the sense that the endemic fractions (y _{1}, y _{2}, and y _{3}) will go to 0 in time. However, if the random screening does not detect a sufficiently high percentage of infectives, the implementation of the same program will not prevent the disease from becoming endemic.
5. In terms of dynamics, the model with only one level of contact tracing (Hsieh et al. 2005) exhibits very simple dynamics where all solutions either approach the DFE or the endemic equilibrium. When secondary level of contact tracing is included, bistability might occur where the solution could approach either the endemic steady state or any of the DFEs (E _{0} or any one of the diseasefree boundary steady states), as illustrated in rows 3 and 4 of Table 1 (and Fig. 4), depending on the initial proportions of population sizes. In other words, whether the contact tracing is effective, or how effective the program is, might depend on the initial state of the epidemic in which a particular country or region is when intervention by contact tracing is first being implemented. This has important public health implications as it has been suggested in [4] that with antiretroviral therapy being more widely available, intervention policy based on intensive HIV testing and tracing of partners could be considered as a possible policy to control HIV/AIDS epidemics in other countries (also see [27]). Although we note that due to our assumption of low HIV prevalence and simulations using an initial proportion of undiagnosed infectives at 30% (or 70% detection rate of the HIVpositive individuals), the results might not be applicable to countries with low detection rate.
6. Finally, the linear force of infection term employed in this work can be considered as a linear approximation of the more commonly used standard incidence, λXS/(S + X + Y_{1} + Y_{2} + Y_{3}), when the total number of HIVinfectives, X + Y_{1} + Y_{2} + Y_{3}, is small compared with S.
A stochastic Markov process model was proposed recently to describe the contact tracing detection in Cuba, in order to compare it with the usual random screening detection [28]. In this model, the component accounting for the contacttracing feature is assumed to be valued in a space of point measures in order to take the time since detection into account. When an individual is identified as infected, it may contribute to detecting other infectious individuals by providing information related to persons with whom they have had possibly infectious contacts. This was expected to reinforce standard randomscreening. Three models with three forms of the detection by contacttracing have been considered by [28], where statistical estimation of the two rates of detection (by random detection and by contact tracing) using Cuban HIV/AIDS data shows that whatever the model used, the contact tracing detection reaches an efficiency equivalent to that of random screening detection approximately four years after the beginning of the epidemic. This fact, which had not been noticed by epidemiologists and health practitioners, is noteworthy since it goes countercurrent against the classical view on contacttracing that it is efficient at the very beginning of its implementation and then less and less useful, insofar as the information it provides becomes redundant and obsolete as the contact tracing within a sexual network continues. More recently, approximate Bayesian computation is also applied to the Cuban HIV/AIDS data with contacttracing and unobserved infectious population to make further observations concerning the efficiency of HIV detection system in Cuba [29].
Our results show that even if the secondary level of contact tracing is less effective in detection, with a substantial plunge in the percentage of detection from the first level contact tracing as noted previously, it still plays an important role in determining whether the disease can eventually be eradicated. This conclusion corroborates with the above conclusion by [28] on the effectiveness of contact tracing even a period of time after its implementation and even when the information it provides might become redundant and obsolete as we continue through the tertiary stages of contact tracing. Our results further indicate that contact tracing for HIV could be suitably used to remedy inadequacies in a universal testing program and in design of timely and effective intervention. Perhaps some combination of periodic (universal) random screening combined with contact tracing would be the most ideal and effective program to achieve the goal of quick detection and to increase early treatment, in order to attain eventual elimination [27].
As a final remark, the Family Doctor Program in Cuba, introduced in 1984, also provides a strong communityoriented primary care network which contributes to improving population health overall [30] and to respond to emerging and reemerging diseases (including HIV), which also provide a frontline for prevention efforts. The family doctors' offices and community polyclinics play a pivotal role in prevention by providing counseling services, active screening services, promoting safer sex practices, and distributing educational materials, condoms, and lubricants, while also serving as the primary care providers for people with HIV participating in the ambulatory care system. Taken together, these initiatives profoundly affected how HIV manifested in Cuba [31] where in recent years the proportion of detections by contact tracing has fallen to 15% of those detected and some persons were detected through the family doctors before they could be traced by the contact tracing system. This is due to the fact that the number of reported contacts per detection has remained on the same level as before but as detection has almost doubled, the system could not effectively and swiftly cope with the larger number of contacts. Moreover, the family doctors are not really a random search or screening, since typically the family doctors counsel someone to test for HIV because that they suspect he/she might be infected. Therefore this element of detection is not readily modeled in our present model. In [28], it is found that a screening through contact tracing could reach peak efficiency after 4 years, and this dropoff conceivably could also apply to the family doctors program, as it is similarly not a random search. Perhaps a worthwhile future model extension, requiring detailed and updated data on the detection through the family doctor program, would include such consideration for understanding the detection and surveillance of HIV.
Declarations
Acknowledgements
YHH and YSW are supported by National Science Council (NSC) of Taiwan under grants NSC972115M039002MY3 and 972314B039 013MY3. HdeA is supported by the Spanish Agency for International Cooperation (AECID) projects A01010807 and D984207. RL and HdeA are supported by the Agence Nationale de la Recherche of France (ANR) under grant ANR08SYSC01604. The authors are grateful to the reviewers for their constructive and insightful comments which significantly improved this paper.
Authors’ Affiliations
References
 Granich RM, Gilks CF, Dye C, De Cock KM, Williams BG: Universal voluntary HIV testing with immediate antiretroviral therapy as a strategy for elimination of HIV transmission: a mathematical model. Lancet. 2008, 373: 4857. 10.1016/S01406736(08)616979.View ArticlePubMedGoogle Scholar
 Inciardi JA, Syvertsen JL, Surratt HL: HIV/AIDS in the Caribbean Basin. AIDS Care. 2005, 17 (Suppl 1): S925. 10.1080/09540120500120955.View ArticlePubMedGoogle Scholar
 UNAIDS: Report on the global AIDS epidemic, August 2008. 2009, Geneva: UNAIDS, Accessed January 5, [http://www.unaids.org/en/KnowledgeCentre/HIVData/GlobalReport/2008/2008_Global_report.asp]Google Scholar
 de Arazoza H, Joanes J, Lounes R, Legeai C, Clemencon S, Perez J, Auvert B: The HIV/AIDS epidemic in Cuba: description and tentative explanation of its low HIV prevalence. BMC Infec Dis. 2007, 7: 13010.1186/147123347130.View ArticleGoogle Scholar
 Hsieh YH, de Arazoza H, Lee SM, Chen CW: Estimating the number of Cubans infected sexually by human immunodeficiency virus using contact tracing data. Int J Epidemiol. 2002, 31 (3): 679683. 10.1093/ije/31.3.679.View ArticlePubMedGoogle Scholar
 de Arazoza H, Lounes R, Pérez J, Hoang T: What percentage of the Cuban HIVAIDS Epidemic is known?. Revista Cubana de Medicina Tropical. 2003, 55 (1): 3037.PubMedGoogle Scholar
 Hsieh YH, Wang HC, de Arazoza H, Lounes R, Twu SJ, Hsu HS: Ascertaining HIV underreporting in low HIV prevalence settings. J Biol Systems. 2005, 13 (4): 441454. 10.1142/S0218339005001616.View ArticleGoogle Scholar
 Lounes R, de Arazoza H, Hsieh YH, Joanes J: Deterministic modeling on the size of the HIV/AIDS epidemic in Cuba. Parametric Mathematical Methods in Survival Analysis Reliability and Quality of Life. Edited by: Huber C, Limnios N, Mesbah M, Nikulin M. 2008, ISTE and J.Wiley London, Chap 20: 317330.Google Scholar
 Hsieh YH, Lee SM, Chen CWS, de Arazoza H: On the recent sharp increase of HIV infections in Cuba. AIDS. 2001, 13 (3): 425428.Google Scholar
 Lounes R, de Arazoza H: A twotype model for the Cuban national programme on HIV/AIDS. IMA J Math Appl Med Biol. 1999, 16: 143154. 10.1093/imammb/16.2.143.View ArticlePubMedGoogle Scholar
 de Arazoza H, Lounes R, Hoang T, Interlan Y: Modeling HIV epidemic under contact tracingthe Cuban case. J Theor Med. 2000, 2: 267274. 10.1080/10273660008833053.View ArticleGoogle Scholar
 de Arazoza H, Lounes R: A nonlinear model for a sexually transmitted disease with contact tracing. Math Med Bio. 2002, 19: 221234. 10.1093/imammb/19.3.221.View ArticleGoogle Scholar
 Hsieh YH, de Arazoza Rodríguez H, Lounes R: A Class of Models for HIV Contact Tracing in Cuba: Implications for Intervention and Treatment. Deterministic and Stochastic Models for AIDS Epidemics and HIV Infection with Interventions. Edited by: Tan WY, Wu HL. 2005, Singapore: World ScientificGoogle Scholar
 Hsieh YH: An AIDS model with screening. Math Comp Mod. 1990, 14: 640643. 10.1016/08957177(90)90260T.View ArticleGoogle Scholar
 Hsieh YH: Modelling the effect of screening in HIV transmission dynamics. Differential Equations Models in Biology, Epiemiology and Ecology, Proc. Inter. Conf. on Diff. Eq. Claremont Lec. Notes in Biomath. 1991, SpringerVerlag, New York, 99120.View ArticleGoogle Scholar
 VelasoHernanadez JX, Hsieh YH: Modeling the effect of treatment and behavioral change in HIV transmission dynamics. J Math Biol. 1994, 32: 233429. 10.1007/BF00163880.View ArticleGoogle Scholar
 Hsieh YH, VelasoHernanadez J: Community treatment of HIV1: Initial and asymptotic dynamics. BioSystems. 1995, 35 (1): 7581. 10.1016/03032647(94)01482M.View ArticlePubMedGoogle Scholar
 Kretzschmar M: Sexual network structure and sexually transmitted disease prevention: a modeling perspective. Sex Transm Dis. 2000, 27 (10): 62735. 10.1097/0000743520001100000011.View ArticlePubMedGoogle Scholar
 Ferguson N, Garnett G: More realistic model of sexually transmitted disease transmission dynamics: sexual partnership networks pair models and moment closure. Sex Transm Dis. 2000, 27 (10): 600609. 10.1097/0000743520001100000008.View ArticlePubMedGoogle Scholar
 Kretzschmar M, Welte R, van den Hoek A, Postma MJ: Comparative Modelbased Analysis of Screening Programs for Chlamydia trachomatis Infections. Am J Epidemiol. 2001, 153 (1): 90101. 10.1093/aje/153.1.90.View ArticlePubMedGoogle Scholar
 Bauch CT: A versatile ODE approximation to a network model for the spread of sexually transmitted diseases. J Math Bio. 2002, 45 (5): 375395. 10.1007/s002850200153.View ArticleGoogle Scholar
 Keeling MJ, Eames KTD: Networks and epidemic models. J Roy Soc Interface. 2005, 2 (4): 295307. 10.1098/rsif.2005.0051.View ArticleGoogle Scholar
 Turner KME, Adams EJ, Gay N, Ghani AG, Mercer C, Edmunds WJ: Developing a realistic sexual network model of chlamydia transmission in Britain. Theoretical Biology and Medical Modelling. 2006, 3: 310.1186/1742468233.View ArticlePubMedPubMed CentralGoogle Scholar
 Bansal S, Grenfell BT, Meyers LA: When individual behaviour matters: homogeneous and network models in epidemiology. J Roy Soc Interface. 2007, 4 (16): 876891. 10.1098/rsif.2007.1100.View ArticleGoogle Scholar
 Eames KTD: Contact tracing strategies in heterogeneous populations. Epidemiology and Infection. 2007, 135 (3): 443454. 10.1017/S0950268806006923.View ArticlePubMedGoogle Scholar
 Arino O, Abdllaoui AE, Mikaram J, Chattopadhyay J: Infection in prey population may act as a biological control in ratiodependent predatorprey models. Nonlinearity. 2004, 17: 11011116. 10.1088/09517715/17/3/018.View ArticleGoogle Scholar
 Hsieh YH, de Arazoza H: Correspondence to "Universal voluntary HIV testing and immediate antiretroviral therapy". Lancet. 2009, 373: 10791080. 10.1016/S01406736(09)606450.View ArticlePubMedGoogle Scholar
 Clemençon S, Tran VC, de Arazoza H: A stochastic SIR model with contacttracing: large population limit with statistical applications. J Bio Dynamics. 2008Google Scholar
 Blum MGB, Tran VC: Approximate Bayesian Computation for epidemiological models: Application to the Cuban HIVAIDS epidemic with contacttracing and unobserved infectious population. Arxiv. 2008, arXiv: 0810.0896,Google Scholar
 Presno Labrador C, Sansó Soberat F: 20 years of family medicine in Cuba. MEDICC Review. 2004, 6 (2):
 Gorry C: Cuba's HIV/AIDS Strategy: An Integrated RightsBased Approach. 2008, Oxfam International and Programa Conjunto de Oxfam, [http://www.medicc.org/ns/assets/documents/Cuban%20HIV%20Strategy.pdf]Google Scholar
 The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712334/10/194/prepub
Prepublication history
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.