PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of pubhealthrepPublic Health Reports
 
Public Health Rep. 2016 May-Jun; 131(3): 491–499.
PMCID: PMC4869089

The Impact of Benzodiazepine Use on Mortality Among Polysubstance Users in Vancouver, Canada

Geoffrey R.T. Walton, MD,a,b Kanna Hayashi, PhD,a,b Paxton Bach, MD,b Huiru Dong, MSc,a Thomas Kerr, PhD,a Keith Ahamad, MD,a MJ Milloy, PhD,a Julio Montaner, MD,a and Evan Wood, MD, PhDcorresponding authora,b

Abstract

Objective

Illicit drug use is a well-established risk factor for increased morbidity and mortality. However, little is known about the impact of benzodiazepine use on mortality among populations of polysubstance users. This study aimed to identify the effect of benzodiazepine use on mortality among a prospective cohort of people in Canada who inject drugs (PWID).

Methods

A cohort of PWID in Vancouver, Canada, was prospectively followed from May 1996 through November 2013. Data on participants were linked to the provincial vital statistics registry to ascertain mortality rates and causes of death. Multivariable extended Cox regression with time-dependent variables was used to investigate the relationship between benzodiazepine use and time to all-cause mortality.

Results

During the study period, 2,802 participants were followed for a median of 67 months (interquartile range: 25–107). In total, 527 (18.8%) participants died, for an incidence density of mortality of 2.9 (95% confidence interval [CI] 2.7, 3.2) deaths per 100 person-years. After adjusting for HIV infection and other potential confounders, benzodiazepine use was independently associated with increased all-cause mortality (adjusted hazard ratio = 1.86, 95% CI 1.38, 2.51) and had a higher risk for mortality than all other traditional substances of abuse among this population. Results were consistent when non-overdose mortality was considered.

Conclusion

In this setting, benzodiazepine use was more strongly associated with mortality than any other substance of abuse. Greater recognition of the safety concerns related to benzodiazepines and strategies to prevent diversion to illicit use are needed.

Substance abuse is a growing burden on health worldwide, even as its prevalence and contribution to morbidity and mortality is increasingly recognized and studied.1 However, an accurate estimate of the prevalence and harms of patterns of drug use is difficult to ascertain. The problem is widespread, with recent United Nations estimates of 15–39 million problem drug users and 11–21 million intravenous drug users globally.2 A 2013 review estimated a pooled crude mortality ratio associated with intravenous drug use of 2.35 per 100 person-years globally and 2.64 per 100 person-years in North America.3 Although illicit drug use is associated with many preventable adverse health effects, including premature mortality, most research has focused on traditional illegal drugs of abuse, such as heroin, cocaine, and amphetamine-type stimulants.46

In recent years, however, concern about prescription drug abuse has increased. Studies examining data from large U.S. surveys have described the burden of disease from illicit drug use, including prescription drugs.7,8 These data demonstrated that, in 2010, almost 3% of the U.S. household population reported nonmedical prescription drug use in the previous month, and that more than 15% of 12th graders reported misusing prescription drugs in the previous year.7,8 Canadian studies published in 2009 and 2012 highlighted emerging concerns about the harms of prescription opioid misuse, including an epidemic of prescription opioid-related overdose deaths.9,10 Nonmedical prescription opioid use has surpassed heroin use among street drug users in some settings.9

Benzodiazepines have been used since chlordiazepoxide was approved in 1960.11 Side effects of benzodiazepine use include psychomotor retardation, memory impairment, and development of physiologic dependence.11,12 In Canada, concern is rising about the diversion and misuse of sedatives, including benzodiazepines.13 Sedatives and hypnotic drugs are among the most common causes of visits to the emergency department for adverse drug events.14 An examination of Canadian data to determine the relationship between benzodiazepines and opioid use found that most opioid-related deaths were associated with non-opioid central nervous depressants, including benzodiazepines.10 Benzodiazepines were also recently associated with opioid-related death during methadone therapy in Canada.15 Studies on benzodiazepine use among people who inject drugs (PWID) in Canada demonstrated that approximately 35% of PWID reported benzodiazepine use in the previous six months16 and that benzodiazepine use was associated with an increased risk of self-reported nonfatal overdose.17 We examined the prevalence of benzodiazepine use among a prospective cohort of PWID in Vancouver, British Columbia, Canada, to determine its association with mortality.

METHODS

The Vancouver Injection Drug Users Study (VIDUS) and AIDS Care Cohort to Evaluate Access to Survival Services (ACCESS) are open prospective cohorts of drug users in Vancouver. The recruitment and follow-up procedures for the two studies are identical to allow for analyses of merged data; the only substantive differences are (1) human immunodeficiency virus (HIV)-positive participants are followed in ACCESS, whereas HIV-negative participants are followed in VIDUS and (2) ACCESS includes non-injecting (e.g., crack cocaine) drug users, whereas VIDUS only enrolls people who inject drugs in the month prior to enrollment. In both studies, the primary modes of enrollment are self-referral, word of mouth, and street outreach. Detailed sampling and recruitment procedures for these two cohorts are described elsewhere.18,19 To be eligible, participants must be $18 years of age, reside in the greater Vancouver region, and have injected drugs in the month prior to enrollment. All participants provide written informed consent. Participants are given a stipend ($30 Canadian) at each study visit for their time and transportation.

At baseline and semianually thereafter, participants complete an interviewer-administered questionnaire on a range of topics, including demographic characteristics, injection and non-injection drug use, and sexual risk behaviors. Venous blood samples drawn at each visit are used for serologic analyses, including testing for HIV and hepatitis C virus (HCV) antibodies. All participants have private interviews and are offered pretest and posttest counseling with trained nurses. Referral for free HIV/acquired immunodeficiency syndrome care is provided to those found to be HIV positive, and these participants are subsequently followed in ACCESS.

Our study included individuals who were recruited from May 1996 through November 2013 and who reported a history of injection drug use at baseline. Fifty-two of 2,854 (1.8%) participants who did not report a history of injection drug use at baseline were excluded; these participants were part of the ACCESS cohort. We collected mortality data, including data on underlying cause of death among participants, through a confidential record linkage with the British Columbia Vital Statistics Agency and ongoing follow-up with contacts (e.g., friends) provided by participants. The Vital Statistics database records cause of death -according to the International Classification of Diseases, 10th Revision.

The primary endpoint in this analysis was all-cause mortality, with overdose and non-overdose mortality considered in subanalyses. The primary explanatory variable was any benzodiazepine use in the previous six months. We also considered demographic, behavioral, clinical, and other exposure characteristics as potential confounders of the association between benzodiazepine use and mortality, based on prior studies of health and social harms among PWID in our setting.2022 These confounders were substance-using behaviors in the previous six months, including at least daily cocaine injection, at least daily heroin injection, at least daily amphetamine injection, at least daily crack cocaine smoking, at least daily speedball (cocaine and heroin) injection, and at least daily alcohol use. Other variables were age, sex (male and female), time since first injection, race (white and nonwhite), HIV serostatus (positive and negative), HCV antibody serostatus (positive and negative), methadone maintenance therapy, sex work, and living in unstable housing in the previous six months. We defined unstable housing as living in one of the city's Downtown Eastside single-room-occupancy hotels, shelters, or other transitional housing, or living on the street. All variables were identical to prior reports on injection drug users in Vancouver.23

As an initial analysis, we calculated the frequency and prevalence of benzodiazepine use at baseline, stratified by each binary explanatory variable. We compared the prevalences using the Wald c2 test. For continuous explanatory variables, we calculated the median and interquartile range (IQR) and compared them using the Wilcoxon rank sum test. We then calculated all-cause mortality rates and overdose mortality rates, including 95% confidence intervals (CIs), using the Poisson distribution.

Next, we used extended Cox regression to examine the bivariable relationship between each explanatory variable and time to all-cause mortality.24 We treated all behavioral variables as time-varying variables. To fit the multivariable model, we used a conservative stepwise backward selection approach. We included all variables found to be significantly associated with time to all-cause mortality in bivariable analyses at p<0.10 in a multivariable model, and we used a stepwise approach to fit a series of reduced models. After comparing the value of the coefficient associated with benzodiazepine use in the full model with the value of the coefficient in each of the reduced models, we dropped the secondary variable associated with the smallest relative change. We continued this iterative process until the minimum change exceeded 5%. We considered remaining variables as potential confounders in a final multivariable model. We used this procedure previously to estimate the effect of a primary explanatory variable on an outcome of interest after adjusting for secondary variables.25

We also conducted subanalyses in which we restricted the dependent variable to deaths resulting from overdose and non-overdose to examine whether benzodiazepine use was independently associated with overdose deaths or non-overdose deaths or both. We performed all statistical analyses using SAS® version 9.3.26 All p-values were two-sided.

RESULTS

A total of 2,802 individuals were eligible for the present study and were followed for a median of 67.0 months (IQR: 24.8–107.0). Among the study sample, 1,855 (66.2%) participants were men, 1,699 (60.6%) participants were white, and 862 (30.8%) participants were HIV positive at baseline. The median age was 37.1 years (IQR: 29.2–43.8), and the median time since first injection was 14.0 years (IQR: 5.8–23.6). At baseline, 734 (26.2%) participants reported benzodiazepine use in the previous six months, 1,028 (36.7%) injected heroin daily, 802 (28.6%) injected cocaine daily, 683 (24.4%) smoked crack cocaine daily, 524 (18.7%) consumed alcohol daily, and 80 (2.9%) injected amphetamines daily (Table 1).

Table 1.
Characteristics of people who inject drugs enrolled in the VIDUS and ACCESS cohorts, stratified by any benzodiazepine use in the past 6 months, Vancouver, Canada, November 1996 to May 2013

Benzodiazepine use was significantly more prevalent among participants who reported daily heroin injection (prevalence ratio [PR] = 1.36, 95% CI 1.20, 1.54; p<0.001) and daily cocaine injection (PR=1.94, 95% CI 1.72, 2.19, p<0.001) than among those who did not report daily heroin or cocaine injection (Table 1).

During the study period, 527 (18.8%) participants died, yielding an incidence density of 2.9 (95% CI 2.7, 3.2) deaths per 100 person-years. In the bivariable analyses, benzodiazepine use was significantly and positively associated with time to all-cause mortality; the unadjusted hazard ratio (HR) was 1.74 (95% CI 1.29, 2.33; p<0.001). Daily cocaine injection (HR=1.46, 95% CI 1.17, 1.82; p<0.001) and daily alcohol use (HR=1.27, 95% CI 1.00, 1.61; p=0.047) were also significantly and positively associated with all-cause mortality. In the multivariable analysis, after adjusting for potential confounders, including HIV serostatus and daily cocaine injection, benzodiazepine use remained independently and positively associated with time to all-cause mortality (adjusted hazard ratio [AHR] = 1.86, 95% CI 1.38, 2.51; p<0.001) (Table 2).

Table 2.
Bivariable and multivariable Cox regression analysis of factors associated with all-cause mortality (n=527) among people who inject drugs enrolled in the VIDUS and ACCESS cohorts (n=2,802), Vancouver, Canada, November 1996 to May 2013

The subanalysis on non-overdose mortality showed that benzodiazepine use had an HR of 1.87 (95% CI 1.34, 2.62) and an AHR of 2.04 (95% CI 1.45, 2.87) (both p<0.001). The other substance use variable significantly associated with non-overdose mortality was daily cocaine use (HR=1.56, 95% CI 1.21, 2.00; p<0.001) (Table 3). The subanalysis on overdose mortality showed that benzodiazepine use was not significantly associated with overdose mortality (HR=1.40; 95% CI 0.76, 2.57; p=0.284; AHR=1.48, 95% CI 0.80, 2.72; p=0.210). No other substance examined was significantly associated with overdose mortality (Table 4).

Table 3.
Bivariable and multivariable Cox regression analyses of factors associated with non-overdose mortality (n=403) among people who inject drugs enrolled in the VIDUS and ACCESS cohorts (n=2,802), Vancouver, Canada, November 1996 to May 2013
Table 4.
Bivariable and multivariable Cox regression analyses of factors associated with overdose mortality (n=124) among people who inject drugs enrolled in the VIDUS and ACCESS cohorts (n=2,802), Vancouver, Canada, November 1996 to May 2013

Among participants who reported benzodiazepine use during the study period, 70 of 294 deaths (23.8%) were overdoses, whereas among participants who did not report benzodiazepine use, 54 of 233 deaths (23.2%) were overdoses. HIV-related causes accounted for 70 of 294 deaths (23.8%) among those who reported benzodiazepine use and 54 of 233 deaths (23.2%) among those who did not report benzodiazepine use. Among the 294 deaths of participants who reported benzodiazepine use and the 233 deaths of participants who did not report benzodiazepine use, liver-related causes accounted for 18 (6.1%) and 14 (6.0%) deaths, suicide accounted for 16 (5.4%) and two (0.9%) deaths, homicide accounted for eight (2.7%) and nine (3.9%) deaths, accidental causes accounted for 13 (4.4%) and 13 (5.6%) deaths, other causes accounted for 77 (26.2%) and 45 (19.3%) deaths, and unspecificied causes accounted for 22 (7.5%) and 42 (18.0%) deaths, respectively.

DISCUSSION

A substantial proportion of PWID in Vancouver reported some benzodiazepine use in the previous six months. We found a strong independent association between benzodiazepine use and mortality in both an analysis on all-cause mortality and a subanalysis on non-overdose mortality. In our cohort of polysubstance-using PWID, benzodiazepine use was more strongly associated with mortality than any other drug of abuse considered.

Previous studies on prescription drug misuse largely focused on opioids, while the nonmedical use of benzodiazepines received relatively little attention.79,13 Examining substance-using (predominantly opioid-using) populations, reviews of coroners' cases27 and other sources9,28 showed an association between benzodiazepine use and overdose mortality. One Canadian study found that 92% of opioid-related deaths were associated with the concomitant use of benzodiazepines or alcohol.9 The evidence associating benzodiazepines with mortality in the polysubstance-using population is limited to a small number of studies focusing on mostly opioid-using populations, and the association was made only with overdose deaths.9,2729 The participants in our cohort reported a broad range of primary drugs of abuse, with only 36% reporting daily heroin use at baseline. Yet, we still showed an association between benzodiazepine use and mortality.

Benzodiazepines are often abused concomitantly with other substances because they are perceived to heighten the effects of co-administered drugs such as opioids.30 A 2012 study examining the pharmacodynamics of co-administering diazepam with either methadone or buprenorphine found increases in the intensity in subjective drug effects and decreases in the performance on psychological testing.19 It has been hypothesized that during concomitant use of benzodiazepines and injection drugs, benzodiazepine-induced psychomotor symptoms contribute to harms involved in the preparation and injection of drugs that could lead to infection and overdose.11 Concerns about these pharmacodynamic interactions that potentially lead to risky behaviors are well founded. A 2012 study on a Thai population showed that daily users of intravenous midazolam reported higher rates of polysubstance use, higher rates of heroin use, higher rates of overdose, and higher rates of binge drug use compared with people who did not use intravenous midazolam daily.31 Benzodiazepines may also lessen the comedown from stimulants such as cocaine or amphetamines.32 The cognitive and psychomotor effects of benzodiazepine use could contribute to the mortality shown in other studies associated with benzodiazepine use in polysubstance users, particularly because the effects produce measurable differences in risk-taking substance-abuse behaviors.

Our study did not show a significant association between benzodiazepine use and overdose mortality, although it did indicate a nonsignificant trend toward an association. The reason for the lack of association is unclear, but it may be that participants in our study used fewer opioids than participants in previous studies. It is also possible that rates of death caused by overdose in our cohort were lower than rates in previous studies because participants in our study had greater access to safe places to inject drugs, such as the supervised injection facility Insite. Another possible reason is that our sample size was insufficient to detect a significant effect. Additionally, as with any mortality data, potentially overdose deaths may have been misclassified as non-overdose deaths. The association between benzodiazepine use and non-overdose mortality could be explained by outcomes that were not measured in this study; for example, as a result of falls, motor vehicle accidents, and other factors that may be consequences of benzodiazepine use, all of which could increase the risk of deaths that would not be classified as overdose related.

Aside from the well-documented psychomotor and addiction liability caused by benzodiazepine use,11,32,33 other concerns about the safety and appropriateness of benzodiazepines exist for many common clinical indications. Evidence points to the association between benzodiazepine use and long-term cognitive effects in multiple cognitive domains,34 an increased risk of falls and fractures in the elderly, and motor vehicle accidents.29 In addition to the dangers of benzodiazepine use, research shows that the efficacy of benzodiazepines for common clinical indications (e.g., anxiety, insomnia) is limited, especially for chronic use.35 Our data suggest clinicians should seek alternative medications when patients, especially polysubstance-using patients, seek benzodiazepines.

Limitations

This study had several limitations. First, this study was limited by its observational design, although it would not be ethical or possible to study this population in a controlled fashion. Second, the results of this study may not be generalizable to other PWID populations because of the unique environment of our setting and the unique behaviors of our study participants. Nevertheless, the diversity of drugs used by participants in our study was greater than the diversity of drugs used by participants in previous studies.28 Third, we did not routinely assess mental illness in this cohort, although individuals with severe mental illness (e.g., psychosis) were not eligible to participate if they could not provide informed consent. Fourth, because of self-reported drug use, our results may have been affected by socially desirable responses, recall bias, or other reporting biases; however, previous studies have shown self-reporting of drug-use behavior to be valid.36 Finally, our primary endpoint was mortality, an objective measure obtained through the provincial database. Our results did not account for people who died in other provinces; however, it has been shown that migration from the province is low in this population.37

CONCLUSION

We found a strong independent association between benzodiazepine use and all-cause mortality in a sample of PWID in Vancouver, demonstrating a greater risk for mortality than all other drugs of abuse considered. This association was robust and persisted when we investigated the subgroup of non-overdose mortality. Our results add to the literature identifying the role of benzodiazepine use in overdose mortality among opioid users. This study was unique in that it was a prospective analysis of polysubstance users, was adjusted for other drugs of abuse as well as HIV-seropositive status, and considered the effect of benzodiazepine use on overdose and non-overdose mortality. Future studies should further examine mortality risk of more specific drug-use patterns in the PWID population, including potentially the co-usage of prescription drugs with other traditional substances of abuse. Given the known harms of benzodiazepines,9,11,12,2733 these findings highlight the importance of physician education aimed at promoting alternative medications to benzodiazepines and safer prescription strategies to avoid diversion to the street.

Footnotes

The authors thank the study participants for their contributions to the research, as well as current and past researchers and staff members. This study was supported by the U.S. National Institutes of Health (NIH) (Vancouver Injection Drug Users Study: U01DA038886; AIDS Care Cohort to Evaluate Access to Survival Services: R01DA021525) and the Canadian Institutes of Health Research through the Canadian Research Initiative on Substance Misuse (SMN-139148). Dr. Milloy is supported in part by NIH (R01-DA051525). Dr. Hayashi is supported by the Canadian Institutes of Health Research Fellowship. This research was undertaken, in part, with funding from a Tier 1 Canada Research Chair in Inner City Medicine, which supports Dr. Wood. Dr. Montaner is supported with grants paid to his institution by the British Columbia Ministry of Health and by NIH (R01DA036307). He has also received limited unrestricted funding, paid to his institution, from Abbvie, Bristol-Myers Squibb, Gilead Sciences, Janssen, Merck, and ViiV Healthcare. Dr. Ahamad is supported by the National Institute on Drug Abuse (R25DA037756).

The viewpoints in this article are those of the authors and do not represent those of the funding agencies. The two cohort studies have been approved by the University of British Columbia/Providence Healthcare Research Ethics Board.

REFERENCES

1. Whiteford HA, Degenhardt L, Rehm J, Baxter AJ, Ferrari AJ, Erskine HE, et al. Global burden of disease attributable to mental and substance use disorders: findings from the Global Burden of Disease Study 2010. Lancet. 2013;382:1575–86. [PubMed]
2. Degenhardt L, Hall W. Extent of illicit drug use and dependence, and their contribution to the global burden of disease. Lancet. 2012;379:55–70. [PubMed]
3. Mathers BM, Degenhardt L, Bucello C, Lemon J, Wiessing L, Hickman M. Mortality among people who inject drugs: a systematic review and meta-analysis. Bull World Health Organ. 2013;91:102–23. [PubMed]
4. Hulse GK, English DR, Milne E, Holman CD. The quantification of mortality resulting from the regular use of illicit opiates. Addiction. 1999;94:221–9. [PubMed]
5. Degenhardt L, Singleton J, Calabria B, McLaren J, Kerr T, Mehta S, et al. Mortality among cocaine users: a systematic review of cohort studies. Drug Alcohol Depend. 2011;113:88–95. [PubMed]
6. Singleton J, Degenhardt L, Hall W, Zabransky T. Mortality among amphetamine users: a systematic review of cohort studies. Drug Alcohol Depend. 2009;105:1–8. [PubMed]
7. Johnston LD, O'Malley PM, Bachman JG, Schulenberg JE, Miech RA. Ann Arbor (MI): University of Michigan Institute for Social Research; 2014. Monitoring the future: national survey results on drug use, 1975–2013: volume I, secondary school students.
8. Schulden JD, Lopez MF, Compton WM. Clinical implications of drug abuse epidemiology. Psychiatr Clin North Am. 2012;35:411–23. [PMC free article] [PubMed]
9. Fischer B, Argento E. Prescription opioid related misuse, harms, diversion and interventions in Canada: a review. Pain Physician. 2012;15(3 Suppl):ES191–203. [PubMed]
10. Dhalla IA, Mamdani MM, Sivilotti ML, Kopp A, Qureshi O, Juurlink DN. Prescribing of opioid analgesics and related mortality before and after the introduction of long-acting oxycodone. CMAJ. 2009;181:891–6. [PMC free article] [PubMed]
11. Lader M. Benzodiazepines revisited—will we ever learn? Addiction. 2011;106:2086–109. [PubMed]
12. Longo LP, Johnson B. Addiction: part I. Benzodiazepines—side effects, abuse risk and alternatives. Am Fam Physician. 2000;61:2121–8. [PubMed]
13. National Advisory Council of Prescription Drug Misuse. Ottawa (ON): National Advisory Council of Prescription Drug Misuse, Canadian Centre on Substance Abuse; 2013. First do no harm: responding to Canada's prescription drug crisis.
14. Hampton LM, Daubresse M, Chang HY, Alexander GC, Budnitz DS. Emergency department visits by adults for psychiatric medication adverse events. JAMA Psychiatry. 2014;71:1006–14. [PMC free article] [PubMed]
15. Leece P, Cavacuiti C, Macdonald EM, Gomes T, Kahan M, Srivastava A, et al. Predictors of opioid-related death during methadone therapy. J Subst Abuse Treat. 2015;57:30–5. [PubMed]
16. Ottawa (ON): Public Health Agency of Canada, Center for Communicable Disease and Infection Control, Infectious Disease Prevention and Control Branch; 2013. I-track: enhanced surveillance of HIV, hepatitis C and associated risk behaviours among people who inject drugs in Canada: phase 2 report.
17. Fischer B, Brissette S, Brochu S, Bruneau J, el-Guebaly N, Noel L, et al. Determinants of overdose incidents among illicit opioid users in 5 Canadian cities. CMAJ. 2004;171:235–9. [PMC free article] [PubMed]
18. Tyndall MW, Currie S, Spittal P, Li K, Wood E, O'Shaughnessy MV, et al. Intensive injection cocaine use as the primary risk factor in the Vancouver HIV-1 epidemic. AIDS. 2003;17:887–93. [PubMed]
19. Milloy MJ, Kerr T, Buxton J, Rhodes T, Krusi A, Guillemi S, et al. Social and environmental predictors of plasma HIV RNA rebound among injection drug users treated with antiretroviral therapy. J Acquir Immune Defic Syndr. 2012;59:393–9. [PMC free article] [PubMed]
20. Hayden A, Hayashi K, Dong H, Milloy MJ, Kerr T, Montaner JS, et al. The impact of drug use patterns on mortality among polysubstance users in a Canadian setting: a prospective cohort study. BMC Public Health. 2014;14:1153. [PMC free article] [PubMed]
21. Nolan S, Hayashi K, Milloy MJ, Kerr T, Dong H, Lima VD, et al. The impact of low-threshold methadone maintenance treatment on mortality in a Canadian setting. Drug Alcohol Depend. 2015;156:57–61. [PMC free article] [PubMed]
22. Zivanovic R, Milloy MJ, Hayashi K, Dong H, Sutherland C, Kerr T, et al. Impact of unstable housing on all-cause mortality among persons who inject drugs. BMC Public Health. 2015;15:106. [PMC free article] [PubMed]
23. Spittal PM, Craib KJ, Wood E, Laliberte N, Li K, Tyndall MW, et al. Risk factors for elevated HIV incidence rates among female injection drug users in Vancouver. CMAJ. 2002;166:894–9. [PMC free article] [PubMed]
24. Kleinbaum DG, Klein M. 3rd ed. New York: Springer-Verlag; 1996. Survival analysis: a self-learning text.
25. Milloy MJ, Kerr T, Buxton J, Rhodes T, Guillemi S, Hogg R, et al. Dose-response effect of incarceration events on nonadherence to HIV antiretroviral therapy among injection drug users. J Infect Dis. 2011;203:1215–21. [PMC free article] [PubMed]
26. SAS Institute, Inc. Cary (NC): SAS Institute, Inc.; 2012. SAS®: Version 9.3.
27. Zamparutti G, Schifano F, Corkery JM, Oyefeso A, Ghodse AH. Deaths of opiate/opioid misusers involving dihydrocodeine, UK, 1997–2007. Br J Clin Pharmacol. 2011;72:330–7. [PMC free article] [PubMed]
28. Gossop M, Stewart D, Treacy S, Marsden J. A prospective study of mortality among drug misusers during a 4-year period after seeking treatment. Addiction. 2002;97:39–47. [PubMed]
29. Charlson F, Degenhardt L, McLaren J, Hall W, Lynskey M. A systematic review of research examining benzodiazepine-related mortality. Pharmacoepidemiol Drug Saf. 2009;18:93–103. [PubMed]
30. Lintzeris N, Mitchell TB, Bond AJ, Nestor L, Strang J. Pharmacodynamics of diazepam co-administered with methadone or buprenorphine under high dose conditions in opioid dependent patients. Drug Alcohol Depend. 2007;91:187–94. [PubMed]
31. Hayashi K, Suwannawong P, Ti L, Kaplan K, Wood E, Kerr T. High rates of midazolam injection and associated harms in Bangkok, Thailand. Addiction. 2013;108:944–52. [PubMed]
32. Lader MH. Limitations on the use of benzodiazepines in anxiety and insomnia: are they justified? Eur Neuropsychopharmacol. 1999;9(Suppl 6):S399–405. [PubMed]
33. Licata SC, Rowlett JK. Abuse and dependence liability of benzodiazepine-type drugs: GABA(A) receptor modulation and beyond. Pharmacol Biochem Behav. 2008;90:74–89. [PMC free article] [PubMed]
34. Barker MJ, Greenwood KM, Jackson M, Crowe SF. Persistence of cognitive effects after withdrawal from long-term benzodiazepine use: a meta-analysis. Arch Clin Neuropsychol. 2004;19:437–54. [PubMed]
35. Dell'osso B, Lader M. Do benzodiazepines still deserve a major role in the treatment of psychiatric disorders? A critical reappraisal. Eur Psychiatry. 2013;28:7–20. [PubMed]
36. Darke S. Self-report among injecting drug users: a review. Drug Alcohol Depend. 1998;51:253–63. discussion 67-8. [PubMed]
37. Rachlis BS, Hogg RS, Wood E, Li K, Kerr T. Factors associated with geographic migration among a cohort of injection drug users. Health Place. 2008;14:536–43. [PubMed]

Articles from Public Health Reports are provided here courtesy of SAGE Publications