Importance
The unprecedented increase in unintentional overdose events that has occurred in tandem with escalating sales of prescription opioids over the past 2 decades has raised concerns about whether the therapeutic use of opioids has contributed to increases in overdose injury. Few controlled studies have examined the extent to which ecologic measures of increases in opioid prescribing and overdose injuries reflect risk among patients prescribed opioids, let alone whether some opioid regimens are safer than others.
Objective
To examine whether the risk of unintentional overdose injury is associated with the duration of opioid action (ie, long-acting vs short-acting formulations).
Design, Setting, and Participants
A propensity score–adjusted cohort study was conducted using population-based health care utilization data from the Veterans Administration Healthcare System. The patients were veterans with chronic painful conditions who began therapy with opioid analgesics between January 1, 2000, and December 31, 2009.
Main Outcomes and Measures
Unintentional overdoses that are explicitly coded using International Classification of Disease, Ninth Revision codes as drug or medication poisonings of accidental intent (E850.x-860.x) or undetermined intent (E980.x or drug poisoning [960.x-980.x] without an accompanying external cause of injury code).
Results
A total of 319 unintentional overdose events were observed. Patients initiating therapy with long-acting opioids were more than twice as likely to overdose compared with persons initiating therapy with short-acting opioids. After adjustment for age, sex, opioid dose, and other clinical characteristics, patients receiving long-acting opioids had a significantly higher rate of overdose injury than did those receiving short-acting opioids (hazard ratio [HR], 2.33; 95% CI, 1.26-4.32). The risk associated with long-acting agents was particularly marked during the first 2 weeks after initiation of treatment (HR, 5.25; 1.88-14.72).
Conclusions and Relevance
To our knowledge, the findings of the present study provide the first evidence that the risk of unintentional overdose injury is related to the prescribed opioid’s duration of action. If replicated in other cohorts, our findings suggest that clinicians weighing the benefits and risks of initiating different opioid regimens should consider not only the daily dose prescribed but also the duration of opioid action, favoring short-acting agents whenever possible, especially during the first 2 weeks of therapy.
The unprecedented increase in unintentional overdose events, both fatal and nonfatal, that has occurred in tandem with escalating sales of prescription opioids over the past 2 decades has raised concerns about whether the therapeutic use of opioids to treat moderate to severe pain has contributed to increases in overdose injury.1-6 However, few controlled studies have examined the extent to which observations from aggregate-level data about parallel increases in opioid prescribing and overdose injuries reflect risk among patients receiving opioids,7-11 let alone whether some opioid regimens are safer than others.
The relationship between particular opioid regimens and the risk of overdose remains poorly understood for several reasons. First, randomized clinical trials have been too small to produce meaningful risk estimates for overdose events (either fatal or nonfatal). Although unintentional drug overdose claimed approximately 33 000 lives in 201012 (nearly the number of individuals who died in motor vehicle crashes) and an estimated 10 times as many emergency department visits each year,13 event rates are still too rare to assess with randomized trials. Second, most nonrandomized (ie, observational) studies1-6 that have examined this issue have relied on ecologic data and thus are subject to the ecologic fallacy (ie, drawing inferences about individuals based on analyses of group data). Third, the few individual-level observational studies that have most thoughtfully investigated the relationship between opioid prescription regimens and unintentional overdose risk have focused on opioid dose8-11 rather than more readily modified attributes of prescriptions, such as whether the particular opioid agent or type of opioid chosen matters when the dose is controlled. Drawing actionable inferences about preventing unintentional overdose from existing studies is further compromised by limitations in design and analysis, including failure to distinguish between intentional and unintentional overdose outcomes7,9,10,14 (and other sources of measurement error, such as ascertaining potential confounders after initiation of opioid therapy7,9,10,15,16) and selection bias, especially because of the inclusion of prevalent users in the analytic sample.7-11
The present cohort study used clinical and pharmacy data from the Veterans Administration (VA) Healthcare System and appropriate pharmacoepidemiologic techniques to examine whether the risk of unintentional overdose injury is associated with the duration of opioid action (ie, long-acting vs short-acting opioid formulations) and, if so, whether the overdose risk is modified by the duration of continuous opioid use. To our knowledge, the present study is the first to address these questions using an incident-user design (ie, a study that sets the cohort’s inception date according to patients’ new use of opioids, in contrast with studies that follow a prevalent user design and enroll patients who were already receiving opioids when follow-up began).
Study Design and Data Collection
The present cohort study examined the relationship between opioid prescription use and nonfatal overdose among veterans receiving care within the US Department of Veterans Affairs Healthcare System between January 1, 2000, and December 31, 2009. All analyses were conducted at the Massachusetts Veterans Epidemiology Research and Information Center, VA Cooperative Studies Coordinating Center, VA Boston Healthcare System. The institutional review board of the VA Boston Healthcare System approved the study. Since health care utilization data stripped of identifiers were used, the need for informed consent was waived.
National pharmacy and administrative data from the Veterans Health Administration (VHA) were linked. Data from the National Patient Care Database includes patient demographics, the location and clinic where services were provided, International Classification of Diseases, Ninth Revision (ICD-9) diagnosis codes, Current Procedural Terminology codes (for outpatient procedures), and ICD-9 procedure codes (for inpatient procedures) for all years since 1997. We obtained deidentified outpatient prescription medication data for opioid agents and concomitant medications from the Pharmacy Benefits Management Database. Pharmacy data include all prescriptions, dosages, days’ supply, and quantities dispensed as well as the date dispensed, in a VA facility and outpatient prescription orders filled at a VA pharmacy or consolidated mail outpatient pharmacy.
We identified all patients who filled an opioid analgesic prescription between January 1, 2000, and December 31, 2009 (Figure). The index date was defined as the first day on which an opioid prescription was filled during this period. Our sample was restricted to new users of opioids, defined as persons having a minimum of 6 months without the use of an opioid before the index opioid prescription (ie, we required a 6-month washout period). We further limited the sample to new users with chronic noncancer pain as the likely indication for treatment prior to or including the index date. Chronic pain was characterized using the definition of Bohnert et al11 in their study of fatal opioid overdose among VHA beneficiaries, including ICD-9 diagnostic codes for general chronic pain; headache (other headache syndromes; migraines; tension, not otherwise specified; and atypical face pain); back and neck pain; arthritis, arthropathies, and other bodily pain; and neuropathies. We excluded patients who were not eligible for VHA benefits, did not have at least 1 inpatient or outpatient encounter in the year prior to their index opioid prescription, or were enrolled in hospice in the year before their index date.
Our study included the most commonly prescribed short-acting and long-acting opioids in the VHA system during the 10-year study period. Long-acting agents have either an inherently long serum half-life or a delivery vehicle that allows less-frequent administration (ie, at most, twice daily). The long-acting opioids in the present study included orally administered sustained-release morphine sulfate, methadone hydrochloride, controlled-release oxycodone hydrochloride, levorphanol tartrate, and fentanyl patches. Liquid methadone hydrochloride was excluded because this is the formulation used by the VHA in treating opiate addiction. Short-acting opioids (ie, codeine phosphate, hydrocodone, and oxycodone) used in our study were orally administered and were packaged as single agents or as combination products containing acetaminophen or aspirin. Opioids that were combined with other agents in elixirs (eg, promethazine/codeine syrup) and injectable opioids were excluded.
Patients began to contribute information to analyses from their index opioid prescription date and continued to do so until they experienced an overdose event, died, switched opioid agents, discontinued opioid treatment, entered hospice, became ineligible for VHA benefits, or reached the end of the study period, whichever came first. When a refill for an opioid prescription was dispensed on a date that occurred before the calculated end date of the prior prescription, the new prescription was assumed to begin the day after the calculated end date of the old prescription, and the days’ supply was accumulated. If a patient accumulated more than 180 days’ supply on a given day, the supply was truncated at 180 days. We provided a grace period of 1.5 times the days’ supply for any given prescription before censoring a patient for discontinuing their opioid medication (ie, patients must have refilled their opioid prescription within this grace period or they were classified as discontinuers at the end of the grace period).
To assess and control for the effect of the opioid dose, we converted each opioid agent to the morphine-equivalent dose following the method of Von Korff et al.17 We computed the morphine-equivalent mean daily dose by dividing the total quantity prescribed by days’ supply and converted the daily dose thus calculated into a corresponding morphine-equivalent dose. After the conversion, prescriptions in morphine-equivalent mean daily doses were categorized as 1 mg to less than 20 mg, 20 mg to less than 50 mg, 50 mg to less than 100 mg, and 100 mg or greater.
Our primary outcome of interest was unintentional overdoses that were coded as drug or medication poisonings of accidental intent using ICD-9 codes (E850.x-860.x) or undetermined intent (E980.x or drug poisoning [960.x-980.x] without an accompanying external cause of injury code). If an e-code indicated that the poisoning was self-inflicted (E950.x) or assault-related (E962.x), it was not counted as an event. Events were recorded at VA institutions and via billing to the VA from outside hospitals.
We obtained baseline demographic and clinical characteristics pertaining to the 12 months prior to and including the index date. Demographic variables included age, sex, race (white, black, other, and unknown or missing), and the percentage of any service-connected disability. Service-connected disability percentage is a measure of the severity of a disability, ranging from 0% to 100% and assigned in 10% increments. Zero percent service-connected disability indicates that an individual did not have such a disability. Clinical characteristics included prior falls and fractures, other medical diagnoses, and psychiatric diagnoses. We categorized VA health care utilization as the use of general mental health clinic services, services provided in the posttraumatic stress disorder clinic, and use of specific therapies, including intensive therapy, rehabilitation, and substance abuse disorder treatment; emergency department and urgent care visits; and inpatient hospitalizations. We characterized comedication with nonopioid agents, including selective cyclooxygenase 2 inhibitors and other nonsteroidal anti-inflammatory drugs.
Propensity score adjustment modeled the probability of receiving a long-acting opioid and included baseline covariates and significant 2-way interaction terms. We trimmed the top and bottom 1% in the propensity score distribution to eliminate patients who almost always would and those who would almost always not receive long-acting opioids. Cox proportional hazards regression models with multivariate and stabilized inverse probability treatment weights were used to assess risk adjusted for potential confounders.18 To correct the standard errors in the weighted analysis, we used the robust sandwich estimate of the covariance matrix (ie, the covs[aggregate] option and identification statement in the phreg procedure in SAS, version 9.2). Survival models were conducted in a piecewise fashion to account for variable hazards over time.
Sensitivity analyses excluded people with (1) methadone as the index opioid owing to concern that providers may preferentially prescribe methadone for pain control when unmeasured factors lead them to suspect a patient is at increased risk of abusing their medication and (2) overdose events for which an e-code was unavailable. All analyses were performed using SAS, version 9.2 (SAS Institute Inc).
Of the 840 606 eligible patients with chronic pain diagnoses within 90 days of their first opioid prescription, the majority of those who filled index prescriptions for opioid monotherapy received short-acting agents (n = 801 729); 18 887 filled prescriptions for long-acting opioids (Figure). New users of long-acting opioid monotherapy in our study represent 17.5% of all new long-acting opioid prescriptions received by VHA beneficiaries during our study period (ie, 83.2% of patients receiving long-acting opioids were given concomitant short-acting opioids). Thus, our cohort of incident users of monotherapy with long-acting opioids constituted 2.3% of patients who used the VA during the year before their index opioid prescription fill (and had a chronic pain indication within 90 days of their index dose).
Selected characteristics of the cohort are described in Table 1. Complete data are provided in eTable 1 in the Supplement. Most of the patients were men. The most common chronic pain diagnoses included osteoarthritis, back or neck pain, and other arthropathies. Hydrocodone was the most commonly prescribed opioid. Patients who were given prescriptions for long-acting opioids were more likely to receive higher daily doses than were those receiving short-acting opioids as well as to have back and neck pain, depression, anxiety, posttraumatic stress disorder, and substance use disorders. The long-acting opioid group also was more likely to be receiving concomitant antidepressants and benzodiazepines.
A total of 319 unintentional overdose events were observed during the study period among patients analyzed in an as-treated fashion (ie, according to our primary censoring criteria). Of these events, 282 occurred among patients initiating therapy with short-acting opioids and 37 among patients initiating therapy with long-acting opioids (Table 2). Approximately half of all the events occurred within the first 60 days after the start of opioid therapy. The crude rate of overdose events observed for both short-acting and long-acting opioids was higher during the first 2 weeks after opioid initiation than thereafter, but the heightened risk immediately after initiation therapy was far more marked among patients initiating therapy with long-acting opioids than for patents initiating therapy with short-acting opioids (Table 2).
The crude hazard ratio (HR) of unintentional overdose events during the study period was more than 2.5 times higher for persons initiating therapy with long-acting opioids (35 per 10 000 person-years) compared with persons initiating therapy with short-acting opioids (14 per 10 000 person-years) (HR, 2.84; 95% CI, 2.01-4.02). After adjustment for age, sex, and opioid dose, the patients receiving long-acting opioids still had a 2.5-fold higher risk of overdose (HR, 2.56; 95% CI, 1.67-3.93). Adjustment using all available covariates slightly reduced the relative risk, but the risk remained significantly higher among patients who initiated therapy with long-acting agents (HR, 2.33; 95% CI, 1.26-4.32). Overdose risk during the first 2 weeks after treatment initiation was more than 5-fold higher for patients who began receiving long-acting opioids compared with those receiving short-acting opioids (HR, 5.25; 1.88-14.72). Relative risk decreased to approximately 2-fold thereafter (higher with long-acting opioids) (Table 2). These general findings remained in sensitivity analyses that excluded patients who initiated opioid therapy with methadone and when analyses were restricted to patients with definitive e-coded event outcomes (eTable 2 and eTable 3 in the Supplement). Of the overdose events in our study, 23.2% were coded as opioid overdose events per se. Overdose risk was greater for patients initiating higher-dose therapy, with the risk among those receiving therapy with more than 50-mg equivalents of morphine being at more than twice the risk of overdose events compared with those receiving opioids at 1- to 20-mg equivalents.
To our knowledge, this is the first study to provide estimates of the relative risk of unintentional overdose events in relation to equianalgesic doses of long-acting vs short-acting opioids. In our study, which was restricted to patients with chronic medical conditions to reduce confounding by indication, patients who initiated opioid therapy with long-acting agents were at significantly higher risk of unintentional overdose events compared with those given prescriptions for shorter-acting agents. Risk was especially high shortly after opioid therapy began (5-fold higher) and remained elevated throughout the 1-year study period (relative risk was approximately 50% higher for patients receiving long-acting opioids beyond 60 days of continuous opioid use). The exceptionally high relative risk observed during the first 2 weeks of therapy was driven by a very high rate of unintentional overdose injury among patients receiving long-acting agents rather than a particularly low rate of injury among those receiving short-acting agents, making it less likely that differentially poor initial adherence to opioid therapy by patients receiving short-acting agents accounts for our findings.
The only other large observational study to examine opioid therapy in relation to overdose risk in a cohort of patients with chronic painful conditions, conducted by Dunn et al,10 did not distinguish between unintentional and intentional overdose events. In that study, which focused on patients who used long-acting opioids and assessed opioid overdose risk for patients who continued treatment for at least 90 days after initiating therapy, overdose risk was elevated with increasing mean daily doses. Consistent with the findings of Dunn et al10 and the observational study that separated prevalent from incident opioid users,11 we found an increasing risk of overdose among patients initiating higher doses of opioids. Similar to Paulozzi et al,19 we found that opioid overdoses represented approximately 1 of every 3 to 4 unintentional overdose events.
Our decision to restrict analyses to patients with documented diagnoses of chronic painful conditions who were beginning new opioid therapy confers several advantages over designs that include prevalent opioid users. The primary advantages include the ability to detect adverse events that occur soon after drug therapy is started, assess risks over time, and control for selection bias with baseline patient characteristics that are not influenced by the effects of opioid treatment. In addition, incident-user designs also mitigate potential selection bias owing to a drug-related history that might affect current treatment assignment. Because we were interested in adverse outcomes most directly linkable to the opioid regimen a patient was receiving, our primary analytic strategy was to censor people at the time of opioid regimen changes. Our decision to censor data for patients at treatment discontinuation and to use a proportional hazards analysis adjusts for differences in treatment persistence. By censoring patient follow-up as soon as an individual switched opioid agents or augmented therapy with a different opioid, we avoided the problematic comparison of patients who change treatment in response to adverse effects, refractory pain, or worsening symptoms (any of these factors might be indications of elevated overdose risk) with patients who do not change treatment.
Several caveats should be considered when interpreting the findings from the present study. First, as in all analyses relying on claims databases, we had limited ability to adjust for the severity of substance use disorders, physical illness, and psychiatric illness, all of which might place patients at heightened risk of unintentional overdose without regard to the opioid regimen. Patients who received long-acting opioids in our study were, in fact, more likely to have baseline risk factors for unintentional overdose injury (eg, more substance use disorders, psychiatric morbidity, and comedication with barbiturates, antidepressants, and sedative-hypnotics)20-30 compared with patients who received shorter-acting opioids. Nevertheless, after adjustment for these differences, unintentional overdose risk remained substantially higher, especially soon after the initiation of long-acting opioid therapy. Propensity scores offer an advantage in studies of rare outcomes (eg, unintentional overdose events) because propensity scores model the relationship of covariates and their interactions with the drug exposure (which is relatively frequent) and not directly with the study outcome (which is often rare), thereby mitigating the risk of overfitting in a traditional outcome model.31,32 As is the case for all observational studies, however, our ability to adjust fully for underlying risk at baseline depends on our ability to accurately classify baseline confounders and is compromised to the extent that measurement and reporting of conditions coded in claims are misclassified.33 However, if unmeasured baseline risk factors were responsible for the differential risk we observed, it is not clear why such risk would be so much higher soon after initiation of therapy than thereafter.
Second, we used administrative data and therefore did not directly measure opioid adherence. Using automated prescription data may, however, more accurately measure use of the medication than studies that rely on data from self-report surveys.34-36 A related point is that we defined drug exposure in our primary analysis in a way that seeks to capture how patients fill their prescriptions (ie, analyses are “as treated”), but in so doing admit possible selection bias owing to censoring.37 Nevertheless, our findings were robust to analyses in which exposure was defined using first treatment carried forward, which would tend to bias findings toward the null. In addition, it is possible that some patients we classified as incident users were, in fact, prevalent users of opioids if they were taking opioids prescribed outside the VA or illicitly at the time our data suggest incident use.
Third, the event rates that we report are necessarily underestimates of the actual unintentional overdose event rate because events may go unreported if patients with unintentional overdose events (1) die without coming to medical attention for their fatal event, (2) do not seek medical attention for a nonfatal event, or (3) receive medical attention that is neither treated nor reimbursed by the VHA system. Since we do not have any reason, a priori, to believe that ascertainment of outcomes would be biased with respect to duration of opioid action, the relative risks that we report are unlikely to be affected by this consideration.
Fourth, our findings apply to patients with chronic painful conditions who received care within the VHA system, the vast majority of whom are male and older than 50 years. The applicability of our findings to females, patients without chronic medical conditions, and younger patients is not known. Given the elevated prevalence of chronic pain and opioid use among a recent sample of active-duty infantry soldiers who are not seeking treatment, many of whom will become eligible for VA services upon separation from the military,38 future work focusing on opioid duration of action and overdose risk among recently separated military personnel may be warranted.
Despite the study’s limitations, we believe that our findings provide the first evidence that the risk of unintentional opioid overdose injury is related to the prescribed drug’s duration of action. If replicated in other cohorts, our findings suggest that clinicians weighing the benefits and risks of different opioid regimens should take into account not only the daily dose prescribed but also the duration of opioid action, favoring short-acting opioids whenever possible, especially during the first 2 weeks after initiation of therapy.
Accepted for Publication: November 4, 2014.
Corresponding Author: Matthew Miller, MD, ScD, Department of Health Sciences and Epidemiology, Northeastern University, 360 Huntington Ave, Room 316 RB, Boston, MA 02115 (ma.miller@neu.edu).
Published Online: February 16, 2015. doi:10.1001/jamainternmed.2014.8071.
Author Contributions: Dr Miller had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Miller, Barber, Hermos, Cho, Gagnon.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Miller, Hermos, Gagnon.
Critical revision of the manuscript for important intellectual content: Barber, Leatherman, Fonda, Hermos, Cho, Gagnon.
Statistical analysis: Leatherman, Fonda, Cho, Gagnon.
Obtained funding: Miller, Barber Cho, Gagnon.
Administrative, technical, or material support: Barber, Fonda, Hermos, Cho, Gagnon.
Study supervision: Miller, Hermos Cho, Gagnon.
Conflict of Interest Disclosures: None reported.
Funding/Support: Dr Miller and Ms Barber received support for this work from an investigator-initiated research grant (R21CE001605) from the Centers for Disease Control and Prevention (principal investigator, Dr Miller).
Role of the Funder/Sponsor: The Centers for Disease Control and Prevention had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
2.Webster
LR, Cochella
S, Dasgupta
N,
et al. An analysis of the root causes for opioid-related overdose deaths in the United States.
Pain Med. 2011;12(suppl 2):S26-S35.
PubMedGoogle ScholarCrossref 3.Modarai
F, Mack
K, Hicks
P,
et al. Relationship of opioid prescription sales and overdoses, North Carolina.
Drug Alcohol Depend. 2013;132(1-2):81-86.
PubMedGoogle ScholarCrossref 4.Franklin
GM, Mai
J, Wickizer
T, Turner
JA, Fulton-Kehoe
D, Grant
L. Opioid dosing trends and mortality in Washington State workers’ compensation, 1996-2002.
Am J Ind Med. 2005;48(2):91-99.
PubMedGoogle ScholarCrossref 5.Paulozzi
LJ, Budnitz
DS, Xi
Y. Increasing deaths from opioid analgesics in the United States.
Pharmacoepidemiol Drug Saf. 2006;15(9):618-627.
PubMedGoogle ScholarCrossref 6.Paulozzi
LJ, Ryan
GW. Opioid analgesics and rates of fatal drug poisoning in the United States.
Am J Prev Med. 2006;31(6):506-511.
PubMedGoogle ScholarCrossref 7.Hartung
DM, Middleton
L, Haxby
DG, Koder
M, Ketchum
KL, Chou
R. Rates of adverse events of long-acting opioids in a state Medicaid program.
Ann Pharmacother. 2007;41(6):921-928.
PubMedGoogle ScholarCrossref 8.Paulozzi
LJ, Kilbourne
EM, Shah
NG,
et al. A history of being prescribed controlled substances and risk of drug overdose death.
Pain Med. 2012;13(1):87-95.
PubMedGoogle ScholarCrossref 9.Gomes
T, Mamdani
MM, Dhalla
IA, Paterson
JM, Juurlink
DN. Opioid dose and drug-related mortality in patients with nonmalignant pain.
Arch Intern Med. 2011;171(7):686-691.
PubMedGoogle ScholarCrossref 10.Dunn
KM, Saunders
KW, Rutter
CM,
et al. Opioid prescriptions for chronic pain and overdose: a cohort study.
Ann Intern Med. 2010;152(2):85-92.
PubMedGoogle ScholarCrossref 11.Bohnert
AS, Valenstein
M, Bair
MJ,
et al. Association between opioid prescribing patterns and opioid overdose-related deaths.
JAMA. 2011;305(13):1315-1321.
PubMedGoogle ScholarCrossref 13.Xiang
Y, Zhao
W, Xiang
H, Smith
GA. ED visits for drug-related poisoning in the United States, 2007.
Am J Emerg Med. 2012;30(2):293-301.
PubMedGoogle ScholarCrossref 14.Gomes
T, Mamdani
MM, Dhalla
IA, Paterson
JM, Juurlink
DN. Opioids and dose-related deaths—association or causation?
Arch Intern Med. 2011;171(18):1688-1689.
Google ScholarCrossref 15.Bohnert
ASB, Ilgen
MA, Ignacio
RV, McCarthy
JF, Valenstein
M, Blow
FC. Risk of death from accidental overdose associated with psychiatric and substance use disorders.
Am J Psychiatry. 2012;169(1):64-70.
PubMedGoogle ScholarCrossref 16.Paulozzi
LJ, Kilbourne
EM, Desai
HA. Prescription drug monitoring programs and death rates from drug overdose.
Pain Med. 2011;12(5):747-754.
PubMedGoogle ScholarCrossref 17.Von Korff
M, Saunders
K, Thomas Ray
G,
et al. De facto long-term opioid therapy for noncancer pain.
Clin J Pain. 2008;24(6):521-527.
PubMedGoogle ScholarCrossref 18.Xu
S, Ross
C, Raebel
MA, Shetterly
S, Blanchette
C, Smith
D. Use of stabilized inverse propensity scores as weights to directly estimate relative risk and its confidence intervals.
Value Health. 2010;13(2):273-277.
PubMedGoogle ScholarCrossref 19.Paulozzi
LJ, Zhang
K, Jones
CM, Mack
KA. Risk of adverse health outcomes with increasing duration and regularity of opioid therapy.
J Am Board Fam Med. 2014;27(3):329-338.
PubMedGoogle ScholarCrossref 20.Toblin
RL, Paulozzi
LJ, Logan
JE, Hall
AJ, Kaplan
JA. Mental illness and psychotropic drug use among prescription drug overdose deaths: a medical examiner chart review.
J Clin Psychiatry. 2010;71(4):491-496.
PubMedGoogle ScholarCrossref 21.Darke
S, Hall
W, Kaye
S, Ross
J, Duflou
J. Hair morphine concentrations of fatal heroin overdose cases and living heroin users.
Addiction. 2002;97(8):977-984.
PubMedGoogle ScholarCrossref 22.McGregor
C, Darke
S, Ali
R, Christie
P. Experience of non-fatal overdose among heroin users in Adelaide, Australia: circumstances and risk perceptions.
Addiction. 1998;93(5):701-711.
PubMedGoogle ScholarCrossref 23.Darke
S, Ross
J, Hall
W. Overdose among heroin users in Sydney, Australia, II: responses to overdose.
Addiction. 1996;91(3):413-417.
PubMedGoogle ScholarCrossref 24.Darke
S, Ross
J, Hall
W. Overdose among heroin users in Sydney, Australia, I: prevalence and correlates of non-fatal overdose.
Addiction. 1996;91(3):405-411.
PubMedGoogle ScholarCrossref 25.Bohnert
AS, Prescott
MR, Vlahov
D, Tardiff
KJ, Galea
S. Ambient temperature and risk of death from accidental drug overdose in New York City, 1990-2006.
Addiction. 2010;105(6):1049-1054.
PubMedGoogle ScholarCrossref 26.Galea
S, Nandi
A, Coffin
PO,
et al. Heroin and cocaine dependence and the risk of accidental non-fatal drug overdose.
J Addict Dis. 2006;25(3):79-87.
PubMedGoogle ScholarCrossref 27.Coffin
PO, Galea
S, Ahern
J, Leon
AC, Vlahov
D, Tardiff
K. Opiates, cocaine and alcohol combinations in accidental drug overdose deaths in New York City, 1990-98.
Addiction. 2003;98(6):739-747.
PubMedGoogle ScholarCrossref 28.Kaye
S, Darke
S. Non-fatal cocaine overdose among injecting and non-injecting cocaine users in Sydney, Australia.
Addiction. 2004;99(10):1315-1322.
PubMedGoogle ScholarCrossref 29.Maloney
E, Degenhardt
L, Darke
S, Nelson
EC. Are non-fatal opioid overdoses misclassified suicide attempts? comparing the associated correlates.
Addict Behav. 2009;34(9):723-729.
PubMedGoogle ScholarCrossref 30.Tobin
KE, Latkin
CA. The relationship between depressive symptoms and nonfatal overdose among a sample of drug users in Baltimore, Maryland.
J Urban Health. 2003;80(2):220-229.
PubMedGoogle ScholarCrossref 31.Braitman
LE, Rosenbaum
PR. Rare outcomes, common treatments: analytic strategies using propensity scores.
Ann Intern Med. 2002;137(8):693-695.
PubMedGoogle ScholarCrossref 32.D’Agostino
RB
Jr. Propensity score methods for bias reduction in the comparison of a treatment to a non-randomized control group.
Stat Med. 1998;17(19):2265-2281.
PubMedGoogle ScholarCrossref 33.Grobbee
DE, Hoes
AW. Confounding and indication for treatment in evaluation of drug treatment for hypertension.
BMJ. 1997;315(7116):1151-1154.
PubMedGoogle ScholarCrossref 34.Piper
JM, Ray
WA, Griffin
MR. Effects of Medicaid eligibility expansion on prenatal care and pregnancy outcome in Tennessee.
JAMA. 1990;264(17):2219-2223.
PubMedGoogle ScholarCrossref 35.Strom
BL, Carson
JL. Use of automated databases for pharmacoepidemiology research.
Epidemiol Rev. 1990;12:87-107.
PubMedGoogle Scholar 36.West
SL, Savitz
DA, Koch
G, Strom
BL, Guess
HA, Hartzema
A. Recall accuracy for prescription medications: self-report compared with database information.
Am J Epidemiol. 1995;142(10):1103-1112.
PubMedGoogle Scholar 38.Toblin
RL, Quartana
PJ, Riviere
LA, Walper
KC, Hoge
CW. Chronic pain and opioid use in US soldiers after combat deployment.
JAMA Intern Med. 2014;174(8):1400-1401.
PubMedGoogle ScholarCrossref