Advertisement
Clinical research study| Volume 130, ISSUE 6, P678-687.e7, June 2017

Trends of Cannabis Use Disorder in the Inpatient: 2002 to 2011

Published:February 01, 2017DOI:https://doi.org/10.1016/j.amjmed.2016.12.035

      Abstract

      Objectives

      The nationwide prevalence of cannabis use/abuse has more than doubled from 2002 to 2011. Whether the outpatient trend is reflected in the inpatient setting is unknown. We examined the prevalence and incidence of cannabis abuse/dependence as determined by discharge coding in a 10-year (2002-2011) National Inpatient Sample, as well as various trends among demographics, comorbidities, and hospitalization outcomes.

      Methods

      Cannabis abuse/dependence was identified on the basis of International Classification of Diseases, 9th Revision, Clinical Modification codes 304.3* and 305.2* in adults aged 18 years or more. We excluded cases coded “in remission.” National estimates of trends and matched-regression analyses were conducted.

      Results

      Overall, 2,833,567 (0.91%) admissions with documented cannabis abuse/dependence were identified, patients had a mean age of 35.12 ± 0.06 years, 62% were male, and there was an increasing trend in prevalence from 0.52% to 1.34% (P <.001). The mean Charlson Comorbidity Index was 0.47 ± 0.006, and inpatient mortality was 0.41%. All of the above demonstrated an increasing trend (P <.001). Mean length of stay was 6.23 ± 0.06 days. The top primary discharge diagnoses were schizoaffective/mood disorders, followed by psychotic disorders and alcoholism. Asthma prevalence in nontobacco smokers had a steeper increase in the cannabis subgroup than in the noncannabis subgroup (P = .002). Among acute pancreatitis admissions, cannabis abusers had a shorter length of stay (−11%) and lower hospitalization costs (−7%) than nonabusers.

      Conclusion

      Cannabis abuse/dependence is on the rise in the inpatient population, with an increasing trend toward older and sicker patients with increasing rates of moderate to severe disability. Psychiatric disorders and alcoholism are the main associated primary conditions. Cannabis abuse is associated with increased asthma incidence in nontobacco smokers and decreased hospital resource use in acute pancreatitis admissions.

      Keywords

      Clinical Significance
      • Cannabis abuse/dependence prevalence in the inpatient increased from 0.5% to 1.3%.
      • Abuse/dependence has an increasing trend in older, sicker patients.
      • The most common coded comorbidities are psychiatric disease and alcoholism.
      • Cannabis abuse is associated with an increase in asthma numbers in nontobacco smokers.
      • Cannabis abuse is associated with shorter length of stay (−11%) and hospital costs (−7%) in acute pancreatitis admissions.
      Cannabis has been a highly debated and controversial issue nationwide. It is the most prevalent illicit substance in the United States and the second most widely smoked substance after tobacco, ranging from simple one-off recreational use to abuse and overt dependence with associated psychiatric illnesses and socioeconomic implications for individuals.
      • Budney A.J.
      • Moore B.A.
      Development and consequences of cannabis dependence.
      In the United States alone, conservative estimates reach 11 million users.
      • Moore B.A.
      • Augustson E.M.
      • Moser R.P.
      • Budney A.J.
      Respiratory effects of marijuana and tobacco use in a U.S. sample.
      Indications for its medicinal use and legalization vary from state to state, potentially influencing the geographic distribution of cannabis and its abuse. By 2016, 25 states had already passed medical marijuana laws, 16 of which were enacted between 1996 and 2011.
      Although there is scarce evidence of the prevalence of cannabis abuse and related medical comorbidities in the adult inpatient population, its association with psychiatric diseases is well known.
      • Brady K.
      • Casto S.
      • Lydiard R.B.
      • Malcolm R.
      • Arana G.
      Substance abuse in an inpatient psychiatric sample.
      The overall prevalence of cannabis use across the nation has more than doubled from 2002 (4.1%) to 2012 (9.5%), with up to a 31.3% prevalence among those aged 18 to 34 years. In addition, cannabis use disorder has increased to 2.9% in 2012 from 1.5% in 2001.
      • Hasin D.S.
      • Saha T.D.
      • Kerridge B.T.
      • et al.
      Prevalence of marijuana use disorders in the United States Between 2001-2002 and 2012-2013.
      Whether the cannabis use disorder trend is reflected in the inpatient setting is not known.
      In this analysis, we examine the prevalence of cannabis abuse over 10 years (2002-2011) among inpatient admissions using the Nationwide Inpatient Sample (NIS). We studied the trends and association of cannabis with various demographic factors, clinically related comorbidities, and hospital-associated outcomes, such as disability on discharge, length of stay, inpatient mortality, and costs.

      Materials and Methods

      Design and Data Source

      We performed a retrospective observational cohort study, using the Healthcare Cost and Utilization Project NIS database, sponsored by the Agency for Healthcare Research and Quality (AHRQ).
      HCUP Databases
      Healthcare Cost and Utilization Project (HCUP).
      The Healthcare Cost and Utilization Project NIS is the largest publicly available database in the United States and includes a 20% sample of all inpatient admissions from 1000 hospitals in more than 40 states. This sample averages 7 million discharges every year, representing more than 95% of the US population. Details regarding the design and validity of the NIS have been described.
      • Patel N.J.
      • Deshmukh A.
      • Pant S.
      • et al.
      Contemporary trends of hospitalization for atrial fibrillation in the United States, 2000 through 2010: implications for healthcare planning.
      We used the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes to identify discharges with cannabis abuse/dependence–coded diagnoses between 2002 and 2011.

      Study Population

      Within the inclusion criteria were all adults aged 18 years or more who had a diagnosis of cannabis abuse or dependence (excluding remission) in any of the discharge diagnosis fields (DX1-DX25). The cannabis abuse/dependence diagnoses (“continuous,” “episodic,” “unspecified”) were identified using the following ICD-9-CM codes, excluding those “in remission”: 304.30, 304.31, 304.32, 305.20, 305.21, and 305.22.
      • Rumalla K.
      • Reddy A.Y.
      • Mittal M.K.
      Association of recreational marijuana use with aneurysmal subarachnoid hemorrhage.
      For trend analysis, observations with missing data in age, gender, or inpatient mortality were excluded. Similar methodology was used in previous studies of chronic obstructive pulmonary disease/bronchiectasis, asthma, and tobacco smoking.
      • Patel N.J.
      • Deshmukh A.
      • Pant S.
      • et al.
      Contemporary trends of hospitalization for atrial fibrillation in the United States, 2000 through 2010: implications for healthcare planning.
      • Rumalla K.
      • Reddy A.Y.
      • Mittal M.K.
      Association of recreational marijuana use with aneurysmal subarachnoid hemorrhage.
      • Tan W.C.
      • Lo C.
      • Jong A.
      • et al.
      Marijuana and chronic obstructive lung disease: a population-based study.
      • Nanchal R.
      • Kumar G.
      • Majumdar T.
      • et al.
      Utilization of mechanical ventilation for asthma exacerbations: analysis of a national database.
      • Roelands J.
      • Jamison M.G.
      • Lyerly A.D.
      • James A.H.
      Consequences of smoking during pregnancy on maternal health.
      Variables (Supplementary Tables 1 and 2, available online) and statistical analyses are shown in the Online Supplement.

      Results

      The total number of admissions with cannabis abuse/dependence over the period 2002 to 2011 was 2,833,567 (0.91% of total admissions), of which 0.6% was coded as the primary discharge diagnosis. The proportion of all admissions involving cannabis abuse/dependence had an increasing trend (P <.001) from 0.52% (2002) to 1.34% (2011), representing a 260% increase over 10 years, along with a concurrent increasing incidence from 737 per million of total US population in 2002 to 1773 per million in 2011 (Figure 1). Cannabis dependence disorder overall comprised 17.3% of the cannabis subpopulation.
      Figure thumbnail gr1
      Figure 1Admissions with cannabis abuse/dependence trends overview. Left Y-axis: number of admissions with 95% CIs from population estimates. Right Y-axis: admissions per million of population (grey, continuous line), hospitalization costs (grey, dotted line) with 95% CIs from population estimates.

      Admissions with Cannabis Abuse/Dependence

      The mean age of cannabis abusers (35.12 ± 0.06 years) has increased by 3.2 years (P <.001). This is due to both decreasing trends in younger age groups (ages <50 years) and increasing trends in older groups (Table 1). Men were the predominant gender (62.5%) with an increasing trend (P <.001). Differences among races showed weaker linear trends, with whites (55.2%) and Hispanics having a decreasing trend, and blacks and Asians showing an increasing trend (P <.001). Native Americans had a strong linear increasing trend (P <.001), with a 300% increase (0.5% in 2002 to 1.5% in 2011). In terms of median income levels, the lowest and highest quartiles showed an increasing percentage of patients discharged with a cannabis abuse/dependence diagnosis, with middle-income quartiles steadily decreasing (P <.001).
      Table 1Main Trends in Cannabis Abuse/Dependence Disorder in the Inpatient Population
      VariableYearOverall TotalP Value for TrendTrend Direction
      2002200320042005200620072008200920102011
      No. of admissions (n ± SE)
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      158,218 ± 34,681175,550 ± 36,761225,742 ± 48,804233,771 ± 48,637272,267 ± 58,431300,954 ± 66,297303,503 ± 69,479351,495 ± 76,711390,582 ± 92,018421,484 ± 94,4432,833,567 ± 74,678
      Admissions prevalence
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      0.52% ± 0.04%0.57% ± 0.03%0.73% ± 0.03%0.75% ± 0.04%0.86% ± 0.04%0.95% ± 0.05%0.95% ± 0.06%1.11% ± 0.05%1.25% ± 0.07%1.34% ± 0.07%0.91% ± 0.02%<.001Increasing
      Admissions incidence per million of population every year
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      737 ± 161809 ± 1691028 ± 2221053 ± 2191212 ± 2601324 ± 2921320 ± 3021512 ± 3301661 ± 3911773 ± 397-<.001Increasing
      Cannabis dependence23.0%18.0%19.6%18.7%18.3%16.7%16.1%16.5%15.8%15.4%17.3%<.001Decreasing
      Cannabis abuse as primary discharge diagnosis0.69%0.36%0.58%0.48%0.43%0.41%0.38%0.38%0.43%0.24%0.4%<.001Decreasing
      Cannabis dependence as Primary discharge diagnosis0.24%0.20%0.17%0.17%0.18%0.17%0.13%0.14%0.13%0.17%0.2%<.001Decreasing
      Age
       Age (y)
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      33.37 ± 0.1433.33 ± 0.1533.83 ± 0.1434.18 ± 0.1534.69 ± 0.1735.05 ± 0.1635.22 ± 0.1735.82 ± 0.1736 ± 0.1636.59 ± 0.1735.12 ± 0.06<.001Increasing
       Age 18-29 y38.8%40.7%40.0%39.8%38.7%39.7%37.7%38.4%36.5%38.7%38.7%<.001Decreasing
       Age 30-49 y53.1%50.7%50.3%49.6%47.4%45.2%45.0%43.3%42.8%46.7%46.7%<.001Decreasing
       Age 50-69 y7.8%8.4%9.3%10.3%13.6%14.7%16.8%17.8%20.0%14.2%14.2%<.001Increasing
       Age ≥70 y0.21%0.24%0.27%0.27%0.33%0.37%0.49%0.52%0.64%0.40%0.4%<.001Increasing
      Female gender36.9%39.4%37.9%38.1%38.2%38.3%38.0%37.2%36.2%36.6%37.5%<.001Decreasing
      Weak linear trend (linear regression P >.05).
      Race
       White55.0%57.7%55.2%58.7%55.4%53.6%58.1%55.2%53.1%54.0%55.2%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Black31.0%29.3%32.7%27.4%30.0%32.8%29.3%30.5%33.9%30.3%31.0%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Hispanic9.3%10.0%8.3%9.1%9.6%8.7%8.4%8.9%8.5%9.5%9.0%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Asian/Pacific Islander0.8%0.8%0.6%0.7%0.5%0.6%0.7%0.8%0.8%0.8%0.7%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Native American0.5%0.3%0.7%0.8%1.4%1.6%0.9%1.1%0.9%1.5%1.1%<.001Increasing
       Other3.3%1.9%2.6%3.3%3.0%2.6%2.7%3.5%2.7%3.9%3.0%<.001Increasing
      Weak linear trend (linear regression P >.05).
      Charlson Comorbidity Index0.313 ± 0.0090.344 ± 0.0120.364 ± 0.010.379 ± 0.0110.414 ± 0.0120.441 ± 0.0160.477 ± 0.0140.523 ± 0.0150.57 ± 0.0180.6 ± 0.0170.47 ± 0.005<.001Increasing
      In-hospital death0.33%0.39%0.36%0.34%0.32%0.39%0.43%0.45%0.43%0.54%0.41%<.001Increasing
      Moderate-severe disability on discharge12.7%11.9%12.3%11.8%11.7%12.2%11.9%11.6%12.9%13.9%12.4%<.001Increasing
      Weak linear trend (linear regression P >.05).
      Length of stay (excludes inpatient deaths)
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      5.79 ± 0.344.90 ± 0.135.71 ± 0.385.26 ± 0.225.11 ± 0.175.17 ± 0.205.00 ± 0.145.18 ± 0.195.24 ± 0.205.20 ± 0.155.23 ± 0.06<.001Decreasing
      Weak linear trend (linear regression P >.05).
      Cost ($USD)
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      7828 ± 8596333 ± 2207506 ± 7826918 ± 4157022 ± 2607373 ± 3727144 ± 1997480 ± 2488203 ± 2668294 ± 2757515 ± 114<.001Increasing
      Income level
       0-25th percentileN/A37.5%42.1%40.5%39.8%41.9%39.8%42.3%41.1%41.4%40.9%<.001Increasing
      Weak linear trend (linear regression P >.05).
       26th-50th percentileN/A27.9%27.6%25.9%25.6%24.5%26.7%26.6%25.6%24.6%25.9%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       51th-75th PercentileN/A21.7%18.4%20.7%20.6%20.1%19.9%19.0%20.3%20.2%20.0%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       76th-100th percentileN/A13.0%11.9%13.0%14.0%13.5%13.7%12.1%13.0%13.8%13.1%<.001Increasing
      Weak linear trend (linear regression P >.05).
      Insurance
       Medicare11.9%11.6%12.4%12.9%13.6%12.8%14.0%14.7%14.9%15.7%13.8%<.001Increasing
       Medicaid35.1%35.4%36.7%37.4%36.4%33.5%33.9%36.9%36.6%35.9%35.8%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Private26.5%26.6%22.8%21.8%20.9%21.7%24.4%19.8%19.6%21.2%22.0%<.001Decreasing
       Self-pay19.3%18.3%19.5%20.2%20.7%22.5%18.9%20.1%20.5%19.3%20.0%<.001Increasing
      Weak linear trend (linear regression P >.05).
       No charge1.6%1.8%2.0%1.7%1.8%2.3%2.2%2.1%2.2%1.6%2.0%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Other5.6%6.3%6.5%5.9%6.6%7.2%6.6%6.4%6.3%6.3%6.4%<.001Increasing
      Weak linear trend (linear regression P >.05).
      AHRQ comorbidities
       Alcoholism37.3%28.3%28.4%28.0%27.6%27.1%28.8%27.7%27.2%26.6%28.1%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Chronic lung disease9.7%10.5%11.6%12.4%13.9%13.4%14.3%15.2%15.3%16.0%13.8%<.001Increasing
       Psychosis8.8%9.2%8.8%10.0%10.6%10.2%10.9%11.2%11.7%12.3%10.7%<.001Increasing
       Depression9.0%9.1%9.1%10.1%9.8%9.8%9.8%10.1%11.1%11.5%10.2%<.001Increasing
       Obesity3.1%3.5%4.0%4.1%4.5%4.8%5.9%6.5%6.6%7.5%5.5%<.001Increasing
       Liver disease3.0%3.3%3.3%3.0%3.0%3.0%3.4%3.6%3.7%4.4%3.5%<.001Increasing
       Renal failure0.6%0.8%0.8%1.0%1.4%1.9%2.2%2.6%2.8%3.3%2.0%<.001Increasing
       Congestive heart failure0.6%0.8%0.8%0.9%1.0%1.3%1.3%1.6%1.7%2.1%1.3%<.001Increasing
       AIDS0.9%0.7%1.1%0.9%0.9%0.7%0.9%1.0%0.9%0.8%0.9%.4834No trend
       Pulmonary hypertension0.2%0.1%0.2%0.2%0.2%0.3%0.4%0.5%0.6%0.8%0.4%<.001Increasing
       Cancer0.8%0.2%0.2%0.2%0.2%0.3%0.3%0.4%0.5%0.6%0.4%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Metastatic cancer0.2%0.2%0.2%0.2%0.2%0.3%0.3%0.4%0.5%0.7%0.4%<.001Increasing
       Diabetes mellitus without complications3.9%4.1%4.6%5.2%5.4%5.9%6.0%6.5%7.0%7.3%5.9%<.001Increasing
       Diabetes mellitus with complications0.6%0.7%0.8%0.8%0.9%1.0%1.2%1.4%1.5%2.0%1.2%<.001Increasing
      Select lung comorbidities
      These comorbidities were extracted per the ICD-9-CM diagnoses codes as outlined in Supplementary Table 1 (available online).
       Any COPD/bronchiectasis3.6%3.8%4.1%4.4%5.2%5.1%5.5%5.9%6.0%6.6%5.3%<.001Increasing
       Asthma7.1%7.9%8.5%9.3%10.0%9.6%10.1%10.8%10.9%11.0%9.9%<.001Increasing
       Tobacco smoking31.4%35.1%34.3%42.8%45.0%44.6%46.4%48.8%50.2%50.5%44.7%<.001Increasing
       PTSD2.59%2.46%2.97%3.07%2.70%2.83%3.47%4.10%3.84%4.18%3.38%<.001Increasing
       Chronic pancreatitis0.51%0.51%0.60%0.57%0.67%0.68%0.81%0.88%1.08%1.05%0.79%<.001Increasing
       Teaching hospital47.7%51.8%46.6%46.3%53.7%53.2%54.9%52.6%55.2%53.8%52.2%<.001Increasing
       Urban hospital87.7%87.2%86.9%87.6%88.5%88.9%88.9%89.0%90.0%89.1%88.6%<.001Increasing
      Hospital region
       Northeast23.3%19.5%24.0%23.6%22.8%22.1%21.6%23.6%23.9%21.8%22.7%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Midwest28.2%23.6%25.3%26.8%27.7%27.6%29.1%28.1%24.5%26.5%26.8%.1518No trend
       South32.3%35.8%34.6%32.1%33.3%35.6%34.9%32.2%34.8%33.3%33.9%.0669No trend
       West16.1%21.2%16.1%17.5%16.2%14.7%14.4%16.1%16.8%18.4%16.6%.0032Decreasing
      Weak linear trend (linear regression P >.05).
      Hospital bed size
       Small7.1%10.6%7.5%8.2%10.7%9.9%10.4%9.7%8.7%8.8%9.3%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Medium28.3%26.8%25.1%27.6%25.9%24.3%22.4%25.7%22.7%25.5%25.1%<.001Decreasing
       Large64.5%62.5%67.4%64.3%63.4%65.7%67.1%64.6%68.6%65.7%65.6%<.001Increasing
      Weak linear trend (linear regression P >.05).
      AHRQ = Agency for Healthcare Research and Quality; AIDS = acquired immunodeficiency syndrome; COPD = chronic obstructive pulmonary disease; N/A = not available; PTSD = post-traumatic stress disorder; SE = standard error.
      Weak linear trend (linear regression P >.05).
      Cannabis abuse/dependence admissions and % from total admissions, incidence per million of population, age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      These comorbidities were extracted per the ICD-9-CM diagnoses codes as outlined in Supplementary Table 1 (available online).
      Discharges with moderate-to-severe disability made up 12.4% of total cannabis abuse/dependence admissions over the 10-year period and showed a weak linear increasing trend (P <.001); further trends within this subset are described in Supplementary Table 3 (available online).
      Mean length of stay (excluding inpatient deaths) was 5.23 ± 0.06 days and had a weak decreasing linear trend, whereas hospitalization costs have been increasing (P <.001), averaging $7514 ± $114. Medicare and Medicaid as payers (49.6% combined) showed an increasing trend (P <.001).
      Clinically relevant, associated, chronic conditions were assessed for trends using AHRQ comorbidity indicators. Alcoholism (28.1%), chronic lung disease (13.8%), psychoses (10.7%), and depression (10.2%) constituted the majority of these comorbidities, with alcoholism showing a decreasing trend (decreased by 10.7%, P <.001), whereas depression and psychoses showed an increasing trend (P <.001). Post-traumatic stress disorder increased from 2.59% to 4.18% over 10 years (P <.001), with an overall prevalence of 3.38%. Of note, chronic lung disease
      • Elixhauser A.
      • Steiner C.
      • Kruzikas D.
      had a 1.7 times increase in the period of 10 years, with a strong increasing trend (P <.001). Other chronic conditions such as congestive heart failure, renal failure, pulmonary hypertension, and metastatic cancer also showed increasing trends (P <.001). The Charlson Comorbidity Index followed in concordance with these findings, almost doubling during the study period (0.313 in 2002 vs 0.600 in 2011; P <.001).
      The 20 most common primary discharge diagnoses among the patients who had cannabis abuse/dependence as a secondary diagnosis are outlined in Supplementary Figure (available online). More than 5500 unique primary discharge diagnoses (ICD-9-CM codes) were recorded in this subpopulation of patients. As shown in Supplementary Figure (available online), approximately 1 of 5 patients (19.3%) was admitted for a primary diagnosis of mood/schizoaffective disorder. Psychosis and other schizophrenia accounted for 5.7% of the primary discharge diagnoses, followed by drug-induced mental disorders (3.5%), alcohol-induced mental disorders (2.2%), and opioid dependence (1.5%). Acute pancreatitis (ICD-9-CM: 577.0) showed an increasing trend from 2002, entering the Top 20 list in 2007 (1.13%), continuing its incremental trend (P <.001), and reaching position number 9 in 2011 (1.75%). Overall, acute pancreatitis was the primary discharge diagnosis for 1.18% of all admissions with documented cannabis abuse/dependence. Given this finding, we conducted a secondary analysis of the cholelithiasis prevalence trend (ICD-9-CM: 574.0-574.9)
      • Kilic A.
      • Sheer A.
      • Shah A.S.
      • Russell S.D.
      • Gourin C.G.
      • Lidor A.O.
      Outcomes of cholecystectomy in US heart transplant recipients.
      in the cannabis abuse/dependence subpopulation and found a doubling in prevalence (0.4% in 2002 vs 0.9% in 2011; P <.001), with the trend having a similar slope as that of acute pancreatitis (0.00052 vs 0.00050; slope Δ P = .816). Furthermore, we tested the hypothesis that among all patients with a primary discharge diagnosis of acute pancreatitis, cannabis abuse/dependence could be more common among the chronic pancreatitis users. The findings were consistent with our hypothesis that cannabis use was 58% more likely in patients with chronic pancreatitis than in patients without the disease (odds ratio [OR], 1.58; 95% confidence interval [CI], 1.48-1.69; P <.001), after adjusting for Charlson Comorbidity Index, gender, age, alcoholism, depression, psychiatric disease, and metastatic cancer. Of note, the prevalence of chronic pancreatitis in the cannabis abuse/dependence population more than doubled over the study period (0.51% in 2002 to 1.05% in 2011, P <.001).

      Cannabis and Prevalence Trends of Respiratory Disease

      After exact matching (outlined in the Online Supplement), linear regression comparing the slopes of trends in prevalence of asthma and chronic obstructive pulmonary disease separately between cannabis and noncannabis users was done. There were no significant differences in chronic obstructive pulmonary disease trend prevalence (Figure 2A). However, in asthma trend prevalence, the nonsmokers group showed a significant difference in slopes (1.7 times steeper increase, P = .002) (Figure 2B), whereas the significant difference was lost (P = .230) (Figure 2B) in the tobacco smokers subgroup, suggesting that cannabis has an effect on asthma prevalence only among nontobacco smokers.
      Figure thumbnail gr2
      Figure 2Bivariate comparison among cannabis and tobacco smoking in lung disease trends. (A) Age- and gender-matched chronic obstructive pulmonary disease/bronchiectasis prevalence trends. (B) Chronic obstructive pulmonary disease/bronchiectasis prevalence trends are plotted for each year, and a linear regression line is fitted for each subgroup. There was no significant difference between the slopes of cannabis versus noncannabis users in both the tobacco and nontobacco groups. (B) Age- and gender-matched asthma prevalence trends. (A) Asthma prevalence trends are plotted for each year, and a linear regression line is fitted for each subgroup. Among nontobacco smokers, cannabis users had a steeper increase of asthma prevalence than noncannabis users. Among tobacco smokers, there was no significant difference between cannabis and noncannabis users.

      Length of Stay and Hospitalization Costs

      Matched regression analysis between cannabis and noncannabis abuse/dependence users was conducted separately for asthma and acute pancreatitis as the primary discharge diagnoses (Table 2). In the asthma model, length of stay was slightly decreased by 6% (time ratio, 0.94; 95% CI, 0.894-0.996; P = .036), whereas costs were not affected (OR, 1.01; 95%, CI, 0.96-1.06; P = .815) by cannabis abuse/dependence after correcting for alcoholism, tobacco smoking, race, and hospital-level variables. Among the acute pancreatitis admissions, both length of stay (time ratio, 0.89; 95% CI, 0.84-0.94; P <.001) and costs (OR, 0.93; 95% CI, 0.88-0.97; P = .003) were significantly decreased in the cannabis abuse/dependence group, while also adjusting for history of chronic pancreatitis.
      Table 2Effect of Cannabis Abuse/Dependence on Hospital Outcomes
      Primary Discharge DiagnosisLength of StayHospitalization Costs
      Time Ratio95% CIP ValueOR95% CIP Value
      Asthma0.9440.894-0.996.0361.010.96-1.06.815
      Chronic obstructive pulmonary disease0.980.90-1.07.6151.050.97-1.13.271
      Acute pancreatitis0.890.84-0.94<.0010.930.88-0.97.003
      Controls were matched on exact age, exact Charlson comorbidity index, gender, specific hospital, and year. Average number of matched controls per cannabis abuse/dependence patient: 3.34. Additional adjustment by multivariable regression was done for tobacco smoking, white race, alcoholism, hospital bed size, hospital region, urban hospital location, and hospital teaching status. For the acute pancreatitis model, adjustment for chronic pancreatitis was made. Note: A time ratio <1 is interpreted as an earlier discharge (shorter length of stay) and vice versa. The OR for hospitalization costs indicates percent change in the presence of cannabis abuse/dependence compared with its absence.
      CI = confidence interval; OR = odds ratio.

      Discussion

      Our analysis demonstrated an approximately 2-fold increase in the prevalence of cannabis abuse/dependence-related admissions, with increasing trends in older age groups and their associated comorbidities and increasing hospitalization costs. Furthermore, a substantial proportion of these admissions are related primarily to psychiatric disease (including substance abuse), as previously described.
      • Brady K.
      • Casto S.
      • Lydiard R.B.
      • Malcolm R.
      • Arana G.
      Substance abuse in an inpatient psychiatric sample.
      In addition, there was a trend toward abuse/dependence in both the poorest and most wealthy quartiles of the study population. Of note, the data also reveal an increased prevalence of acute pancreatitis as a primary discharge diagnosis over this 10-year period. Of note, nontobacco smokers who abuse cannabis showed a steeper incidence of asthma compared with noncannabis users. Finally, cannabis abusers being admitted for acute pancreatitis are independently associated with shorter length of stay and lower hospitalization costs.
      A clearly increasing trend of admissions involving cannabis abuse/dependence is evident. This could be explained partly from the legalization,
      • Hoffmann D.E.
      • Weber E.
      Medical marijuana and the law.
      • Cerda M.
      • Wall M.
      • Keyes K.M.
      • Galea S.
      • Hasin D.
      Medical marijuana laws in 50 states: investigating the relationship between state legalization of medical marijuana and marijuana use, abuse and dependence.
      to various extents, of cannabis as a medical or recreational substance, which in turn increases the prevalence of users in the general population and by extension the chance to observe them in the inpatient setting. This finding correlates with the results from a recent outpatient interview study.
      • Hasin D.S.
      • Saha T.D.
      • Kerridge B.T.
      • et al.
      Prevalence of marijuana use disorders in the United States Between 2001-2002 and 2012-2013.
      We observed that cannabis users who are admitted to the hospital are becoming older, accompanied by a larger comorbidity burden as evident by the Charlson Comorbidity Index score and various other chronic AHRQ comorbidities, a finding that is expected with increasing age. Because patients who are admitted are increasingly more clinically complicated, we can see a correlating increasing trend of hospitalization costs and a weakly increasing trend of moderate-to-severe disability-on-discharge rates over time, with ages >50 years tripling in percentage within this discharge disability subgroup. Medical use of cannabis is potentially another explanation for the aging population of inpatient cannabis users. Legalization of medical cannabis over the last 15 years has led to extended research and subsequently its use in a wider spectrum of medicine, such as chemotherapy-associated nausea and vomiting, multiple sclerosis, neuropathic pain, acquired immunodeficiency syndrome–related anorexia and cachexia, and advanced cancer pain.
      • Hazekamp A.
      • Ware M.A.
      • Muller-Vahl K.R.
      • Abrams D.
      • Grotenhermen F.
      The medicinal use of cannabis and cannabinoids–an international cross-sectional survey on administration forms.
      Indeed, the latter 2 have shown increasing trends in our analysis. Because these conditions are more likely to be seen in patients aged >50 years, this could explain the doubling of the 50- to 69-year age group over the last 10 years. Since there is no way of capturing “medical cannabis use” in the NIS while differentiating it from abuse and dependence, and given the simultaneous increase of conditions that cannabis is indicated as an adjunctive therapy,
      • Hazekamp A.
      • Ware M.A.
      • Muller-Vahl K.R.
      • Abrams D.
      • Grotenhermen F.
      The medicinal use of cannabis and cannabinoids–an international cross-sectional survey on administration forms.
      there could be a link between abuse and medicinal indication.
      • Caplan G.
      Medical marijuana: a study of unintended consequences.
      Such a link could be explored further in patients with these conditions who have been prescribed medicinal cannabis by following them for development of abuse/dependence over a period of years.
      Moreover, a socioeconomic/racial picture seems to emerge, with white men maintaining predominance, and the poor and the rich with increasing abuse/dependence versus the middle-class. A significantly increasing trend (tripling in 10 years) in the proportion of the Native American minority could reflect an overall increase in their outpatient use of cannabis. The West region has the lowest percent of cannabis abusers, even though the state of California legalized cannabis in 1996, 6 years before the beginning of our study period. Future studies that examine possible confounders, such as prevailing age and comorbidity difference between geographic sections, may help explain this finding.
      Approximately 1 of 3 inpatient cannabis users have a history of alcoholism. It is the most common AHRQ comorbidity we analyzed, but, in fact, this is decreasing over the years. Its popularity with cannabis use has been a known social issue
      • Gledhill-Hoyt J.
      • Lee H.
      • Strote J.
      • Wechsler H.
      Increased use of marijuana and other illicit drugs at US colleges in the 1990s: results of three national surveys.
      • Mohler-Kuo M.
      • Lee J.E.
      • Wechsler H.
      Trends in marijuana and other illicit drug use among college students: results from 4 Harvard School of Public Health College Alcohol Study surveys: 1993-2001.
      ; however, there is no clear reason why there is a decreasing trend. This could be due to a possible “substitution” effect, with cannabis replacing alcohol, among other substances.
      • Lucas P.
      • Walsh Z.
      • Crosby K.
      • et al.
      Substituting cannabis for prescription drugs, alcohol and other substances among medical cannabis patients: the impact of contextual factors.
      • Lucas P.
      • Reiman A.
      • Earleywine M.
      • et al.
      Cannabis as a substitute for alcohol and other drugs: a dispensary-based survey of substitution effect in Canadian medical cannabis patients.
      The increasing incidence of acute pancreatitis admissions within the cannabis abuse/dependence inpatient subpopulation was an unforeseen finding, considering the fact that alcoholism (a common cause of acute pancreatitis) has a decreasing trend. In an effort to evaluate another possible cause for the increasing trend of acute pancreatitis, we secondarily examined the trend of gallstone disease as a possible cause of pancreatitis. We discovered a significant doubling over 10 years (0.4%-0.9%) and a similar trend slope. Furthermore, it seems that an increasing proportion of patients among the cannabis group have chronic pancreatitis, which could explain the use of cannabis as a pain relief adjunct. In fact, patients with acute-on-chronic pancreatitis are 58% more likely to abuse cannabis, compared with those with nonchronic pancreatitis. Whether this increase can be explained by cannabis use being associated with higher body mass index and subsequently a higher chance for gallstone disease,
      • Bonfrate L.
      • Wang D.Q.
      • Garruti G.
      • Portincasa P.
      Obesity and the risk and prognosis of gallstone disease and pancreatitis.
      or simply by the general increase of obesity in the general population,
      • Ng M.
      • Fleming T.
      • Robinson M.
      • et al.
      Global, regional, and national prevalence of overweight and obesity in children and adults during 1980-2013: a systematic analysis for the Global Burden of Disease Study 2013.
      • Flegal K.M.
      • Carroll M.D.
      • Kit B.K.
      • Ogden C.L.
      Prevalence of obesity and trends in the distribution of body mass index among US adults, 1999-2010.
      is unknown and would nevertheless merit further analysis in future studies.
      Among the admissions of patients with asthma, we observed a nonclinically significant, yet statistically significant decrease of 6% in length of stay in patients with cannabis/dependence abuse. We expect this finding to be secondary to increased power of the study, rather than a true observation, especially considering the nonsignificant effect on hospitalization costs (Table 2). Similar findings are shown for chronic obstructive pulmonary disease admissions. On the contrary, among the acute pancreatitis admissions, results were more striking: Length of stay was significantly shorter (−11%), and hospitalization costs were significantly less (−7%) in patients with cannabis abuse/dependence, even after adjusting for history of chronic pancreatitis. We hypothesize 2 possible explanations for this finding: First, the half-life of cannabis in the plasma in chronic users is approximately 4 days
      • Johansson E.
      • Agurell S.
      • Hollister L.E.
      • Halldin M.M.
      Prolonged apparent half-life of delta 1-tetrahydrocannabinol in plasma of chronic marijuana users.
      ; therefore, this finding could be explained by the chronic use of cannabis, resulting in a more seamless hospital stay with easier pain control. Furthermore, we hypothesize that chronic cannabis users would rather try and go home earlier, by subjectively reporting less pain, to use their cannabis at home for further pain relief. Nevertheless, this finding would merit of further research in the future.
      Our study exhibits evidence of the association between cannabis use and increasing asthma prevalence among nontobacco smokers from the general population. Among tobacco smokers, the cannabis abuse/dependence did not have any significant contribution to increasing asthma prevalence. This suggests that tobacco smoking has a strong enough effect that any cannabis effect would be masked, although it was shown in other studies that cannabis and tobacco have an additive synergy in predisposing asthmatic attacks.
      • Tan W.C.
      • Lo C.
      • Jong A.
      • et al.
      Marijuana and chronic obstructive lung disease: a population-based study.
      • Aldington S.
      • Williams M.
      • Nowitz M.
      • et al.
      Effects of cannabis on pulmonary structure, function and symptoms.
      Cannabis did not seem to affect the prevalence trend of chronic obstructive pulmonary disease, regardless of tobacco smoking status. This could be secondary to tobacco smoking being a more chronic influence on a chronic disease such as chronic obstructive pulmonary disease, because tobacco has been readily available to the masses for a longer time than cannabis.

      Study Limitations

      The use of administrative databases has been increasing in clinical research, but as with any such database, they are prone to coding inaccuracies and deficiencies, and mainly rely on the clinician's documentation and the coder's experience.
      • O'Malley K.J.
      • Cook K.F.
      • Price M.D.
      • Wildes K.R.
      • Hurdle J.F.
      • Ashton C.M.
      Measuring diagnoses: ICD code accuracy.
      Specifically, although the cannabis ICD-9-CM codes were used previously in the literature, it has been shown that, in general, illicit drug abuse/dependence is suffering from under-reporting (ie, low sensitivity), with the positive predictive value ranging from 48% to 65%.
      • Kim H.M.
      • Smith E.G.
      • Stano C.M.
      • et al.
      Validation of key behaviourally based mental health diagnoses in administrative data: suicide attempt, alcohol abuse, illicit drug abuse and tobacco use.
      This has been attributed to the fact that coding for these diagnoses was occurring only when there was documentation of the clinician advising the patient for cessation or if the patient reported problem behaviors (which also explains the high prevalence of mood disorders in our study).
      • Kim H.M.
      • Smith E.G.
      • Stano C.M.
      • et al.
      Validation of key behaviourally based mental health diagnoses in administrative data: suicide attempt, alcohol abuse, illicit drug abuse and tobacco use.
      Naturally, it is expected that cannabis abuse would be under-reported by patients because of its controversial social and legal nature, making it difficult to accurately extrapolate the findings of this study to the outpatient setting and more generally to the whole population.
      • O'Malley K.J.
      • Cook K.F.
      • Price M.D.
      • Wildes K.R.
      • Hurdle J.F.
      • Ashton C.M.
      Measuring diagnoses: ICD code accuracy.
      The NIS does not provide data on a patient's last use of cannabis before admission; therefore, to increase our specificity, we excluded cases that were coded as “in remission.”

      Conclusions

      Cannabis abuse/dependence is on the rise in the inpatient population, with older and sicker patients becoming the majority with increased rates of moderate to severe disability on discharge. The inpatient cannabis subpopulation follows specific socioeconomic profiles and is mainly associated with mood disorders, whereas there is a signal of a possible relationship with disorders related to obesity, such as acute pancreatitis and gallstone disease. Cannabis abuse may be contributing to increased asthma incidence in nontobacco smokers. Finally, cannabis abuse is associated with decreased hospital resource use in acute pancreatitis admissions. Further studies that analyze subpopulations with these diseases between cannabis and noncannabis users are warranted.

      Appendix: Online Supplement

      Variables

      The Nationwide Inpatient Sample (NIS) variables were used to identify demographic information, including age, gender, and race. Hospitalization-related variables included primary payer (insurance types), bed size, region, teaching versus nonteaching, urban versus rural, and median household income for patient ZIP code (categorized by quartile from 2003 onwards). Select Agency for Healthcare Research and Quality comorbidity variables
      • Elixhauser A.
      • Steiner C.
      • Kruzikas D.
      from the NIS were extracted and analyzed for trends. The variables of chronic obstructive pulmonary disease/bronchiectasis, asthma, tobacco smoking, and post-traumatic stress disorder were identified according to International Classification of Diseases, 9th Revision, Clinical Modification codes outlined in Supplementary Table 1 (available online). A modified version of the Charlson Comorbidity Index was calculated for every observation

      Gordon M. Calculates comorbidity indices based on ICD-9/10 (comorbidities.icd10). 0.6.1 ed. Available at: https://github.com/gforge/comorbidities.icd10 2016. Accessed June 1, 2016.

      • Quan H.
      • Sundararajan V.
      • Halfon P.
      • et al.
      Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data.
      • Quan H.
      • Li B.
      • Couris C.M.
      • et al.
      Updating and validating the Charlson Comorbidity Index and score for risk adjustment in hospital discharge abstracts using data from 6 countries.
      (Supplementary Table 2, available online).
      Length of stay was estimated after excluding patients who died in the hospital. Total charges billed for hospital services per hospital stay were converted to costs using the cost-to-charge tables provided by the Healthcare Cost and Utilization Project (HCUP), weighted for missing data and adjusted for annual inflation.
      Disability-on-discharge severity was computed using the “DISPUniform” variable, as described in Supplementary Table 3 (available online).

      Statistical Analysis

      Categoric and continuous variables were compared using chi-square and t test/Mann-Whitney tests, where appropriate, respectively. Tests for trend included the Cochran-Armitage test
      • Armitage P.
      Tests for linear trends in proportions and frequencies.

      Signorell A. DescTools: Tools for Descriptive Statistics. Vol R package version 0.99.16 2016. Available at: https://cran.r-project.org/package=DescTools. Accessed March 21, 2016.

      in R and linear regression to assess the linearity of proportions for categoric variables, as well as the Cuzick nonparametric test in Stata (StataCorp LP, College Station, Tex) for continuous variables.
      • Cuzick J.
      A Wilcoxon-type test for trend.
      Categoric variables with P >.05 in linear regression of proportions against year were considered “weak linear trends,” even when the Cochran-Armitage test was significant for trend.
      Incidence of hospitalizations with cannabis abuse/dependence was calculated per million US population, by dividing the weighted national estimates of hospitalizations in a given year by the US census population estimate ≥18 years of age for the corresponding year.
      Population estimates - vintage 2009: annual estimates of the resident population by sex and selected age groups for the United States: April 1, 2000 to July 1, 2009 (NC-EST2009-02).
      Population estimates - vintage 2011: annual estimates of the resident population by sex and selected age groups for the United States: April 1, 2010 to July 1, 2011 (NC-EST2011-02).
      These analyses were completed in Stata MP 14.1 (StataCorp LP) and R 3.3.0.
      • Team R.C.
      R: A Language and Environment for Statistical Computing.
      The outcomes of length of stay and hospitalization costs were calculated with lognormal accelerated failure time model and a linear regression model, respectively. For length of stay, the event was defined as variable DIED = 0, meaning that the patient was discharged, whereas DIED = 1 implying in-hospital death was the censored observation. In this way, the survival model would allow to retain power without upfront exclusion of deaths from the analysis. The accelerated failure time model also is not subject to the proportional hazards assumption, in contrast to the Cox proportional hazards model. The time ratio calculated in this accelerated failure time model is essentially a percentage change of length of stay (ie, time ratio = 1.5 means 50% increase in length of stay). For hospitalization costs, given their exponential distribution, they were first log-transformed and the coefficients of the model exponentiated to produce odds ratios (ie, odds ratio = 1.5 means 50% increase in cost). The analysis for both outcomes was stratified by exact matching on gender, exact age in years, exact Charlson Comorbidity Index, and year-hospital-id. The matched controls (drawn from a 15% random sample of noncannabis abuse/dependence patients) in each matching stratum were then weighted as the inverse of the total controls in a stratum, essentially creating arithmetically 1:1-matched clusters. Each matching stratum corresponds to a cluster of matched cases, which was used to correct the regression coefficients' standard errors for using matched observations. Race information was collapsed to white versus nonwhite race to enhance power. Smoking and race were then included as additional confounders in the regression models, together with other hospital-level variables. This analysis was done in Stata 14, using the matching module “cem”.

      Blackwell M, Iacus SM, King G, Porro G. cem: Coarsened Exact Matching in Stata. Available at: http://gking.harvard.edu/files/gking/files/cem-stata.pdf. Accessed December 12, 2016.

      Likewise, a 15% random sample of noncannabis abuse/dependence patients were matched on exact age, gender, year, and tobacco smoking status with cannabis abuse/dependence patients, and the annual prevalence trends of chronic obstructive pulmonary disease/bronchiectasis and asthma were compared between cannabis and noncannabis users. A linear regression line was fitted on the prevalence (% cases) of the selected variable against years. The slopes (β term) of the fitted lines were then compared for significant difference. The latter analysis was done in GraphPad Prism 6.0f (GraphPad Software, La Jolla, Calif).
      All analyses used national estimates that were calculated using the TRENDWT NIS variable, as suggested by the HCUP. Calculation of standard errors for population estimates was done by taking into account possible missing hospitals among the cannabis subpopulation and the complex survey design of the NIS, by using the “SVY” commands in Stata (except for matched analyses), as recommended by the HCUP.
      • Houchens R.
      • Elixhauser A.
      Final Report on Calculating Nationwide Inpatient Sample (NIS) Variances for Data Years 2011 and Earlier.
      P values were a priori designated significant as <.025 for trends and <.05 for bivariate and other comparisons.
      Figure thumbnail fx1
      Supplementary FigureTop 20 primary discharge diagnoses in cannabis abuse/dependence. Graphical illustration of the distribution of various, arbitrarily categorized, primary discharge diagnoses (DX1) over the years 2002 to 2011. Purple: mood disorder excluding schizophrenia/psychoses. Green: alcoholism-related disorders or any drug withdrawal–related disorders (excluding opioids). Yellow: opioid abuse/dependence. Red: schizophrenia-related disorders. Blue: unspecified psychoses. Turquoise: acute Pancreatitis. White: chest pain, pneumonias, or any other conditions not listed.
      Supplementary Table 1International Classification of Diseases, 9th Revision, Clinical Modification Codes Use
      VariableICD-9-CM Codes
      Chronic obstructive pulmonary disease, bronchiectasis490, 491.0, 491.1, 491.2, 491.20, 491.21, 491.22, 491.8, 491.9, 492.0, 492.8, 494, 494.0, 494.1, 496
      Asthma493.xx
      Tobacco smoking305.1, V15.82
      Post-traumatic stress disorder309.81
      Chronic pancreatitis577.1
      ICD-9-CM = International Classification of Diseases, 9th Revision, Clinical Modification.
      Supplementary Table 2Comorbidity Scoring System
      Reference.
      • Brady K.
      • Casto S.
      • Lydiard R.B.
      • Malcolm R.
      • Arana G.
      Substance abuse in an inpatient psychiatric sample.
      ComorbiditiesUpdated Weights (2010)Original Charlson Weights
      Myocardial infarction01
      Congestive heart failure21
      Peripheral vascular disease01
      Cerebrovascular disease01
      Dementia21
      Chronic pulmonary disease11
      Rheumatologic disease11
      Peptic ulcer disease01
      Mild liver disease21
      Diabetes without chronic complications01
      Diabetes with chronic complications12
      Hemiplegia or paraplegia22
      Renal disease12
      Any malignancy, including leukemia and lymphoma22
      Moderate or severe liver disease43
      Metastatic solid tumor66
      AIDS/HIV46
      Maximum Comorbidity Score2429
      AIDS = acquired immunodeficiency syndrome; HIV = human immunodeficiency virus.
      Supplementary Table 3Moderate-to-Severe Disability on Discharge Subgroup
      VariableYearOverall TotalP Value for TrendTrend Direction
      2002200320042005200620072008200920102011
      Moderate-severe disability prevalence12.7%11.9%12.3%11.8%11.7%12.2%11.9%11.6%12.9%13.9%12.4%<.001Increasing
      Weak linear trend (linear regression P >.05).
      Moderate-severe disability subgroup (n)18,72219,64926,11125,84629,91934,42433,95538,41147,54955,21432,9799
      Cannabis dependence26.7%18.9%18.3%19.0%18.2%17.7%17.4%14.0%14.3%15.0%17.1%<.001Decreasing
      Age (y)
      Age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      34.78 ± 0.2835.24 ± 0.2835.9 ± 0.3936.35 ± 0.2537.5 ± 0.2937.56 ± 0.3638.28 ± 0.339.76 ± 0.3140.15 ± 0.2940.75 ± 0.2938.27 ± 0.11<.001Increasing
      Age 18-29 y34.9%36.0%34.7%34.1%30.9%32.7%32.5%29.4%28.5%27.9%31.3%<.001Decreasing
      Age 30-49 y55.0%51.7%50.7%50.6%50.2%47.5%44.6%42.6%42.3%40.3%46.1%<.001Decreasing
      Age 50-69 y9.5%11.5%13.9%14.5%18.0%18.9%21.8%26.5%27.4%29.9%21.4%<.001Increasing
      Age ≥70 y0.6%0.8%0.8%0.8%0.9%0.9%1.0%1.5%1.8%1.9%1.2%<.001Increasing
      Female31.0%33.9%30.7%32.2%33.6%31.0%33.6%32.0%32.9%32.3%32.3%.004Increasing
      Weak linear trend (linear regression P >.05).
      Race
       White55.4%60.7%55.7%60.0%58.8%54.8%63.5%58.8%54.9%54.8%57.3%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Black31.6%29.3%33.4%27.0%28.3%31.9%25.8%29.1%33.6%30.9%30.4%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Hispanic6.8%7.1%6.7%8.2%9.0%7.8%7.0%7.3%7.4%9.2%7.8%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Asian/Pacific Islander0.8%1.2%0.4%0.4%0.5%0.7%0.7%0.7%0.6%0.7%0.7%.4181No trend
       Native American0.6%0.1%0.8%1.0%1.1%1.5%0.9%1.2%0.7%1.2%1.0%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Other4.7%1.6%3.0%3.4%2.2%3.4%2.0%2.8%2.8%3.3%2.9%.4729No trend
      Charlson Comorbidity Index
      Age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      0.444 ± 0.0270.455 ± 0.0230.49 ± 0.0260.526 ± 0.0230.602 ± 0.0250.626 ± 0.0310.684 ± 0.0280.834 ± 0.030.892 ± 0.0340.946 ± 0.0370.707 ± 0.011<.001Increasing
      Length of stay (excluding inpatient deaths)
      Age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      7.98 ± 0.657.24 ± 0.397.69 ± 0.837.05 ± 0.347.07 ± 0.247.83 ± 0.617.07 ± 0.267.41 ± 0.317.82 ± 0.47.58 ± 0.247.49 ± 0.13<.001Decreasing
      Weak linear trend (linear regression P >.05).
      Cost ($USD)
      Age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      12410 ± 158810238 ± 55511946 ± 169410625 ± 73910947 ± 48112769 ± 126112049 ± 50812940 ± 51714949 ± 68014779 ± 63112788 ± 273<.001Increasing
      Income Level
       0-25th percentileN/A37.0%41.7%40.6%38.9%39.6%37.6%40.9%40.6%40.7%39.9%<.001Increasing
      Weak linear trend (linear regression P >.05).
       26th-50th percentileN/A26.5%27.5%25.9%26.5%23.8%27.7%27.5%26.3%24.9%26.2%.0002Decreasing
      Weak linear trend (linear regression P >.05).
       51th-75th percentileN/A22.9%18.4%20.6%20.8%23.0%19.9%19.0%20.2%20.4%20.5%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       76th-100th percentileN/A13.7%12.4%12.9%13.8%13.6%14.8%12.5%12.9%14.0%13.4%.0141Increasing
      Weak linear trend (linear regression P >.05).
      Insurance
       Medicare15.6%16.8%17.0%18.8%20.4%18.6%20.9%24.0%23.4%24.4%20.9%<.001Increasing
       Medicaid35.0%35.1%37.3%38.5%34.7%32.3%32.0%36.7%37.8%37.1%35.8%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Private24.6%27.2%21.4%20.8%21.0%21.3%24.8%18.8%18.7%19.6%21.2%<.001Decreasing
       Self-pay14.8%13.8%15.3%14.9%16.2%17.4%13.9%14.0%12.8%12.5%14.3%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       No charge3.8%1.6%1.5%1.6%2.1%4.7%1.8%1.5%1.5%1.1%2.0%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Other6.2%5.5%7.4%5.4%5.6%5.6%6.6%5.0%5.9%5.3%5.8%<.001Decreasing
      Weak linear trend (linear regression P >.05).
      AHRQ Comorbidities
       Alcoholism42.9%31.6%33.0%32.1%31.5%30.4%32.6%32.4%29.4%30.9%32.0%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Psychiatric12.9%13.5%12.7%15.7%14.9%13.6%15.1%14.2%15.9%15.9%14.7%<.001Increasing
       Chronic lung disease3.0%4.2%4.3%3.9%5.2%5.4%6.6%7.9%7.8%9.0%6.3%<.001Increasing
       Depression18.0%18.8%17.2%20.6%20.4%20.1%22.6%22.6%22.4%22.4%21.0%<.001Increasing
       Obesity10.9%12.1%12.8%14.4%15.0%14.3%15.9%18.4%18.2%19.7%16.1%<.001Increasing
       Liver disease3.7%3.9%4.0%3.4%3.8%3.4%4.4%5.1%5.1%6.2%4.5%<.001Increasing
       Renal failure0.9%0.8%0.9%1.4%2.0%2.5%2.9%3.8%4.2%5.5%3.0%<.001Increasing
       Congestive heart failure0.8%1.0%1.3%1.5%1.7%1.9%2.1%2.9%3.2%3.9%2.4%<.001Increasing
       AIDS1.2%0.8%1.4%1.2%1.2%1.0%1.4%1.4%1.1%1.3%1.2%.068No trend
       Metastatic cancer0.3%0.4%0.4%0.4%0.5%0.7%0.7%1.0%1.3%1.4%0.8%<.001Increasing
       Pulmonary hypertension0.1%0.2%0.3%0.2%0.3%0.6%0.9%1.0%1.2%1.6%0.8%<.001Increasing
       Cancer0.9%0.3%0.4%0.4%0.4%0.4%0.6%0.9%1.1%1.1%0.7%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Diabetes mellitus without complications4.5%5.8%6.2%7.2%7.7%7.8%7.9%9.6%10.0%10.6%8.3%<.001Increasing
       Diabetes mellitus with complications0.7%1.0%1.2%1.4%1.4%1.6%1.7%2.5%3.0%3.4%2.1%<.001Increasing
      Select lung comorbidities
      These comorbidities were extracted according to the ICD-9-CM diagnoses codes as outlined in Supplementary Table 1 (available online).
       Any COPD/bronchiectasis4.9%5.1%5.5%6.3%6.8%6.7%7.2%9.6%9.0%9.9%7.7%<.001Increasing
       Asthma6.4%7.7%8.1%8.8%9.0%8.2%8.8%9.5%10.0%10.3%9.0%<.001Increasing
       Tobacco smoking29.8%34.5%33.1%40.3%44.6%44.1%47.7%50.6%51.5%51.1%45.0%<.001Increasing
       Teaching Hospital49.9%50.4%47.3%43.8%53.1%57.6%53.1%49.0%54.8%56.4%52.3%<.001Increasing
      Weak linear trend (linear regression P >.05).
       Urban hospital89.7%86.5%87.3%86.9%88.1%88.2%88.6%90.4%90.1%89.8%88.8%<.001Increasing
      Weak linear trend (linear regression P >.05).
      Hospital Bed Size
       Small7.6%13.3%8.1%9.9%12.4%10.6%11.4%10.3%9.4%8.8%10.1%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Medium24.4%25.9%22.6%26.6%24.9%21.0%21.7%24.7%23.1%26.2%24.1%.2691No trend
       Large68.0%60.8%69.3%63.5%62.7%68.4%66.8%65.0%67.5%64.9%65.8%.0166Increasing
      Weak linear trend (linear regression P >.05).
      Hospital Region
       Northeast25.4%21.6%28.0%26.7%23.1%28.8%23.6%23.5%25.2%23.8%24.9%<.001Decreasing
      Weak linear trend (linear regression P >.05).
       Midwest30.3%22.1%27.1%27.3%27.1%26.7%28.8%29.5%26.1%26.9%27.2%.0899No trend
       South29.3%33.4%29.2%26.9%32.0%30.5%30.4%29.3%32.2%30.1%30.4%.0019Increasing
      Weak linear trend (linear regression P >.05).
       West15.0%22.8%15.7%19.2%17.8%14.0%17.3%17.7%16.5%19.2%17.5%.1067No trend
      The DISPUNIFORM variable was used to calculate moderate-to-severe disability on discharge cases after excluding inpatient deaths and “against medical advice” discharges. Values 2, 5, and 6 were considered “moderate-severe discharge,” which include transfers to short-term hospital, skilled nursing facility, intermediate care facility, another type of facility, and home health care. Note: Hospital deaths are excluded in table.
      AHRQ = Agency for Healthcare Research and Quality; AIDS = acquired immunodeficiency syndrome; COPD = chronic obstructive pulmonary disease.
      Weak linear trend (linear regression P >.05).
      Age, length of stay, hospitalization costs, and Charlson Comorbidity Index are presented with their linearized standard errors population estimates.
      These comorbidities were extracted according to the ICD-9-CM diagnoses codes as outlined in Supplementary Table 1 (available online).

      References

        • Budney A.J.
        • Moore B.A.
        Development and consequences of cannabis dependence.
        J Clin Pharmacol. 2002; 42: 28S-33S
        • Moore B.A.
        • Augustson E.M.
        • Moser R.P.
        • Budney A.J.
        Respiratory effects of marijuana and tobacco use in a U.S. sample.
        J Gen Intern Med. 2005; 20: 33-37
      1. 25 Legal Medical Marijuana States and DC.
        (Available at:) (Accessed June 26, 2016)
        • Brady K.
        • Casto S.
        • Lydiard R.B.
        • Malcolm R.
        • Arana G.
        Substance abuse in an inpatient psychiatric sample.
        Am J Drug Alcohol Abuse. 1991; 17: 389-397
        • Hasin D.S.
        • Saha T.D.
        • Kerridge B.T.
        • et al.
        Prevalence of marijuana use disorders in the United States Between 2001-2002 and 2012-2013.
        JAMA Psychiatry. 2015; 72: 1235-1242
        • HCUP Databases
        Healthcare Cost and Utilization Project (HCUP).
        (Available at:) (Accessed June 26, 2016)
      2. HCUP NIS Database Documentation.
        (Available at:) (Accessed June 26, 2016)
        • Patel N.J.
        • Deshmukh A.
        • Pant S.
        • et al.
        Contemporary trends of hospitalization for atrial fibrillation in the United States, 2000 through 2010: implications for healthcare planning.
        Circulation. 2014; 129: 2371-2379
        • Rumalla K.
        • Reddy A.Y.
        • Mittal M.K.
        Association of recreational marijuana use with aneurysmal subarachnoid hemorrhage.
        J Stroke Cerebrovasc Dis. 2016; 25: 452-460
        • Tan W.C.
        • Lo C.
        • Jong A.
        • et al.
        Marijuana and chronic obstructive lung disease: a population-based study.
        CMAJ. 2009; 180: 814-820
        • Nanchal R.
        • Kumar G.
        • Majumdar T.
        • et al.
        Utilization of mechanical ventilation for asthma exacerbations: analysis of a national database.
        Respir Care. 2014; 59: 644-653
        • Roelands J.
        • Jamison M.G.
        • Lyerly A.D.
        • James A.H.
        Consequences of smoking during pregnancy on maternal health.
        J Womens Health (Larchmt). 2009; 18: 867-872
        • Elixhauser A.
        • Steiner C.
        • Kruzikas D.
        Comorbidity Software Documentation.
        (Available at:) (Accessed June 26, 2016)
        • Kilic A.
        • Sheer A.
        • Shah A.S.
        • Russell S.D.
        • Gourin C.G.
        • Lidor A.O.
        Outcomes of cholecystectomy in US heart transplant recipients.
        Ann Surg. 2013; 258: 312-317
        • Hoffmann D.E.
        • Weber E.
        Medical marijuana and the law.
        N Engl J Med. 2010; 362: 1453-1457
        • Cerda M.
        • Wall M.
        • Keyes K.M.
        • Galea S.
        • Hasin D.
        Medical marijuana laws in 50 states: investigating the relationship between state legalization of medical marijuana and marijuana use, abuse and dependence.
        Drug Alcohol Depend. 2012; 120: 22-27
        • Hazekamp A.
        • Ware M.A.
        • Muller-Vahl K.R.
        • Abrams D.
        • Grotenhermen F.
        The medicinal use of cannabis and cannabinoids–an international cross-sectional survey on administration forms.
        J Psychoactive Drugs. 2013; 45: 199-210
        • Caplan G.
        Medical marijuana: a study of unintended consequences.
        McGeorge L Rev. 2012; 43: 127
        • Gledhill-Hoyt J.
        • Lee H.
        • Strote J.
        • Wechsler H.
        Increased use of marijuana and other illicit drugs at US colleges in the 1990s: results of three national surveys.
        Addiction. 2000; 95: 1655-1667
        • Mohler-Kuo M.
        • Lee J.E.
        • Wechsler H.
        Trends in marijuana and other illicit drug use among college students: results from 4 Harvard School of Public Health College Alcohol Study surveys: 1993-2001.
        J Am Coll Health. 2003; 52: 17-24
        • Lucas P.
        • Walsh Z.
        • Crosby K.
        • et al.
        Substituting cannabis for prescription drugs, alcohol and other substances among medical cannabis patients: the impact of contextual factors.
        Drug Alcohol Rev. 2016; 35: 326-333
        • Lucas P.
        • Reiman A.
        • Earleywine M.
        • et al.
        Cannabis as a substitute for alcohol and other drugs: a dispensary-based survey of substitution effect in Canadian medical cannabis patients.
        Addict Res Theory. 2013; 21: 435-442
        • Bonfrate L.
        • Wang D.Q.
        • Garruti G.
        • Portincasa P.
        Obesity and the risk and prognosis of gallstone disease and pancreatitis.
        Best Pract Res Clin Gastroenterol. 2014; 28: 623-635
        • Ng M.
        • Fleming T.
        • Robinson M.
        • et al.
        Global, regional, and national prevalence of overweight and obesity in children and adults during 1980-2013: a systematic analysis for the Global Burden of Disease Study 2013.
        Lancet. 2014; 384: 766-781
        • Flegal K.M.
        • Carroll M.D.
        • Kit B.K.
        • Ogden C.L.
        Prevalence of obesity and trends in the distribution of body mass index among US adults, 1999-2010.
        JAMA. 2012; 307: 491-497
        • Johansson E.
        • Agurell S.
        • Hollister L.E.
        • Halldin M.M.
        Prolonged apparent half-life of delta 1-tetrahydrocannabinol in plasma of chronic marijuana users.
        J Pharm Pharmacol. 1988; 40: 374-375
        • Aldington S.
        • Williams M.
        • Nowitz M.
        • et al.
        Effects of cannabis on pulmonary structure, function and symptoms.
        Thorax. 2007; 62: 1058-1063
        • O'Malley K.J.
        • Cook K.F.
        • Price M.D.
        • Wildes K.R.
        • Hurdle J.F.
        • Ashton C.M.
        Measuring diagnoses: ICD code accuracy.
        Health Serv Res. 2005; 40: 1620-1639
        • Kim H.M.
        • Smith E.G.
        • Stano C.M.
        • et al.
        Validation of key behaviourally based mental health diagnoses in administrative data: suicide attempt, alcohol abuse, illicit drug abuse and tobacco use.
        BMC Health Serv Res. 2012; 12: 18