HIV epidemic in Far-Western Nepal: effect of seasonal labor migration to India
© Vaidya and Wu; licensee BioMed Central Ltd. 2011
Received: 24 November 2010
Accepted: 13 May 2011
Published: 13 May 2011
Because of limited work opportunities in Nepal and the open-border provision between Nepal and India, a seasonal labor migration of males from Far-Western Nepal to India is common. Unsafe sexual activities of these migrants in India, such as frequent visits to brothels, lead to a high HIV prevalence among them and to a potential transmission upon their return home to Nepal. The present study aims to evaluate the role of such seasonal labor-migration to India on HIV transmission in Far-Western Nepal and to assess prevention programs.
An HIV epidemic model was developed for a population in Far-Western Nepal. The model was fitted to the data to estimate the back and forth mobility rates of labor-migrants to India, the HIV prevalence among migrants and the HIV transmission rate in Far-Western Nepal. HIV prevalence, new infections, disease deaths and HIV infections recruited from India were calculated. Prevention programs targeting the general population and the migrants were evaluated.
Without any intervention programs, Far-Western Nepal will have about 7,000 HIV infected individuals returning from India by 2015, and 12,000 labor-migrants living with HIV in India. An increase of condom use among the general population from 39% to 80% will reduce new HIV infections due to sexual activity in Far-Western Nepal from 239 to 77. However, such a program loses its effectiveness due to the recruitment of HIV infections via returning migrants from India. The reduction of prevalence among migrants from 2.2% to 1.1% can bring general prevalence down to 0.4% with only 3,500 recruitments of HIV infections from India.
Recruitment of HIV infections from India via seasonal labor-migrants is the key factor contributing to the HIV epidemic in Far-Western Nepal. Prevention programs focused on the general population are ineffective. Our finding highlights the urgency of developing prevention programs which reduce the prevalence of HIV among migrants for a successful control of the HIV epidemic in Far-Western Nepal.
Since the first case in July 1988 [1–3], the reported number of HIV infections in Nepal has gradually increased to 10,546 HIV cases and 1,610 AIDS cases as of December 2007 . Until the late 1990s, Nepal was classified as having a low-level epidemic. However, after 1997, Nepal has been experiencing a concentrated epidemic with rapid spread amongst high-risk groups [2–4]. Even though the HIV prevalence among the general population has been estimated to be low (~0.5%) [4, 5], among the high-risk groups it varies from 1-3% among female sex workers (FSWs) in Kathmandu (Capital of Nepal), 2.8% among migrants returning from Mumbai, India, 3.3% among urban-based males having sex with males (MSMs) and 10% among male migrants returning to Doti (a Far-Western district of Nepal) to 34.7-58% among injection drug users (IDUs) in Kathmandu [2, 4, 6–9]. Among these HIV risk factors, a high rate of seasonal migration of the male population to India has recently been extremely threatening in the Far-Western districts of Nepal such as Doti, Achham, Dang, Kanchanpur and Kailali [4, 7–9]. Nepal's 2007 United Nations General Assembly Special Session (UNGASS) report indicates that labor migrants make up 41% of the total known HIV infections in the country .
Because of limited work opportunities in Nepal and the open-border provision between Nepal and India, 1.5 to 2 million Nepalis have been estimated to migrate to India for seasonal and long-term work. These especially include labor migrants from Far-Western Nepal [4, 8–10]. 50-80% of households in some communities of this region have at least one family member working in India, and most of these migrants seasonally return home [11, 12]. Their major destinations are Mumbai, Panjab and Chenai in India, where HIV prevalence among sex workers is high (for example, there was 71% prevalence among the sex workers of Mumbai in 1997) [7, 9, 13]. Their unsafe sexual activities in India - for example, frequent brothel visits, mostly after drinking alcohol [7, 8] - make for a high HIV prevalence among these migrants group. Upon their return home, in addition to transmitting HIV to their wives, they also may transmit to other women because their new social status affords them the chance to have extramarital sex, for example during frequently occurring local festivals and Deuda (a popular local cultural event) [7, 8, 12].
While there are many existing mathematical models developed to understand the dynamics of HIV/AIDS in different parts of the world [14–19] and there are some field surveys [7, 8, 20, 21], no attempt has been made to mathematically study the epidemic of HIV/AIDS in Far-Western Nepal. In this study, we developed a mathematical model to evaluate effects of seasonal labor migrants to India on the spread of HIV/AIDS in Far-Western Nepal. We estimated key parameters by fitting our model to the estimated data. We also assessed two prevention programs, one targeted at the general population and another at the migrants.
Estimated adult (15-49 age group) HIV prevalence % in Far-Western Nepal during the year from 1990 to 2007.
The parameters α and λ are the rates of migration to India and from India, respectively. Returned migrants enter into the HIV infected and susceptible compartments with rates p and 1-p, respectively. Here, p denotes the portion of HIV positive individuals among migrants in India. μ and δ denote the natural death rate and the disease death rate, respectively. We take Λ as a rate of new entries into the population of interest, i.e., the maturation rate for sexual activity. We model disease transmission and population changes over time using a system of nonlinear differential equations as follows:
Here, the HIV transmission rate, β, is a function of the average number of sexual partners, n, the average condom use rate, c, and infectivity (i.e., the probability per unprotected sexual partnership that an infected individual transmits the disease to a susceptible individual), θ, according to formula: β = n(1-c)θ .
Note that our model can further be extended by dividing the total migrant group, M, into susceptible and infected groups. However, among the male migrants considered in this study, a new HIV infection in India is not due to sexual contacts between susceptible and infected migrants, but is due to unsafe sexual activities such as frequent brothel visits [7, 8]. Therefore, the model having susceptible and infected migrant groups separately will require simply a linear term, rather than non-linear term, representing a flow with a constant rate from the susceptible migrant group to the infected migrant group. In this case, considering a single parameter, as p in our model, to represent the HIV prevalence among migrants can produce equivalent results. Moreover, a representation of HIV prevalence among migrants as a single parameter makes it easier to compare our estimates with the field survey studies.
See the Appendix for the expression of HIV prevalence in a long-term and further model analysis.
Even though Nepal is facing the challenge of controlling HIV among its high-risk groups, the disease prevalence among the general population remains low (less than 0.5% in 2007 [4, 5]). Therefore, the model simplified to a low prevalence scenario can reasonably approximate the dynamics of HIV/AIDS in Far-Western Nepal. In the low prevalence setting, we can solve the model to obtain close form solutions S, I and M (See Appendix).
Nepal's life expectancy at birth is 60.56 years . With 15-year as the age of entering into the sexually active population, a remaining life-time of the healthy population following recruitment to the sexually active group is 60.56-15 = 45.56 years. This gives 1/μ = 45.56 years, i.e. the natural death rate is μ = 0.022 y-1. Realizing that the life-expectancy at age 15 is usually higher than at birth due to a higher infancy and childhood death rate in Nepal, we performed a sensitivity analysis by taking 70-year and 80-year life expectancies. The predicted HIV prevalence increased by 0.03% and 0.05% for 70-year and 80-year life expectancies, respectively. This shows that the predicted prevalence is less sensitive to μ, and we present results with μ = 0.022 y-1.
The estimated mean duration of infection without treatment varies widely from 8.6 years to 19 years [17, 19, 32, 33]. The HIV/AIDS response in Nepal has largely focused on prevention and awareness, with less attention to the treatment, care, and support of persons living with HIV/AIDS . For example, by the end of 2003, estimated 4,000 adults in Nepal needed antiretroviral therapy, but only 77 received the treatment as of June 2004. Since the health care system is poor and the education level is low, the life expectancy of HIV infected people is relatively short. Therefore, we assumed that the mean infection period is 8.6 years, which gives 1/δ = 8.6 years, i.e. δ = 0.116 y-1.
We fitted both the close form solution and the full SIM model to the HIV prevalence data and estimated parameters λ, α, p, β and Λ (See Appendix for a detailed description of data fitting and calculation of 95% confidence intervals).
Model parameters for Far-Western Nepal.
Base Value [Low-High]
Total sexually active population
Estimated, Calculated, 
Migrants in India
Death rate (non-HIV)
HIV+ probability in returned migrants
HIV Disease Parameters
HIV/AIDS death rate
0.116 [0.116 - 0.053]
Sexual Behavior Parameters
Annual no. of sex partners
HIV transmission probability
We estimated the HIV prevalence among the returning migrants as 2.2% (95% CI: 2.0-2.3), consistent with a finding of 2.8% prevalence among the Mid-Far Western migrants returning from Mumbai, India . We note that there are survey reports which indicate higher HIV prevalence among some returning migrants [7–9]. However, those studies were focused on the specific areas of Doti and Achham districts, and so do not represent overall Far-Western Nepal.
According to the World Factbook , the birth rate of Nepal is 30.46 births/1,000 population and the infant mortality rate is 63.66 deaths/1,000 live births. These provide the maturation rate of 2.85% of the population, which is in accordance with our estimate of Λ = 15,321 (95% CI: 15,303-15,339), i.e. 2.40% of the total population.
Among non-migrants surveyed by Poudel et al. , 75% did not have multiple sex partners in the last 5 years. Taking one sex partner for 75% of the population and on average 3 sex partners for the remaining 25% of the population, we get the average number of sex partners to be n = 1.5. Some studies show that the rate of condom use in Far-Western Nepal varies from 12% by married migrants with their wives  and 44% by non-migrants  in pre- or extra-marital sex to 60% by some migrants with sex workers . Taking an average, the rate of condom use is 38.7% i.e. c~0.39. Using this value in our estimate of β = 0.02 (95% CI: 0.01-0.03), we obtained the transmission probability θ = 0.02, which is within the limit of 0.01-0.06 given in [17, 33, 34].
We obtained the rate of mobility to India as α = 0.39 (95% CI: 0.36-0.41) and the returning rate as λ = 0.024 (95% CI: 0.020-0.028). In the early 90s the seasonal migration to India was relatively low, as a 1994 survey conducted in 11 districts of the Far-Western and the Mid-Western regions found that 15% of adults migrated seasonally to India . However, a decade of violent conflict (Maoist insurgency 1996-2006) dramatically exacerbated migration and displacement [36, 37]. As a consequence, a recent estimate suggests an out-migration of up to 90% of adult men in some villages of the Far-Western hills [24, 38]. The high mobility to India is also in agreement with the fact that one record in late December 2004 puts the figure at 200 Nepalis crossing the border every hour .
Disease Outcomes: Prevalence and Infection
Effect of two prevention programs - PTG and PTM - on disease outcomes.
Base case (39% condom use and p = 2.2%)
80% condom use
p = 1.1%
p = 0.5%
Effect of Seasonal Labor Migration on Disease Outcomes
If the out-flow migration to India increases from the current rate of 39% to 60%, HIV prevalence in 2015 increases to 0.91%. In this case, 13,000 HIV infected migrants will be living in India in the year 2015 with the total recruitment of 7,500 HIV infections into Far-Western Nepal. Similarly, an increase in the rate of in-flow migration from India to 5% causes an increase of HIV prevalence to 1.11% with about 14,000 HIV infections entering Far-Western Nepal and 12,000 migrants living with HIV/AIDS in India in the year 2015.
HIV prevalence in the year 2015 among the general population in Far-Western Nepal increases linearly with a rate 0.32% per % increase in p, which represents the unsafe sexual activity level of migrants in their work place in India (Figure 3c). If the prevalence among the returning migrants from India reaches as high as 10% obtained in a 2001 survey conducted in some villages of Doti district [7, 8], the prevalence among the general population will reach about 3.5% in the year 2015. This will project into a scary figure of 56,000 Far-Western migrants living with HIV/AIDS in their work place in India with a recruitment of 32,000 HIV infections to Far-Western Nepal and 159,000 deaths due to AIDS by the year 2015.
Effect of Back Home Sexual Activity on Disease Outcomes
The sexual activity back home in the Far-Western region can be evaluated with the parameter β. In the absence of HIV infection due to migration, β also defines one of the important disease dynamic properties, the basic reproductive number, R 0 = β/(α+ μ+δ) [40, 41] (See Appendix), which represents the effective number of secondary infections caused by a typical infected individual during his/her entire period of infectiousness . The current basic reproductive number in Far-Western Nepal is R 0 = 0.03 < 1 indicating that the disease will theoretically be under control provided that the recruitment of HIV infection via the returning migrants is completely stopped. A change of sexual behavior resulting in an increase of β by 10-fold i.e. β = 0.18 leads to an increase of HIV prevalence to 1.22% in the year 2015 with an increase of the total new infections to 3,200 and the total death due to AIDS to 36,000 (Figure 3d).
Effect of Prevention Programs on Disease Outcomes
We evaluated two illustrative prevention programs: one targeted at the general population (PTG) and another at the migrants (PTM). While PTG is focused on the behavioral interventions in the Far-Western region such as educating villagers to increase condom use and/or reduce the number of sexual partners, PTM is focused on the behavioral interventions in India such as reducing the frequency of brothel visits by migrants, particularly after drinking alcohol, or increasing safe sex practices in India.
In our model, varying condom use is similar to varying the number of sexual partners. If the average condom use increases from the current 39% to 80%, HIV prevalence in the year 2015 can be reduced by 0.02% to the prevalence of 0.82%, which decreases the total new infections due to the sexual activities back home from 239 to 77. Despite a significant reduction of new infections due to the sexual activities back home, PTG remains ineffective in controlling the HIV epidemics. More importantly, PTG does not help suppress the HIV infections due to the recruitment via returning migrants from India. About 12,000 migrants will still be living with HIV/AIDS, with the recruitment of 7,000 HIV infections from India (Table 3).
PTM, on the other hand, is significantly effective in reducing HIV prevalence both among the migrants and the general population. By the effective intervention that reduces p from 2.2% to 1.1%, HIV prevalence among the general population can be reduced to 0.42% as compared to 0.84% of the base case. With p = 1.1% only 6,300 migrants will be living with HIV/AIDS in India in the year 2015. Moreover, the number of HIV infections recruited via the returning migrants, the total death due to AIDS and the total new infections due to the sexual activities back home reduce to 3600, 800 and 130, respectively (Table 3). By reducing p further, the disease epidemic can significantly be mitigated (See Table 3 and Figure 3c).
Among the four major HIV epidemic regions of Nepal, i.e., highway districts (49% of HIV infections), Kathmandu valley (16%), Far-Western region (16%) and remaining hills (19%), the Far-Western region is unique in terms of the recruitment of HIV infections via seasonally returning migrants from India. Here we developed a mathematical model to evaluate the effects of seasonal labor migration to India on exacerbating the HIV burden in the Far-Western region of Nepal. Using the calibrated model, we assessed two prevention strategies - one focusing on back home and another on migrants in their work place in India.
Based on the USAID/WHO data for the whole country of Nepal, we derived HIV prevalence in the Far-Western region. By fitting our model to this derived data, we estimated the mobility rate of seasonal migrants to and from India, the rate of HIV prevalence among migrants, the transmission rate due to sexual activities back home and the maturation rate of the population. We calculated the number of HIV infected Far-Western migrants living in India, new infection due to back home sexual activities, HIV infections recruited from India via returning migrants and deaths due to AIDS by the year 2015.
The estimation from data fitting using our model predicts a high mobility rate of seasonal labor migrants from Far-Western Nepal to India. This high mobility is related to the most economic deteriorated condition of the Far-Western region and 10-year long violent conflicts (Maoist Insurgency 1996-2006) within the country. The Far-Western region has been the most affected region by this economic hardship, leading directly to a huge displacement of adult population to India [24, 36–39], particularly to Mumbai. There is evidence of unsafe sexual activities of these migrants while they are working in India - for example, frequent brothel visits, mostly after drinking alcohol [7, 8]. Such activities put these migrants at a high risk of contracting HIV and subsequently transmitting it to the Nepalese population upon their return. Our model predicts the HIV prevalence among these returning migrants to be 2.2%, which is consistent with a previous study . However, this rate could be as high as 10% if the study is focused on some specific communities such as some Doti villages and Achham villages of the Far-Western region, where the vulnerability of HIV infections due to labor migrants is quite high [7–9]. The HIV prevalence among returning migrants can significantly exacerbate HIV prevalence among the general population (Figure 3c).
Compared to the total number of HIV infections recruited from India via seasonally returned migrants, the total number of new infections generated due to back home sexual activities is significantly low. This can also be reflected by the low reproductive number of 0.03 in the absence of the recruitment of HIV infections from India. This explains the current condition of low HIV prevalence of Nepal among the general population (less than 0.5% [4, 5]) even though the high HIV prevalence among high-risk groups such as male migrants is one of the biggest concerns of the country's public health.
Our results show that HIV cannot be eradicated from the Far-Western region of Nepal unless the recruitment of HIV infections from India is stopped. Since the major source of HIV infections is the returning migrants rather than back home sexual activities, the prevention programs that educate villagers to promote safe-sex practices such as condom use are ineffective in suppressing the HIV burden in the region. Moreover, given the high mobility across the open-borders between Nepal and India, implementing a screening program at the borders is an infeasible immediate policy. Our model suggests that prevention programs focused on migrants to promote safe-sex activities in their workplace will be most effective on suppressing the HIV epidemic in Far-Western Nepal. In a survey study by Poudel et al. , the practice of safer sex in Mumbai by migrant workers has been reported after they received a package containing a letter with pamphlets about HIV/AIDS and condom packets. Such prevention programs that reduce the HIV prevalence among migrants can significantly lower the HIV prevalence among the migrants as well as the general population, new infections and disease deaths (Table 3). Therefore, programs that reduce the prevalence of HIV among migrants are urgently needed for a successful control of the HIV epidemic in Far-Western Nepal.
It is worth mentioning that population mobility to and from a region has been a fundamental characteristic of today's globalized world. The effect of such mobility on increasing the HIV burden has been evidenced in many parts of the world such as South Africa [43–45] and Australia . Therefore, our model and analysis with necessary modifications might be applied to evaluate the HIV dynamics in other parts of the world where migration is playing a role in HIV epidemics.
Mathematical models developed based on the disease characteristic in a particular region help better understanding the disease dynamics of that region and guide policy-makers in allocating resources for prevention and control of the infectious disease epidemic. We developed a mathematical model to explain the dynamic of HIV epidemics in Far-Western Nepal. Our SIM model suggests that the seasonal labor migration to India and their sexual activities in their workplace are the key factors contributing to the dynamics of the HIV epidemic in Far-Western Nepal. Any prevention strategy focused on activities of labor migrants in their workplace could provide an effective way of fighting against the HIV epidemic in Far-Western Nepal.
Stable Equilibrium Prevalence
For p = 0 (i.e. in the absence of the recruitment of HIV infection from India via seasonally returned migrants), R 0 represents the basic reproductive number, which is defined as the effective number of secondary infections caused by a typical infected individual during his/her entire period of infectiousness.
Data Fitting and Confidence Interval
Here, is the HIV prevalence data at time point t i given in Table 1. We used t 1990 = 0, the initial time for the computation. The parameter Λ can be approximated by using the mean-value, where t f is the final time-point considered.
Furthermore, we estimated the parameters λ, α, p, β and Λ by fitting the full SIM model (1) - (3) to the data (Table 1). Since the estimated parameters (Table 1) and the predicted curves (Figure 2) by both SIM model (1) - (3) and the solution (20) - (22) are almost the same, the expressions (20) - (22) provide a reasonably good approximation to the solution of the SIM model.
This work began when NKV was at York University. The work of NKV has been partially supported by Susan Mann Scholarship Award of York University, Canada and Postdoctoral Research Program at Los Alamos National Laboratory, USA. The work of JW has been supported in part by The Canada Research Chair Program, by Natural Sciences and Engineering Research Council of Canada, and by Mathematics for Information Technology and Complex Systems. Authors would like to thank the reviewers for useful comments.
- Suvedi BK: Transition of HIV epidemic in Nepal. Kathmandu Univ Med J (KUMJ). 2006, 4 (1): 115-118.Google Scholar
- Health profile: Nepal, HIV/AIDS. 2005, USAID
- HIV/AIDS in Nepal. 2008, World Bank
- HIV/AIDS health profile: Nepal. 2008, USAID
- Epidemiological factsheet on HIV and AIDS. Core data on epidemiology and response. Fact Sheet: Nepal. http://www.searo.who.int/LinkFiles/Facts_and_Figures_EFS2008_Nepal.pdf, http://www.unaids.org/en/regionscountries/countries/nepal/ accessed: July 10, 2009
- Integrated Bio-Behavioral Survey among Male Labor Migrants in 11 Districts in Western, and Mid-Far Western Regions of Nepal. 2006, New Era
- Poudel KC, Okumura J, Sherchand JB, Jimba M, Murakami I, Wakai S: Mumbai disease in far western Nepal: HIV infection and syphilis among male migrant-returnees and non-migrants. Trop Med Int Health. 2003, 8 (10): 933-939. 10.1046/j.1365-3156.2003.01110.x.View ArticlePubMedGoogle Scholar
- Poudel KC, Jimba M, Okumura J, Joshi AB, Wakai S: Migrants' risky sexual behaviours in India and at home in far western Nepal. Trop Med Int Health. 2004, 9 (8): 897-903. 10.1111/j.1365-3156.2004.01276.x.View ArticlePubMedGoogle Scholar
- National report: Nepal. 2005, UNGASS
- Pick B, Rogosch J, Thompson C, Wolf C: Two year review of the USAID/Nepal HIV/AIDS strategy 2001-2006. 2004Google Scholar
- Furber AS, Newell JN, Lubben MM: A systematic review of current knowledge of HIV epidemiology and of sexual behaviour in Nepal. Trop Med Int Health. 2002, 7 (2): 140-148. 10.1046/j.1365-3156.2002.00844.x.View ArticlePubMedGoogle Scholar
- Poudel KC, Jimba M, Okumura J, Sharma M, Poudel Tandukar K, Wakai S: Migration in far western Nepal: a time bomb for a future HIV/AIDS epidemic?. Trop Doct. 2004, 34 (1): 30-31.PubMedGoogle Scholar
- Schwartlander B, Garnett G, Walker N, Anderson R: AIDS in a new millennium. Science. 2000, 289 (5476): 64-66. 10.1126/science.289.5476.64.View ArticlePubMedGoogle Scholar
- Gani J, Yakowitz S, Blount M: The spread and quarantine of HIV infection in a prison system. SIAM J Appl Math. 1997, 57: 1510-1530. 10.1137/S0036139995283237.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): 679-683. 10.1093/ije/31.3.679.View ArticlePubMedGoogle Scholar
- Kgosimore M, Lungu EM: The effects of vertical transmission on the spread of hiv/aids in the presence of treatment. Math Biosci Eng. 2006, 3 (2): 297-312.View ArticlePubMedGoogle Scholar
- Long EF, Vaidya NK, Brandeau ML: Controlling Co-Epidemics: Analysis of HIV and Tuberculosis Infection Dynamics. Oper Res. 2008, 56 (6): 1366-1381. 10.1287/opre.1080.0571.View ArticlePubMedPubMed CentralGoogle Scholar
- Rapataski B, Klepac P, Dueck S, Liu M, Weiss LI: Mathematical epidemiology of hiv/aids in cuba during the period 1986-2000. Math Biosci Eng. 2006, 3 (3): 545-556.View ArticleGoogle Scholar
- Kribs-Zaleta CM, Lee M, Roman C, Wiley S, Hernandez-Suarez CM: The effect of the hiv/aids epidemic on africa's truck drivers. Math Biosci Eng. 2005, 2 (4): 771-788.View ArticlePubMedGoogle Scholar
- Poudel KC, Jimba M, Joshi AB, Poudel-Tandukar K, Sharma M, Wakai S: Retention and effectiveness of HIV/AIDS training of traditional healers in far western Nepal. Trop Med Int Health. 2005, 10 (7): 640-646. 10.1111/j.1365-3156.2005.01443.x.View ArticlePubMedGoogle Scholar
- Suvedi BK: HIV-TB co-infection in Nepal. J Inst Med. 2003, 25: 19-21.Google Scholar
- Epidemiological fact sheet on HIV/AIDS and sexually transmitted infections. 2006, Fact Sheet: Nepal
- UNGASS Country Progress Report [Nepal] (Reporting Period: January 2006-December 2007). 2008, UNGASS, (accessed: January 17, 2010), [http://data.unaids.org/pub/Report/2008/nepal_2008_country_progress_report_en.pdf]
- Nepal B: Population mobility and spread of HIV across the Indo-Nepal border. J Health Popul Nutr. 2007, 25 (3): 267-277.PubMedPubMed CentralGoogle Scholar
- Prasai SB: Stop AIDS, halt trafficking in Nepalese women. Am Chron. 2008, (Access date: November 14, 2008), [http://www.americanchronicle.com/articles/view/53039]Google Scholar
- Regions of Nepal. 2008, Statoids, (accessed: November 17, 2008), [http://www.statoids.com/unp.html]
- CBS: Central Bureau of Statistics, Nepal. 2001Google Scholar
- Keyfitz N, Flieger W: World Population Growth and Aging. Demographic Trends in the Late Twentieth Century. 1990, The University of Chicago Press, USAGoogle Scholar
- CBS: Central Bureau of Statistics, Nepal. 1995Google Scholar
- A survey of drug use patterns in western Nepal. [http://clinmed.netprints.org/cgi/content/full/2002050003v1]
- The World Factbook 2007 Nepal. 2007, United States Central Intelligence Agency (USCIA), Washington D.C., United States, (accessed: April 19, 2007), [https://www.cia.gov/library/publications/the-world-factbook/geos/np.html]
- Hyman JM, Li J, Stanley EA: The differential infectivity and staged progression models for the transmission of HIV. Math Biosci. 1999, 155 (2): 77-109. 10.1016/S0025-5564(98)10057-3.View ArticlePubMedGoogle Scholar
- Zaric GS, Barnett PG, Brandeau ML: HIV transmission and the cost-effectiveness of methadone maintenance. Am J Public Health. 2000, 90 (7): 1100-1111. 10.2105/AJPH.90.7.1100.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanders GD, Bayoumi AM, Sundaram V, Bilir SP, Neukermans CP, Rydzak CE, Douglass LR, Lazzeroni LC, Holodniy M, Owens DK: Cost-effectiveness of screening for HIV in the era of highly active antiretroviral therapy. N Engl J Med. 2005, 352 (6): 570-585. 10.1056/NEJMsa042657.View ArticlePubMedGoogle Scholar
- Suvedi BK, Gurubacharya VL, Thapa K: Seasonal migration and its relation to HIV transmission in western Nepal. International Conference on AIDS 10, 60 (abstract no 192D). 1994Google Scholar
- Singh S: Nepal's war and conflict-sensitive development. PLoS Med. 2005, 2 (1): e29-10.1371/journal.pmed.0020029.View ArticlePubMedPubMed CentralGoogle Scholar
- Singh S, Mills E, Honeyman S, Suvedi BK, Pant NP: HIV in Nepal: is the violent conflict fuelling the epidemic?. PLoS Med. 2005, 2 (8): e216-10.1371/journal.pmed.0020216.View ArticlePubMedPubMed CentralGoogle Scholar
- Global IDP Project: Profile of internal displacement: Nepal, compilation of the information available in the global IDP database of the Norwegian Refugee Council. 2005Google Scholar
- Bhattrai R: Open borders, closed citizenships: Nepali labor migrants in Delhi. International migration, multi-local livelihoods and human security: Perspectives from Europe, Asia and Africa: 2007; Institue of Social Studies, The Netherlands. 2007Google Scholar
- van den Driessche P, Watmough J: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002, 180: 29-48. 10.1016/S0025-5564(02)00108-6.View ArticlePubMedGoogle Scholar
- Heffernan JM, Smith RJ, Wahl LM: Perspectives on the basic reproductive ratio. J R Soc Interface. 2005, 2 (4): 281-293. 10.1098/rsif.2005.0042.View ArticlePubMedPubMed CentralGoogle Scholar
- Diekmann O, Heesterbeek JA, Metz JA: On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J Math Biol. 1990, 28 (4): 365-382.View ArticlePubMedGoogle Scholar
- Camlin CS, Hosegood V, Newell ML, McGrath N, Barnighausen T, Snow RC: Gender, migration and HIV in rural KwaZulu-Natal, South Africa. PLoS One. 5 (7): e11539-
- Jochelson K, Mothibeli M, Leger JP: Human immunodeficiency virus and migrant labor in South Africa. Int J Health Serv. 1991, 21 (1): 157-173. 10.2190/11UE-L88J-46HN-HR0K.View ArticlePubMedGoogle Scholar
- Nunn AJ, Wagner HU, Kamali A, Kengeya-Kayondo JF, Mulder DW: Migration and HIV-1 seroprevalence in a rural Ugandan population. Aids. 1995, 9 (5): 503-506.View ArticlePubMedGoogle Scholar
- Horyniak D, Stoove M, Yohannes K, Breschkin A, Carter T, Hatch B, Tomnay J, Hellard M, Guy R: The impact of immigration on the burden of HIV infection in Victoria, Australia. Sex Health. 2009, 6 (2): 123-128. 10.1071/SH08088.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/11/310/prepub