Factors associated with needle sharing among people who inject drugs in Yunnan, China: a combined network and regression analysis
© The Author(s). 2016
Received: 5 October 2015
Accepted: 7 July 2016
Published: 9 August 2016
Network analyses have been widely utilized to evaluate large datasets, but have not yet been used to explore factors associated with risk behaviours. In combination with traditional regression analysis, network analyses may provide useful information and highlight key factors for reducing needle sharing behaviours among people who inject drugs (PWID).
Sociodemographic data, and information on injection behaviour and sexual practices were collected from a cross-sectional survey that was conducted with PWID in five prefectures of Yunnan province, China. A combination of logistic regression and correlation network analyses were used to explore key factors for reducing needle-sharing behaviours among PWID.
In a total of 1 049 PWID, 37.5 % had a history of needle or syringe sharing. The logistic analysis showed that Zhaotong, Qujing, Dehong, or Lincang residents, diazepam use, longer injection duration, needle reuse, and infection with HIV, viral hepatitis, tuberculosis and/or malaria were independently associated with needle sharing. The correlation network analyses showed that, compared to PWID who had never shared needles, PWID who did share needles would achieve harm reduction goals faster and more permanently. HIV serostatus and marital status were found to be closely associated with other risk factors. By combining regression analyses with network analyses, it was shown that PWID who are HIV seropositive will be an ideal target group for harm reduction programs.
Needle-sharing behaviours are common among PWID in Yunnan, and harm reduction programs may help PWID who are HIV seropositive reduce risk behaviours and prevent blood borne diseases.
Please see Additional file 1 for translations of the abstract into the five official working languages of the United Nations.
Bordering north-eastern Myanmar, northern Laos, and northern Vietnam, Yunnan is a region that plays a crucial role in the transmission of blood borne infectious diseases in the Lancang-Mekong subregion. According to recent studies, drugs produced in the ‘Golden Triangle’ region, which consists of eastern Myanmar, western Laos, and northern Thailand, are exported from Myanmar to Yunnan . Along this drug route, subtypes B and C of the human immunodeficiency virus type 1 (HIV-1) are transmitted; subtypes B, C, and B/C recombinants are transmitted through the northwest of China, particularly through the Xinjiang province [1–4]. Studies have also shown that the circulating form of HIV, CRF01_AE, was introduced into Yunnan from northern Vietnam, likely through cross-border commercial sexual activities . In 2012, a study showed that unique recombinant forms of HIV-1 were transmitted to Myanmar from Yunnan in the early 1990s . These findings suggest that Yunnan may be a ‘transfer station’ of multidirectional transmission of HIV between China and its bordering countries.
People who inject drugs (PWID) are at a high risk of contracting blood borne diseases, and needle-sharing behaviours are believed to be the main route of transmission. The percentage of those who share needles among PWID ranges from 16 to 76 % in different countries, and is approximately 37 % in some countries [7–11]. Therefore, higher rates of needle or paraphernalia sharing can lead to an increased prevalence of blood borne diseases, including HIV, and hepatitis B and C virus infections [12, 13]. Although sharing of injecting equipment has been shown among PWID in Yunnan, there is a lack of up-to-date data in this field. As blood borne virus prevalence can be an indicator of needle and syringe sharing, a new evaluation of HIV prevalence in this population can be used as a measure of levels of needle and syringe sharing. Furthermore, reliable data on factors associated with needle and syringe sharing among PWID in Yunnan are scarce. Accurate information in this field is crucial for developing effective harm reduction strategies, reducing needle and syringe sharing among PWID, and lowering the disease burden of blood borne viruses.
Most previous studies have used traditional methods (including logistic regression) to investigate factors associated with needle and syringe sharing, and findings of these studies have contributed greatly to reducing needle-sharing behaviours among PWID. However, the interaction between each factor cannot be explored by traditional methods. Network analysis, which has been widely utilized in larger datasets in fields such as public health, social sciences, international relations, and protein dynamics, has contributed greatly to dealing with large datasets [14–17]. Network analyses can thoroughly and elegantly interpret large datasets into visual figures, show the inherent complexity of factor interactions, and explore the rate of influence that a given factor has in a correlation network. Therefore, using network analysis to analyse needle-sharing behaviours among PWID may provide useful information and highlight key factors for harm reduction programs.
In the present study, 1 049 PWID were recruited from Yunnan province in 2009; 20 sociodemographic and drug use behaviour factors were assessed. A combination of regression and network analyses was performed to evaluate needle-sharing characteristics, examine the possible predictors of needle sharing, compare the influence of each factor on PWID who partake and do not partake in needle-sharing behaviours, and explore key factors necessary for reducing needle-sharing behaviours.
Study participants and data collection
Basic characteristics of needle-sharing behaviors among PWID in Yunnan, China
Those who shared needles
Univariate logistic regression
Multivariate logistic regression
OR (95 % CI)
OR (95 % CI)
P - value
Place of residence
Divorced or widowed
> 10 years
Drug of choice
Heroin + otherse
Other drug use
Diazepam + othersf
Method of drug intake
Injection + oral
Duration of drug abuse
< 5 years
Drug use frequency (per day)
Most recent drug use
< 1 year
> 2 years
Duration of drug use (injection)
< 5 years
> 11 years
Source of syringe
Average reuse times per needle
Number of sexual partners
Diagnosis of other infectious diseaseh
Five milliliters of venous blood was collected from each study subject using vacutainer tubes containing EDTA-2 K. Blood samples were centrifuged and serum was stored at −80 °C in the freezer until HIV antibody testing was performed. The HIV serostatus was assessed using a diagnostic kit for antibody to HIV (1 + 2) (Colloidal Gold) (Alere Medical Co., Ltd, Chiba, Japan), and those who tested positive were confirmed using a diagnostic kit for antibodies to HIV (1 + 2) (enzyme-linked immunosorbent assay) (Wantai Medical Co., Ltd, Beijing, China), according to the manufacturers’ instructions.
All the data obtained from the questionnaires and previous laboratory tests were used in the following analyses. To explore the key factors for reducing needle-sharing behaviors among PWID, participants were divided into two subgroups: those who had never shared needles were allocated to the never sharing needles group (NSNG) and those who had a history of needle sharing were allocated to the sharing needles group (SNG) (see Fig. 1).
To compare the differences between PWID in different groups (NSNG and SNG), Student’s t-tests (two-tailed) were used for the factor of age and Mann–Whitney U tests (two-tailed) were used for the factors of gender, ethnicity, occupation, marital status, and education level.
A univariate logistic regression analysis was conducted to compare the percentage of those PWID who shared needles between the NSNG and the SNG. Variates with P-values less than 0.05 were included in the multivariate logistic regression model to explore independent predictors of needle sharing.
All analyses mentioned above were performed using the statistical software Statistical Package for Social Sciences (SPSS, version 22.0, IBM corporation, Armonk, USA) and a P-value of <0.05 was considered statistically significant.
The correlation coefficient between each factor was calculated using SPSS v22.0, and network analyses of the NSNG and the SNG were performed using Cytoscapev 2.8.3 (http://www.cytoscape.org/) [18, 19]. Based on correlation coefficients above 0.1 or below −0.1, and P-values less than 0.05, a network of correlation profiles consisting of the connections between each risk factor were generated, with nodes corresponding to risk factors and edges to pairs of risk factors that were significantly correlated. Network statistics that were used to describe global and local network properties were calculated and edges were treated as undirected as they represent mutual interactions. The calculation of global network properties contained the number of nodes and edges, the average number of neighbours, the network diameter, density, centralization, heterogeneity, clustering coefficient, and the characteristic path length. The calculation of local network properties contained the betweenness centralities (BCs) and closeness centralities (CCs) of each node. Random network statistics were also calculated from 1 000 randomly generated networks that had the same number of nodes and edges as the empirical networks, and were randomly assembled using the Erdös–Rényi model networks. Network module identification was performed using the MCODE v1.3 plugin for Cytoscape .
Sociodemographic characteristics of study participants
Initially, 2 048 drug users were recruited from the Dehong, Baoshan, Lincang, Qujing, and Zhaotong prefectures in Yunnan. Among them, 1 797 were PWID, but only 1 049 successfully completed the questionnaire (see Fig. 1). Most of the participants were male (97.0 %), of the Han ethnicity (81.2 %), aged between 20 and 40 years (88.1 %), and had school experience (93.6 %). Out of all the participants, 18.7 % were HIV seropositive and 15.3 % self-reported to have other infectious diseases, including hepatitis A, B, or C; tuberculosis; and/or malaria (see Table 1).
Among the 1 049 PWID, 37.5 % (n = 393) had a history of needle sharing. Among these, 13.2 % (n = 52) shared needles with those known to be HIV seropositive. No differences were observed between the NSNG and the SNG in terms of age (30.77 ± 7.17 versus 30.95 ± 6.93; P =0.349), gender (P = 0.543), or marital status (P=0.230). Compared to PWID with a needle-sharing history, those without a needle-sharing history were more likely to be of the Han ethnicity, have full-time jobs, and be highly educated (P < 0.01).
Factors associated with needle sharing
Because such a high percentage of needle sharing among PWID was observed, a univariate logistic regression analysis was conducted to compare the different percentage of PWID who partake in needle-sharing behaviours between different substratification of each factor (see Table 1). Results showed that needle sharing was not significantly correlated with age, gender, or marital status of a drug user, however, the PWID who were jobless, belonged to a minor ethnicity, had a low education, and lived in Myanmar-bordering prefectures (Dehong and Lincang) had a greater tendency toward needle-sharing behaviours. A longer duration and greater frequency of drug abuse/injection increased the possibility of sharing needles. Moreover, PWID who combined injecting drugs with oral intake, bought needles from hospital, drugstore or other resources, or reused needles were found to share needles more often. In addition, use of diazepam and heroin facilitated needle sharing, while using other drugs did not. Lastly, PWID who had been diagnosed with sexually transmitted diseases (STDs), HIV, or other infectious diseases (hepatitis, tuberculosis, and/or malaria) were more likely to share needles, suggesting that sharing needles might transmit blood borne diseases.
All 15 factors related to needle sharing were included in the multivariate logistic regression analysis (see Table 1). Of these, living in Zhaotong, Qujing, Dehong or Lincang, using diazepam, injecting for a long time, reusing needles often, and being HIV positive and having other infectious diseases, such as hepatitis, tuberculosis, and/or malaria, were found to be independently associated with needle sharing. These results showed that the factors affecting needle-sharing behaviours were mainly due to drug use behaviours rather than sociodemographic factors, verifying that contracting blood borne diseases (e.g. HIV, viral hepatitis) might be the main consequence of needle sharing.
The global and local properties of correlation networks of PWID in the NSNG and the SNG
Global properties of empirical and random networks of the NSNG and the SNG
Never shared needles
Network clustering coefficient
Network clustering coefficientRandom
Characteristic path length
Characteristic path lengthRandom
Number of nodes
Number of edges
Average number of neighbours
Network clustering coefficient
Characteristic path length
In the correlation network of the SNG, the number of sexual partners was negatively correlated with residence, age, duration of drug abuse/injection, and HIV serostatus, indicating that long-term, frequent drug use might inhibit sexual desire, and that the predominant route of HIV transmission among PWID could be needle sharing (through blood).
The module in the correlation network of the NSNG and the SNG consisted of four and seven hub nodes, respectively, which were involved in many interactions and represented functionally important factors. Changes in hub nodes could affect the most nodes in their vicinity, as well as the whole structure of the correlation network.
In comparison to the NSNG, the network size, average number of neighbours, network density, and network clustering coefficient of the SNG were larger, while the network heterogeneity was lower (see Table 2), indicating that the correlation among each risk factor in the SNG was higher. In addition, as the diameter of the correlation network in both groups was equal, the characteristic path length of the correlation network of the SNG was lower. Changes in each risk factor could permeate the whole correlation network more quickly in the SNG by minimizing the transition time between the interactions of each factor. Paired with higher network centralization, a lower characteristic path length implied that the resilient ability of changes in the network of the SNG was also lower. Overall, these results suggest that the global properties of correlation networks vary between the NSNG and the SNG, and that harm reduction measures aimed at the SNG could achieve goals in a more rapid and durable manner.
BCs and CCs in the NSNG and the SNG
Never shared needles
BC value (rank)
CC value (rank)
BC value (rank)
CC value (rank)
Other drug usea
Method of drug intake
Source of syringe
Duration of drug injection
Drug use frequency (per day)
Times reusing needle
Duration of drug abuse
Most recent drug use
Number of sexual partners
Ever diagnosed with a STD
Other infectious diseasesc
Drug of choice
Key factors for reducing needle-sharing behaviours
Based on the regression analysis, six factors were found to be independently associated with needle sharing. From the network analysis, seven factors were above average for BCs and CCs for all nodes in the SNG and eight in the NSNG. All factors, a total of 12, could be targets for harm reduction measures. As these 12 factor-associated PWID would almost covered the entire population of PWID, it’s difficult to taking measures among this group.
To narrow down the key factors for reducing needle-sharing behaviours among PWID, combined network and regression analyses were conducted. It was elucidated that taking measures aimed at the SNG would achieve goals quicker and more successfully. The BCs and CCs of HIV serostatus and marital status were above average only in the SNG;HIV seropositivity was found to be associated with needle sharing behaviours, while any subcategory of marital status (unmarried, married, cohabitated, divorced, or widowed) was not.
In summary, HIV seropositivity is a unique factor that is associated with needle sharing and has a strong effect only on PWID who partake in needle-sharing behaviours. Therefore, those who are HIV seropositive would be the ideal target group for harm reduction programs.
Research on needle-sharing behaviours has been conducted across different provinces of China, and consequent findings have contributed tremendously to the development of successful harm-reduction programs. From January to April 2001, J. T. Lau et al.  administered a survey of 262 drug users in Shenzhen, Guangdong province, and found that 60.6 % of males and 45.3 % of females among PWID shared needles with others. From August 2003 to June 2004, Dick Chamla et al.  conducted a survey of 266 PWID in Chengdu, Sichuan, and found that 38.72 % had shared needles with their partners. In the present study, we systematically reported the characteristics of needle-sharing behaviours among PWID in the Yunnan province. Although the total prevalence of needle sharing in Yunnan was not significantly different compared to other provinces in China, the percentage of those who shared needles greatly varied across the different prefectures (15.0–72.4 %). We also found that a small quantity of PWID shared needles with partners who are known to be HIV seropositive. These results provide important information on the educational interventions that are immediately needed, and highlight which regions require the most attention.
Factors associated with needle sharing have been described in previous studies. For example, in 2007, Choi et al. reported that sex, ethnicity, having a primary partner who injects drugs, family support, and social discrimination were all significantly associated with needle sharing, based on a cross-sectional survey of 200 PWID in Sichuan, China . Based on our results, six factors were found to be independently associated with needle sharing. The results did not significantly differ from those of previous studies, though larger numbers of participants and factors were taken into consideration.
By comparing the correlation networks of PWID in the NSNG and the SNG, we found that programs that attempt to reform the behaviours of PWID who partake in needle-sharing behaviours might be more successful than for PWID who do not partake in needle-sharing behaviours. Further, by calculating the BCs and CCs for each node of the NSNG and the SNG, we found that HIV serostatus and marital status have a strong effect on PWID in the SNG, making it possible to distinguish PWID who partake in needle-sharing behaviours from those who do not. Finally, by combining network analysis with regression analysis, HIV seropositivity was found to be a unique factor that had a strong effect only on PWID in the SNG, and was significantly associated with needle sharing. All these findings drew a clearer picture for harm reduction programs: taking actions aimed at PWID who are HIV seropositive may work better to reduce needle-sharing behaviours across all groups of PWID.
According to a study by M. Vazirian et al. , shared use of needles or syringes was significantly lower among PWID who received an estimated ≥7 syringes per week than those who did not [2.9 % versus 35.1 %; adjusted odds ratio (OR) = 14.36, 95 % confidence interval (CI) 2.30–89.56]. In China, the prevalence of HIV among PWID was approximately 10.9 % by 2011 (compared to 18.7 % in the present study), meaning that the HIV seropositive group was a relatively smaller group . Therefore, focusing on PWID who are HIV seropositive can ensure that this group receives clean needles, especially in locations that cannot afford to distribute needles to all PWID. Free antiviral drugs are available in healthcare centres in China and HIV seropositive patients can receive them voluntarily, making it more convenient to provide PWID who are HIV seropositive with clean needles.
This study had several limitations. First, a timeframe was not applied to address PWID who most recently partook in drug injection and needle-sharing behaviours, influencing the timeliness of this study. Second, characteristics of harm reduction were not included in the questionnaire. Third, using hepatitis C virus as a predictor of a blood borne disease may have been better than using HIV, as the seroprevalence of hepatitis C is much higher than HIV among PWID. Despite these limitations, this study had evaluated the needle-sharing characteristics among PWID in Yunnan, and could provide useful information for the following harm reduction programs.
In summary, based on a cross-sectional survey, we found that the percentage of those who shared needles among PWID in Yunnan was as high as reported in previous studies in other provinces in China. By combining regression and network analyses, PWID who are HIV seropositive were found to be the ideal target group for taking measures to reduce needle-sharing behaviours among PWID.
BC, betweenness centrality; CC, closeness centrality; HIV-1, human immunodeficiency virus type 1; NSNG, never sharing needles group; PWID, people who inject drugs; SNG, sharing needles group; STD, sexually transmitted disease.
We thank Dr. Muhammad Shahzad for his assistance in editing the paper. We also thank all the participants who took part in this research. This work was supported in part by grants from the National Natural Science Foundation of China (81271892, U1302224, U1202228), and the Key Scientific and Technological Program of China (2012ZX10001-006, 2012ZX10001-007). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the paper.
YTZ, WP, and XC conceived and designed the study. YHZ, FLL, ZHY, LD, and HL collected the samples. XC performed the experiments. XC and LZ analyzed the data and drafted the paper. YTZ, WP, and MY critically revised the paper. All authors read and approved the final paper.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
The protocol of this study was approved by the Ethics Committee of the Kunming Institute of Zoology, Chinese Academy of Sciences (approval number: SWYX-2009021; date: January 7, 2009). All participants gave written informed consent.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Beyrer C, Razak MH, Lisam K, Chen J, Lui W, Yu XF. Overland heroin trafficking routes and HIV-1 spread in south and south-east Asia. AIDS. 2000;14(1):75–83.View ArticlePubMedGoogle Scholar
- Li Y, Uenishi R, Hase S, Liao H, Li XJ, Tsuchiura T, et al. Explosive HIV-1 subtype B’ epidemics in Asia driven by geographic and risk group founder events. Virology. 2010;402(2):223–7.View ArticlePubMedGoogle Scholar
- Li Z, He X, Wang Z, Xing H, Li F, Yang Y, et al. Tracing the origin and history of HIV-1 subtype B’ epidemic by near full-length genome analyses. AIDS. 2012;26(7):877–84.View ArticlePubMedGoogle Scholar
- Tee KK, Pybus OG, Li XJ, Han X, Shang H, Kamarulzaman A, et al. Temporal and spatial dynamics of human immunodeficiency virus type 1 circulating recombinant forms 08_BC and 07_BC in Asia. J Virol. 2008;82(18):9206–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Liao H, Tee KK, Hase S, Uenishi R, Li XJ, Kusagawa S, et al. Phylodynamic analysis of the dissemination of HIV-1 CRF01_AE in Vietnam. Virology. 2009;391(1):51–6.View ArticlePubMedGoogle Scholar
- Liu J, Jia Y, Xu Q, Zheng YT, Zhang C. Phylodynamics of HIV-1 unique recombinant forms in China-Myanmar border: implication for HIV-1 transmission to Myanmar from Dehong, China. Infect Genet Evol. 2012;12(8):1944–8.View ArticlePubMedGoogle Scholar
- Van Griensven F, Pitisuttithum P, Vanichseni S, Wichienkuer P, Tappero JW, Sangkum U, et al. Trends in the injection of midazolam and other drugs and needle sharing among injection drug users enrolled in the AIDSVAX B/E HIV-1 vaccine trial in Bangkok, Thailand. Int J Drug Policy. 2005;16(3):171–5.View ArticleGoogle Scholar
- Strathdee SA, Zafar T, Brahmbhatt H, Baksh A, Hassan SU. Rise in needle sharing among injection drug users in Pakistan during the Afghanistan war. Drug Alcohol Depen. 2003;71(1):17–24.View ArticleGoogle Scholar
- Rafiey H, Narenjiha H, Shirinbayan P, Noori R, Javadipour M, Roshanpajouh M, et al. Needle and syringe sharing among Iranian drug injectors. Harm Reduct J. 2009;6(21.10):1186. doi:10.1186/1477-7517-6-21.Google Scholar
- Chakrapani V, Newman PA, Shunmugam M, Dubrow R. Social-structural contexts of needle and syringe sharing behaviours of HIV-positive injecting drug users in Manipur, India: a mixed methods investigation. Harm Reduct J. 2011; 8(9). doi: 10.1186/1477-7517-8-9.
- Zhang L, Li J, Lai W, Feng L, Zeng Y, Liu L, et al. Prevalence and correlates of needle-sharing among new and long-term injection drug users in southwest China. Subst Use Misuse. 2010;45(14):2503–23.View ArticlePubMedGoogle Scholar
- Hope VD, Rogers PA, Jordan L, Paine TC, Barnett S, Parry JV, et al. Sustained increase in the sharing of needles and syringes among drug users in England and Wales. AIDS. 2002;16(18):2494–6.View ArticlePubMedGoogle Scholar
- Matheï C, Shkedy Z, Denis B, Kabali C, Aerts M, Molenberghs G, et al. Evidence for a substantial role of sharing of injecting paraphernalia other than syringes/needles to the spread of hepatitis C among injecting drug users. J Viral Hepatitis. 2006;13(8):560–70.View ArticleGoogle Scholar
- Luke DA, Harris JK. Network analysis in public health: history, methods, and applications. Annu Rev Public Health. 2007;28:69–93.View ArticlePubMedGoogle Scholar
- Hafner-Burton EM, Kahler M, Montgomery AH. Network analysis for international relations. Int Organ. 2009;63(03):559–92.View ArticleGoogle Scholar
- Borgatti SP, Mehra A, Brass DJ, Labianca G. Network analysis in the social sciences. Science. 2009;323(5916):892–5.View ArticlePubMedGoogle Scholar
- Böde C, Kovács IA, Szalay MS, Palotai R, Korcsmáros T, Csermely P. Network analysis of protein dynamics. FEBS Lett. 2007;581(15):2776–82.View ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–2.View ArticlePubMedGoogle Scholar
- Parkes E, Duffy B. An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations. Comput Phys Commun. 1996;98(3):288–300.View ArticleGoogle Scholar
- Lau JT, Feng T, Lin X, Wang Q, Tsui H. Needle sharing and sex-related risk behaviours among drug users in Shenzhen, a city in Guangdong, southern China. AIDS Care. 2005;17(2):166–81.View ArticlePubMedGoogle Scholar
- Chamla D, Chamla JH, Dabin W, Delin H, Rennes N. Transition to injection and sharing of needles/syringes: Potential for HIV transmission among heroin users in Chengdu, China. Addict Behav. 2006;31(4):697–701.View ArticlePubMedGoogle Scholar
- Choi SY, Cheung Y, Jiang Z. Ethnicity and risk factors in needle sharing among intravenous drug users in Sichuan Province, China. AIDS Care. 2007;19(1):1–8.View ArticlePubMedGoogle Scholar
- Vazirian M, Nassirimanesh B, Zamani S, Ono-Kihara M, Kihara M, Ravari SM, et al. Needle and syringe sharing practices of injecting drug users participating in an outreach HIV prevention program in Tehran, Iran: a cross-sectional study. Harm Reduct J. 2005;2(1):19. doi:10.1186/1477-7517-2-19.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Guo W, Li D, Ding Z, McGoogan JM, Wang N, et al. HIV epidemic among drug users in China: 1995–2011. Addiction. 2015;110(S1):20–8.View ArticlePubMedPubMed CentralGoogle Scholar