Data Preparation

For this final project proposal, I am planning to leverage the Drug Use by Age dataset from FiveThirtyEight https://github.com/fivethirtyeight/data. The reason I chose this dataset is because many people are addicted to drugs ranging from multiple ages and I wanted to explore which ages and what types of drugs certain people are addicted to.

# Load necessary libraries
#install.packages("psych")  # Install the psych package
library(psych)             # Load the psych package
## Warning: package 'psych' was built under R version 4.3.3
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.0     ✔ readr     2.1.5
## ✔ ggplot2   3.4.4     ✔ stringr   1.5.1
## ✔ lubridate 1.9.3     ✔ tibble    3.2.1
## ✔ purrr     1.0.2     ✔ tidyr     1.3.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ ggplot2::%+%()   masks psych::%+%()
## ✖ ggplot2::alpha() masks psych::alpha()
## ✖ dplyr::filter()  masks stats::filter()
## ✖ dplyr::lag()     masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
# Load the data
drug_data <- read.csv("https://raw.githubusercontent.com/fivethirtyeight/data/master/drug-use-by-age/drug-use-by-age.csv")

# View the structure of the data
str(drug_data)
## 'data.frame':    17 obs. of  28 variables:
##  $ age                    : chr  "12" "13" "14" "15" ...
##  $ n                      : int  2798 2757 2792 2956 3058 3038 2469 2223 2271 2354 ...
##  $ alcohol_use            : num  3.9 8.5 18.1 29.2 40.1 49.3 58.7 64.6 69.7 83.2 ...
##  $ alcohol_frequency      : num  3 6 5 6 10 13 24 36 48 52 ...
##  $ marijuana_use          : num  1.1 3.4 8.7 14.5 22.5 28 33.7 33.4 34 33 ...
##  $ marijuana_frequency    : num  4 15 24 25 30 36 52 60 60 52 ...
##  $ cocaine_use            : num  0.1 0.1 0.1 0.5 1 2 3.2 4.1 4.9 4.8 ...
##  $ cocaine_frequency      : chr  "5.0" "1.0" "5.5" "4.0" ...
##  $ crack_use              : num  0 0 0 0.1 0 0.1 0.4 0.5 0.6 0.5 ...
##  $ crack_frequency        : chr  "-" "3.0" "-" "9.5" ...
##  $ heroin_use             : num  0.1 0 0.1 0.2 0.1 0.1 0.4 0.5 0.9 0.6 ...
##  $ heroin_frequency       : chr  "35.5" "-" "2.0" "1.0" ...
##  $ hallucinogen_use       : num  0.2 0.6 1.6 2.1 3.4 4.8 7 8.6 7.4 6.3 ...
##  $ hallucinogen_frequency : num  52 6 3 4 3 3 4 3 2 4 ...
##  $ inhalant_use           : num  1.6 2.5 2.6 2.5 3 2 1.8 1.4 1.5 1.4 ...
##  $ inhalant_frequency     : chr  "19.0" "12.0" "5.0" "5.5" ...
##  $ pain_releiver_use      : num  2 2.4 3.9 5.5 6.2 8.5 9.2 9.4 10 9 ...
##  $ pain_releiver_frequency: num  36 14 12 10 7 9 12 12 10 15 ...
##  $ oxycontin_use          : num  0.1 0.1 0.4 0.8 1.1 1.4 1.7 1.5 1.7 1.3 ...
##  $ oxycontin_frequency    : chr  "24.5" "41.0" "4.5" "3.0" ...
##  $ tranquilizer_use       : num  0.2 0.3 0.9 2 2.4 3.5 4.9 4.2 5.4 3.9 ...
##  $ tranquilizer_frequency : num  52 25.5 5 4.5 11 7 12 4.5 10 7 ...
##  $ stimulant_use          : num  0.2 0.3 0.8 1.5 1.8 2.8 3 3.3 4 4.1 ...
##  $ stimulant_frequency    : num  2 4 12 6 9.5 9 8 6 12 10 ...
##  $ meth_use               : num  0 0.1 0.1 0.3 0.3 0.6 0.5 0.4 0.9 0.6 ...
##  $ meth_frequency         : chr  "-" "5.0" "24.0" "10.5" ...
##  $ sedative_use           : num  0.2 0.1 0.2 0.4 0.2 0.5 0.4 0.3 0.5 0.3 ...
##  $ sedative_frequency     : num  13 19 16.5 30 3 6.5 10 6 4 9 ...
# Check the first few rows of the data
head(drug_data)
##   age    n alcohol_use alcohol_frequency marijuana_use marijuana_frequency
## 1  12 2798         3.9                 3           1.1                   4
## 2  13 2757         8.5                 6           3.4                  15
## 3  14 2792        18.1                 5           8.7                  24
## 4  15 2956        29.2                 6          14.5                  25
## 5  16 3058        40.1                10          22.5                  30
## 6  17 3038        49.3                13          28.0                  36
##   cocaine_use cocaine_frequency crack_use crack_frequency heroin_use
## 1         0.1               5.0       0.0               -        0.1
## 2         0.1               1.0       0.0             3.0        0.0
## 3         0.1               5.5       0.0               -        0.1
## 4         0.5               4.0       0.1             9.5        0.2
## 5         1.0               7.0       0.0             1.0        0.1
## 6         2.0               5.0       0.1            21.0        0.1
##   heroin_frequency hallucinogen_use hallucinogen_frequency inhalant_use
## 1             35.5              0.2                     52          1.6
## 2                -              0.6                      6          2.5
## 3              2.0              1.6                      3          2.6
## 4              1.0              2.1                      4          2.5
## 5             66.5              3.4                      3          3.0
## 6             64.0              4.8                      3          2.0
##   inhalant_frequency pain_releiver_use pain_releiver_frequency oxycontin_use
## 1               19.0               2.0                      36           0.1
## 2               12.0               2.4                      14           0.1
## 3                5.0               3.9                      12           0.4
## 4                5.5               5.5                      10           0.8
## 5                3.0               6.2                       7           1.1
## 6                4.0               8.5                       9           1.4
##   oxycontin_frequency tranquilizer_use tranquilizer_frequency stimulant_use
## 1                24.5              0.2                   52.0           0.2
## 2                41.0              0.3                   25.5           0.3
## 3                 4.5              0.9                    5.0           0.8
## 4                 3.0              2.0                    4.5           1.5
## 5                 4.0              2.4                   11.0           1.8
## 6                 6.0              3.5                    7.0           2.8
##   stimulant_frequency meth_use meth_frequency sedative_use sedative_frequency
## 1                 2.0      0.0              -          0.2               13.0
## 2                 4.0      0.1            5.0          0.1               19.0
## 3                12.0      0.1           24.0          0.2               16.5
## 4                 6.0      0.3           10.5          0.4               30.0
## 5                 9.5      0.3           36.0          0.2                3.0
## 6                 9.0      0.6           48.0          0.5                6.5
# Summary statistics for quantitative variables
describe(drug_data)
##                         vars  n    mean      sd median trimmed    mad    min
## age*                       1 17    9.00    5.05    9.0    9.00   5.93    1.0
## n                          2 17 3251.06 1297.89 2798.0 3043.60 487.78 2223.0
## alcohol_use                3 17   55.43   26.88   64.6   56.95  23.87    3.9
## alcohol_frequency          4 17   33.35   21.32   48.0   34.13   5.93    3.0
## marijuana_use              5 17   18.92   11.96   20.8   19.11  17.94    1.1
## marijuana_frequency        6 17   42.94   18.36   52.0   43.60  11.86    4.0
## cocaine_use                7 17    2.18    1.82    2.0    2.14   2.82    0.0
## cocaine_frequency*         8 17    6.00    2.52    6.0    6.07   1.48    1.0
## crack_use                  9 17    0.29    0.24    0.4    0.29   0.15    0.0
## crack_frequency*          10 17    6.71    4.16    7.0    6.67   5.93    1.0
## heroin_use                11 17    0.35    0.33    0.2    0.33   0.30    0.0
## heroin_frequency*         12 17    9.00    5.05    9.0    9.00   5.93    1.0
## hallucinogen_use          13 17    3.39    2.79    3.2    3.27   3.85    0.1
## hallucinogen_frequency    14 17    8.41   15.00    3.0    5.93   1.48    2.0
## inhalant_use              15 17    1.39    0.93    1.4    1.37   1.19    0.0
## inhalant_frequency*       16 17    6.76    2.93    7.0    6.87   2.97    1.0
## pain_releiver_use         17 17    6.27    3.17    6.2    6.40   4.15    0.6
## pain_releiver_frequency   18 17   14.71    6.94   12.0   13.80   2.97    7.0
## oxycontin_use             19 17    0.94    0.61    1.1    0.95   0.89    0.0
## oxycontin_frequency*      20 17    7.35    4.57    7.0    7.27   5.93    1.0
## tranquilizer_use          21 17    2.81    1.75    3.5    2.81   2.08    0.2
## tranquilizer_frequency    22 17   11.74   11.49   10.0    9.53   4.45    4.5
## stimulant_use             23 17    1.92    1.41    1.8    1.90   1.78    0.0
## stimulant_frequency       24 17   31.15   85.97   10.0   10.90   2.97    2.0
## meth_use                  25 17    0.38    0.26    0.4    0.37   0.30    0.0
## meth_frequency*           26 17    7.06    4.15    7.0    7.00   4.45    1.0
## sedative_use              27 17    0.28    0.14    0.3    0.29   0.15    0.0
## sedative_frequency        28 17   19.38   24.83   10.0   14.83   8.90    3.0
##                            max  range  skew kurtosis     se
## age*                      17.0   16.0  0.00    -1.41   1.22
## n                       7391.0 5168.0  1.95     3.33 314.78
## alcohol_use               84.2   80.3 -0.63    -1.07   6.52
## alcohol_frequency         52.0   49.0 -0.36    -1.83   5.17
## marijuana_use             34.0   32.9 -0.14    -1.59   2.90
## marijuana_frequency       72.0   68.0 -0.50    -0.84   4.45
## cocaine_use                4.9    4.9  0.20    -1.65   0.44
## cocaine_frequency*        10.0    9.0 -0.22    -0.70   0.61
## crack_use                  0.6    0.6 -0.26    -1.85   0.06
## crack_frequency*          13.0   12.0 -0.08    -1.59   1.01
## heroin_use                 1.1    1.1  0.74    -0.73   0.08
## heroin_frequency*         17.0   16.0  0.00    -1.41   1.22
## hallucinogen_use           8.6    8.5  0.36    -1.35   0.68
## hallucinogen_frequency    52.0   50.0  2.19     3.07   3.64
## inhalant_use               3.0    3.0  0.12    -1.35   0.22
## inhalant_frequency*       11.0   10.0 -0.50    -1.04   0.71
## pain_releiver_use         10.0    9.4 -0.32    -1.51   0.77
## pain_releiver_frequency   36.0   29.0  1.74     2.64   1.68
## oxycontin_use              1.7    1.7 -0.21    -1.59   0.15
## oxycontin_frequency*      15.0   14.0  0.19    -1.48   1.11
## tranquilizer_use           5.4    5.2 -0.22    -1.53   0.43
## tranquilizer_frequency    52.0   47.5  2.56     6.05   2.79
## stimulant_use              4.1    4.1  0.10    -1.56   0.34
## stimulant_frequency      364.0  362.0  3.40    10.22  20.85
## meth_use                   0.9    0.9  0.14    -1.16   0.06
## meth_frequency*           14.0   13.0  0.06    -1.36   1.01
## sedative_use               0.5    0.5 -0.11    -0.89   0.03
## sedative_frequency       104.0  101.0  2.39     5.15   6.02
# Summary statistics by age group
describeBy(drug_data, group = drug_data$age)
## 
##  Descriptive statistics by group 
## group: 12
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## n                          2 1 2798.0 NA 2798.0  2798.0   0 2798.0 2798.0     0
## alcohol_use                3 1    3.9 NA    3.9     3.9   0    3.9    3.9     0
## alcohol_frequency          4 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## marijuana_use              5 1    1.1 NA    1.1     1.1   0    1.1    1.1     0
## marijuana_frequency        6 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## cocaine_use                7 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## crack_frequency           10 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## hallucinogen_use          13 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## hallucinogen_frequency    14 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## inhalant_use              15 1    1.6 NA    1.6     1.6   0    1.6    1.6     0
## inhalant_frequency        16 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## pain_releiver_use         17 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## pain_releiver_frequency   18 1   36.0 NA   36.0    36.0   0   36.0   36.0     0
## oxycontin_use             19 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## oxycontin_frequency       20 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## tranquilizer_use          21 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## tranquilizer_frequency    22 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## stimulant_use             23 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## stimulant_frequency       24 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## meth_use                  25 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## meth_frequency            26 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 13
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## n                          2 1 2757.0 NA 2757.0  2757.0   0 2757.0 2757.0     0
## alcohol_use                3 1    8.5 NA    8.5     8.5   0    8.5    8.5     0
## alcohol_frequency          4 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## marijuana_use              5 1    3.4 NA    3.4     3.4   0    3.4    3.4     0
## marijuana_frequency        6 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## cocaine_use                7 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## cocaine_frequency          8 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## crack_use                  9 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## crack_frequency           10 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## heroin_use                11 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## heroin_frequency          12 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## hallucinogen_use          13 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## hallucinogen_frequency    14 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## inhalant_use              15 1    2.5 NA    2.5     2.5   0    2.5    2.5     0
## inhalant_frequency        16 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## pain_releiver_use         17 1    2.4 NA    2.4     2.4   0    2.4    2.4     0
## pain_releiver_frequency   18 1   14.0 NA   14.0    14.0   0   14.0   14.0     0
## oxycontin_use             19 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## oxycontin_frequency       20 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## tranquilizer_use          21 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## tranquilizer_frequency    22 1   25.5 NA   25.5    25.5   0   25.5   25.5     0
## stimulant_use             23 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## stimulant_frequency       24 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## meth_use                  25 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## meth_frequency            26 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## sedative_use              27 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## sedative_frequency        28 1   19.0 NA   19.0    19.0   0   19.0   19.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 14
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## n                          2 1 2792.0 NA 2792.0  2792.0   0 2792.0 2792.0     0
## alcohol_use                3 1   18.1 NA   18.1    18.1   0   18.1   18.1     0
## alcohol_frequency          4 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## marijuana_use              5 1    8.7 NA    8.7     8.7   0    8.7    8.7     0
## marijuana_frequency        6 1   24.0 NA   24.0    24.0   0   24.0   24.0     0
## cocaine_use                7 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## cocaine_frequency          8 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## crack_use                  9 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## crack_frequency           10 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## hallucinogen_use          13 1    1.6 NA    1.6     1.6   0    1.6    1.6     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    2.6 NA    2.6     2.6   0    2.6    2.6     0
## inhalant_frequency        16 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## pain_releiver_use         17 1    3.9 NA    3.9     3.9   0    3.9    3.9     0
## pain_releiver_frequency   18 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## oxycontin_use             19 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## oxycontin_frequency       20 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## tranquilizer_use          21 1    0.9 NA    0.9     0.9   0    0.9    0.9     0
## tranquilizer_frequency    22 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## stimulant_use             23 1    0.8 NA    0.8     0.8   0    0.8    0.8     0
## stimulant_frequency       24 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## meth_use                  25 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## meth_frequency            26 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1   16.5 NA   16.5    16.5   0   16.5   16.5     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 15
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## n                          2 1 2956.0 NA 2956.0  2956.0   0 2956.0 2956.0     0
## alcohol_use                3 1   29.2 NA   29.2    29.2   0   29.2   29.2     0
## alcohol_frequency          4 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## marijuana_use              5 1   14.5 NA   14.5    14.5   0   14.5   14.5     0
## marijuana_frequency        6 1   25.0 NA   25.0    25.0   0   25.0   25.0     0
## cocaine_use                7 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## cocaine_frequency          8 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## crack_use                  9 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## crack_frequency           10 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## heroin_use                11 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## heroin_frequency          12 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## hallucinogen_use          13 1    2.1 NA    2.1     2.1   0    2.1    2.1     0
## hallucinogen_frequency    14 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## inhalant_use              15 1    2.5 NA    2.5     2.5   0    2.5    2.5     0
## inhalant_frequency        16 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## pain_releiver_use         17 1    5.5 NA    5.5     5.5   0    5.5    5.5     0
## pain_releiver_frequency   18 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## oxycontin_use             19 1    0.8 NA    0.8     0.8   0    0.8    0.8     0
## oxycontin_frequency       20 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## tranquilizer_use          21 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## tranquilizer_frequency    22 1    4.5 NA    4.5     4.5   0    4.5    4.5     0
## stimulant_use             23 1    1.5 NA    1.5     1.5   0    1.5    1.5     0
## stimulant_frequency       24 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## meth_use                  25 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## meth_frequency            26 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## sedative_use              27 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## sedative_frequency        28 1   30.0 NA   30.0    30.0   0   30.0   30.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 16
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## n                          2 1 3058.0 NA 3058.0  3058.0   0 3058.0 3058.0     0
## alcohol_use                3 1   40.1 NA   40.1    40.1   0   40.1   40.1     0
## alcohol_frequency          4 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## marijuana_use              5 1   22.5 NA   22.5    22.5   0   22.5   22.5     0
## marijuana_frequency        6 1   30.0 NA   30.0    30.0   0   30.0   30.0     0
## cocaine_use                7 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## cocaine_frequency          8 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## crack_use                  9 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## crack_frequency           10 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1   16.0 NA   16.0    16.0   0   16.0   16.0     0
## hallucinogen_use          13 1    3.4 NA    3.4     3.4   0    3.4    3.4     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_frequency        16 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## pain_releiver_use         17 1    6.2 NA    6.2     6.2   0    6.2    6.2     0
## pain_releiver_frequency   18 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## oxycontin_use             19 1    1.1 NA    1.1     1.1   0    1.1    1.1     0
## oxycontin_frequency       20 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## tranquilizer_use          21 1    2.4 NA    2.4     2.4   0    2.4    2.4     0
## tranquilizer_frequency    22 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## stimulant_use             23 1    1.8 NA    1.8     1.8   0    1.8    1.8     0
## stimulant_frequency       24 1    9.5 NA    9.5     9.5   0    9.5    9.5     0
## meth_use                  25 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## meth_frequency            26 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 17
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## n                          2 1 3038.0 NA 3038.0  3038.0   0 3038.0 3038.0     0
## alcohol_use                3 1   49.3 NA   49.3    49.3   0   49.3   49.3     0
## alcohol_frequency          4 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## marijuana_use              5 1   28.0 NA   28.0    28.0   0   28.0   28.0     0
## marijuana_frequency        6 1   36.0 NA   36.0    36.0   0   36.0   36.0     0
## cocaine_use                7 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## crack_frequency           10 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1   14.0 NA   14.0    14.0   0   14.0   14.0     0
## hallucinogen_use          13 1    4.8 NA    4.8     4.8   0    4.8    4.8     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## inhalant_frequency        16 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_use         17 1    8.5 NA    8.5     8.5   0    8.5    8.5     0
## pain_releiver_frequency   18 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## oxycontin_use             19 1    1.4 NA    1.4     1.4   0    1.4    1.4     0
## oxycontin_frequency       20 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## tranquilizer_use          21 1    3.5 NA    3.5     3.5   0    3.5    3.5     0
## tranquilizer_frequency    22 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## stimulant_use             23 1    2.8 NA    2.8     2.8   0    2.8    2.8     0
## stimulant_frequency       24 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## meth_use                  25 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## meth_frequency            26 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## sedative_use              27 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## sedative_frequency        28 1    6.5 NA    6.5     6.5   0    6.5    6.5     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 18
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## n                          2 1 2469.0 NA 2469.0  2469.0   0 2469.0 2469.0     0
## alcohol_use                3 1   58.7 NA   58.7    58.7   0   58.7   58.7     0
## alcohol_frequency          4 1   24.0 NA   24.0    24.0   0   24.0   24.0     0
## marijuana_use              5 1   33.7 NA   33.7    33.7   0   33.7   33.7     0
## marijuana_frequency        6 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## cocaine_use                7 1    3.2 NA    3.2     3.2   0    3.2    3.2     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## crack_frequency           10 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## heroin_use                11 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## heroin_frequency          12 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## hallucinogen_use          13 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## hallucinogen_frequency    14 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## inhalant_use              15 1    1.8 NA    1.8     1.8   0    1.8    1.8     0
## inhalant_frequency        16 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_use         17 1    9.2 NA    9.2     9.2   0    9.2    9.2     0
## pain_releiver_frequency   18 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## oxycontin_use             19 1    1.7 NA    1.7     1.7   0    1.7    1.7     0
## oxycontin_frequency       20 1   14.0 NA   14.0    14.0   0   14.0   14.0     0
## tranquilizer_use          21 1    4.9 NA    4.9     4.9   0    4.9    4.9     0
## tranquilizer_frequency    22 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## stimulant_use             23 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## stimulant_frequency       24 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## meth_use                  25 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## meth_frequency            26 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## sedative_use              27 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## sedative_frequency        28 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 19
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## n                          2 1 2223.0 NA 2223.0  2223.0   0 2223.0 2223.0     0
## alcohol_use                3 1   64.6 NA   64.6    64.6   0   64.6   64.6     0
## alcohol_frequency          4 1   36.0 NA   36.0    36.0   0   36.0   36.0     0
## marijuana_use              5 1   33.4 NA   33.4    33.4   0   33.4   33.4     0
## marijuana_frequency        6 1   60.0 NA   60.0    60.0   0   60.0   60.0     0
## cocaine_use                7 1    4.1 NA    4.1     4.1   0    4.1    4.1     0
## cocaine_frequency          8 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## heroin_use                11 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## heroin_frequency          12 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## hallucinogen_use          13 1    8.6 NA    8.6     8.6   0    8.6    8.6     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    1.4 NA    1.4     1.4   0    1.4    1.4     0
## inhalant_frequency        16 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## pain_releiver_use         17 1    9.4 NA    9.4     9.4   0    9.4    9.4     0
## pain_releiver_frequency   18 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## oxycontin_use             19 1    1.5 NA    1.5     1.5   0    1.5    1.5     0
## oxycontin_frequency       20 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## tranquilizer_use          21 1    4.2 NA    4.2     4.2   0    4.2    4.2     0
## tranquilizer_frequency    22 1    4.5 NA    4.5     4.5   0    4.5    4.5     0
## stimulant_use             23 1    3.3 NA    3.3     3.3   0    3.3    3.3     0
## stimulant_frequency       24 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## meth_use                  25 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## meth_frequency            26 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## sedative_use              27 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## sedative_frequency        28 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 20
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## n                          2 1 2271.0 NA 2271.0  2271.0   0 2271.0 2271.0     0
## alcohol_use                3 1   69.7 NA   69.7    69.7   0   69.7   69.7     0
## alcohol_frequency          4 1   48.0 NA   48.0    48.0   0   48.0   48.0     0
## marijuana_use              5 1   34.0 NA   34.0    34.0   0   34.0   34.0     0
## marijuana_frequency        6 1   60.0 NA   60.0    60.0   0   60.0   60.0     0
## cocaine_use                7 1    4.9 NA    4.9     4.9   0    4.9    4.9     0
## cocaine_frequency          8 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## crack_use                  9 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## crack_frequency           10 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## heroin_use                11 1    0.9 NA    0.9     0.9   0    0.9    0.9     0
## heroin_frequency          12 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## hallucinogen_use          13 1    7.4 NA    7.4     7.4   0    7.4    7.4     0
## hallucinogen_frequency    14 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## inhalant_use              15 1    1.5 NA    1.5     1.5   0    1.5    1.5     0
## inhalant_frequency        16 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_use         17 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## pain_releiver_frequency   18 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## oxycontin_use             19 1    1.7 NA    1.7     1.7   0    1.7    1.7     0
## oxycontin_frequency       20 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## tranquilizer_use          21 1    5.4 NA    5.4     5.4   0    5.4    5.4     0
## tranquilizer_frequency    22 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## stimulant_use             23 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## stimulant_frequency       24 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## meth_use                  25 1    0.9 NA    0.9     0.9   0    0.9    0.9     0
## meth_frequency            26 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## sedative_use              27 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## sedative_frequency        28 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 21
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## n                          2 1 2354.0 NA 2354.0  2354.0   0 2354.0 2354.0     0
## alcohol_use                3 1   83.2 NA   83.2    83.2   0   83.2   83.2     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   33.0 NA   33.0    33.0   0   33.0   33.0     0
## marijuana_frequency        6 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## cocaine_use                7 1    4.8 NA    4.8     4.8   0    4.8    4.8     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## heroin_use                11 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## heroin_frequency          12 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## hallucinogen_use          13 1    6.3 NA    6.3     6.3   0    6.3    6.3     0
## hallucinogen_frequency    14 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## inhalant_use              15 1    1.4 NA    1.4     1.4   0    1.4    1.4     0
## inhalant_frequency        16 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## pain_releiver_use         17 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_frequency   18 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## oxycontin_use             19 1    1.3 NA    1.3     1.3   0    1.3    1.3     0
## oxycontin_frequency       20 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## tranquilizer_use          21 1    3.9 NA    3.9     3.9   0    3.9    3.9     0
## tranquilizer_frequency    22 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## stimulant_use             23 1    4.1 NA    4.1     4.1   0    4.1    4.1     0
## stimulant_frequency       24 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## meth_use                  25 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## meth_frequency            26 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## sedative_use              27 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## sedative_frequency        28 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 22-23
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## n                          2 1 4707.0 NA 4707.0  4707.0   0 4707.0 4707.0     0
## alcohol_use                3 1   84.2 NA   84.2    84.2   0   84.2   84.2     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   28.4 NA   28.4    28.4   0   28.4   28.4     0
## marijuana_frequency        6 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## cocaine_use                7 1    4.5 NA    4.5     4.5   0    4.5    4.5     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## heroin_use                11 1    1.1 NA    1.1     1.1   0    1.1    1.1     0
## heroin_frequency          12 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## hallucinogen_use          13 1    5.2 NA    5.2     5.2   0    5.2    5.2     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## inhalant_frequency        16 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_use         17 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## pain_releiver_frequency   18 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## oxycontin_use             19 1    1.7 NA    1.7     1.7   0    1.7    1.7     0
## oxycontin_frequency       20 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## tranquilizer_use          21 1    4.4 NA    4.4     4.4   0    4.4    4.4     0
## tranquilizer_frequency    22 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## stimulant_use             23 1    3.6 NA    3.6     3.6   0    3.6    3.6     0
## stimulant_frequency       24 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## meth_use                  25 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## meth_frequency            26 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 24-25
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## n                          2 1 4591.0 NA 4591.0  4591.0   0 4591.0 4591.0     0
## alcohol_use                3 1   83.1 NA   83.1    83.1   0   83.1   83.1     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   24.9 NA   24.9    24.9   0   24.9   24.9     0
## marijuana_frequency        6 1   60.0 NA   60.0    60.0   0   60.0   60.0     0
## cocaine_use                7 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## cocaine_frequency          8 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## heroin_use                11 1    0.7 NA    0.7     0.7   0    0.7    0.7     0
## heroin_frequency          12 1   17.0 NA   17.0    17.0   0   17.0   17.0     0
## hallucinogen_use          13 1    4.5 NA    4.5     4.5   0    4.5    4.5     0
## hallucinogen_frequency    14 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## inhalant_use              15 1    0.8 NA    0.8     0.8   0    0.8    0.8     0
## inhalant_frequency        16 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## pain_releiver_use         17 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_frequency   18 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## oxycontin_use             19 1    1.3 NA    1.3     1.3   0    1.3    1.3     0
## oxycontin_frequency       20 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## tranquilizer_use          21 1    4.3 NA    4.3     4.3   0    4.3    4.3     0
## tranquilizer_frequency    22 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## stimulant_use             23 1    2.6 NA    2.6     2.6   0    2.6    2.6     0
## stimulant_frequency       24 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## meth_use                  25 1    0.7 NA    0.7     0.7   0    0.7    0.7     0
## meth_frequency            26 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1   17.5 NA   17.5    17.5   0   17.5   17.5     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 26-29
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## n                          2 1 2628.0 NA 2628.0  2628.0   0 2628.0 2628.0     0
## alcohol_use                3 1   80.7 NA   80.7    80.7   0   80.7   80.7     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   20.8 NA   20.8    20.8   0   20.8   20.8     0
## marijuana_frequency        6 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## cocaine_use                7 1    3.2 NA    3.2     3.2   0    3.2    3.2     0
## cocaine_frequency          8 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## crack_use                  9 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## crack_frequency           10 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## heroin_use                11 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## heroin_frequency          12 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## hallucinogen_use          13 1    3.2 NA    3.2     3.2   0    3.2    3.2     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## inhalant_frequency        16 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## pain_releiver_use         17 1    8.3 NA    8.3     8.3   0    8.3    8.3     0
## pain_releiver_frequency   18 1   13.0 NA   13.0    13.0   0   13.0   13.0     0
## oxycontin_use             19 1    1.2 NA    1.2     1.2   0    1.2    1.2     0
## oxycontin_frequency       20 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## tranquilizer_use          21 1    4.2 NA    4.2     4.2   0    4.2    4.2     0
## tranquilizer_frequency    22 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## stimulant_use             23 1    2.3 NA    2.3     2.3   0    2.3    2.3     0
## stimulant_frequency       24 1    7.0 NA    7.0     7.0   0    7.0    7.0     0
## meth_use                  25 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## meth_frequency            26 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## sedative_use              27 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## sedative_frequency        28 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 30-34
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   14.0 NA   14.0    14.0   0   14.0   14.0     0
## n                          2 1 2864.0 NA 2864.0  2864.0   0 2864.0 2864.0     0
## alcohol_use                3 1   77.5 NA   77.5    77.5   0   77.5   77.5     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   16.4 NA   16.4    16.4   0   16.4   16.4     0
## marijuana_frequency        6 1   72.0 NA   72.0    72.0   0   72.0   72.0     0
## cocaine_use                7 1    2.1 NA    2.1     2.1   0    2.1    2.1     0
## cocaine_frequency          8 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## heroin_use                11 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## heroin_frequency          12 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## hallucinogen_use          13 1    1.8 NA    1.8     1.8   0    1.8    1.8     0
## hallucinogen_frequency    14 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## inhalant_use              15 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## inhalant_frequency        16 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## pain_releiver_use         17 1    5.9 NA    5.9     5.9   0    5.9    5.9     0
## pain_releiver_frequency   18 1   22.0 NA   22.0    22.0   0   22.0   22.0     0
## oxycontin_use             19 1    0.9 NA    0.9     0.9   0    0.9    0.9     0
## oxycontin_frequency       20 1   11.0 NA   11.0    11.0   0   11.0   11.0     0
## tranquilizer_use          21 1    3.6 NA    3.6     3.6   0    3.6    3.6     0
## tranquilizer_frequency    22 1    8.0 NA    8.0     8.0   0    8.0    8.0     0
## stimulant_use             23 1    1.4 NA    1.4     1.4   0    1.4    1.4     0
## stimulant_frequency       24 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## meth_use                  25 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## meth_frequency            26 1   14.0 NA   14.0    14.0   0   14.0   14.0     0
## sedative_use              27 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## sedative_frequency        28 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 35-49
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
## n                          2 1 7391.0 NA 7391.0  7391.0   0 7391.0 7391.0     0
## alcohol_use                3 1   75.0 NA   75.0    75.0   0   75.0   75.0     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1   10.4 NA   10.4    10.4   0   10.4   10.4     0
## marijuana_frequency        6 1   48.0 NA   48.0    48.0   0   48.0   48.0     0
## cocaine_use                7 1    1.5 NA    1.5     1.5   0    1.5    1.5     0
## cocaine_frequency          8 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## crack_use                  9 1    0.5 NA    0.5     0.5   0    0.5    0.5     0
## crack_frequency           10 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## hallucinogen_use          13 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## hallucinogen_frequency    14 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## inhalant_use              15 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## inhalant_frequency        16 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## pain_releiver_use         17 1    4.2 NA    4.2     4.2   0    4.2    4.2     0
## pain_releiver_frequency   18 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## oxycontin_use             19 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## oxycontin_frequency       20 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## tranquilizer_use          21 1    1.9 NA    1.9     1.9   0    1.9    1.9     0
## tranquilizer_frequency    22 1    6.0 NA    6.0     6.0   0    6.0    6.0     0
## stimulant_use             23 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## stimulant_frequency       24 1   24.0 NA   24.0    24.0   0   24.0   24.0     0
## meth_use                  25 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## meth_frequency            26 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## sedative_use              27 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## sedative_frequency        28 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 50-64
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   16.0 NA   16.0    16.0   0   16.0   16.0     0
## n                          2 1 3923.0 NA 3923.0  3923.0   0 3923.0 3923.0     0
## alcohol_use                3 1   67.2 NA   67.2    67.2   0   67.2   67.2     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1    7.3 NA    7.3     7.3   0    7.3    7.3     0
## marijuana_frequency        6 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## cocaine_use                7 1    0.9 NA    0.9     0.9   0    0.9    0.9     0
## cocaine_frequency          8 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## crack_use                  9 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## crack_frequency           10 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## heroin_use                11 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## heroin_frequency          12 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## hallucinogen_use          13 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## hallucinogen_frequency    14 1   44.0 NA   44.0    44.0   0   44.0   44.0     0
## inhalant_use              15 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## inhalant_frequency        16 1    4.0 NA    4.0     4.0   0    4.0    4.0     0
## pain_releiver_use         17 1    2.5 NA    2.5     2.5   0    2.5    2.5     0
## pain_releiver_frequency   18 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## oxycontin_use             19 1    0.4 NA    0.4     0.4   0    0.4    0.4     0
## oxycontin_frequency       20 1   12.0 NA   12.0    12.0   0   12.0   12.0     0
## tranquilizer_use          21 1    1.4 NA    1.4     1.4   0    1.4    1.4     0
## tranquilizer_frequency    22 1   10.0 NA   10.0    10.0   0   10.0   10.0     0
## stimulant_use             23 1    0.3 NA    0.3     0.3   0    0.3    0.3     0
## stimulant_frequency       24 1   24.0 NA   24.0    24.0   0   24.0   24.0     0
## meth_use                  25 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## meth_frequency            26 1    9.0 NA    9.0     9.0   0    9.0    9.0     0
## sedative_use              27 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## sedative_frequency        28 1  104.0 NA  104.0   104.0   0  104.0  104.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
## ------------------------------------------------------------ 
## group: 65+
##                         vars n   mean sd median trimmed mad    min    max range
## age                        1 1   17.0 NA   17.0    17.0   0   17.0   17.0     0
## n                          2 1 2448.0 NA 2448.0  2448.0   0 2448.0 2448.0     0
## alcohol_use                3 1   49.3 NA   49.3    49.3   0   49.3   49.3     0
## alcohol_frequency          4 1   52.0 NA   52.0    52.0   0   52.0   52.0     0
## marijuana_use              5 1    1.2 NA    1.2     1.2   0    1.2    1.2     0
## marijuana_frequency        6 1   36.0 NA   36.0    36.0   0   36.0   36.0     0
## cocaine_use                7 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## cocaine_frequency          8 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## crack_use                  9 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## crack_frequency           10 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## heroin_use                11 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## heroin_frequency          12 1    3.0 NA    3.0     3.0   0    3.0    3.0     0
## hallucinogen_use          13 1    0.1 NA    0.1     0.1   0    0.1    0.1     0
## hallucinogen_frequency    14 1    2.0 NA    2.0     2.0   0    2.0    2.0     0
## inhalant_use              15 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## inhalant_frequency        16 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## pain_releiver_use         17 1    0.6 NA    0.6     0.6   0    0.6    0.6     0
## pain_releiver_frequency   18 1   24.0 NA   24.0    24.0   0   24.0   24.0     0
## oxycontin_use             19 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## oxycontin_frequency       20 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## tranquilizer_use          21 1    0.2 NA    0.2     0.2   0    0.2    0.2     0
## tranquilizer_frequency    22 1    5.0 NA    5.0     5.0   0    5.0    5.0     0
## stimulant_use             23 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## stimulant_frequency       24 1  364.0 NA  364.0   364.0   0  364.0  364.0     0
## meth_use                  25 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## meth_frequency            26 1    1.0 NA    1.0     1.0   0    1.0    1.0     0
## sedative_use              27 1    0.0 NA    0.0     0.0   0    0.0    0.0     0
## sedative_frequency        28 1   15.0 NA   15.0    15.0   0   15.0   15.0     0
##                         skew kurtosis se
## age                       NA       NA NA
## n                         NA       NA NA
## alcohol_use               NA       NA NA
## alcohol_frequency         NA       NA NA
## marijuana_use             NA       NA NA
## marijuana_frequency       NA       NA NA
## cocaine_use               NA       NA NA
## cocaine_frequency         NA       NA NA
## crack_use                 NA       NA NA
## crack_frequency           NA       NA NA
## heroin_use                NA       NA NA
## heroin_frequency          NA       NA NA
## hallucinogen_use          NA       NA NA
## hallucinogen_frequency    NA       NA NA
## inhalant_use              NA       NA NA
## inhalant_frequency        NA       NA NA
## pain_releiver_use         NA       NA NA
## pain_releiver_frequency   NA       NA NA
## oxycontin_use             NA       NA NA
## oxycontin_frequency       NA       NA NA
## tranquilizer_use          NA       NA NA
## tranquilizer_frequency    NA       NA NA
## stimulant_use             NA       NA NA
## stimulant_frequency       NA       NA NA
## meth_use                  NA       NA NA
## meth_frequency            NA       NA NA
## sedative_use              NA       NA NA
## sedative_frequency        NA       NA NA
# Frequency table for age group
age_table <- table(drug_data$age)
print(age_table)
## 
##    12    13    14    15    16    17    18    19    20    21 22-23 24-25 26-29 
##     1     1     1     1     1     1     1     1     1     1     1     1     1 
## 30-34 35-49 50-64   65+ 
##     1     1     1     1
# Print column names
colnames(drug_data)
##  [1] "age"                     "n"                      
##  [3] "alcohol_use"             "alcohol_frequency"      
##  [5] "marijuana_use"           "marijuana_frequency"    
##  [7] "cocaine_use"             "cocaine_frequency"      
##  [9] "crack_use"               "crack_frequency"        
## [11] "heroin_use"              "heroin_frequency"       
## [13] "hallucinogen_use"        "hallucinogen_frequency" 
## [15] "inhalant_use"            "inhalant_frequency"     
## [17] "pain_releiver_use"       "pain_releiver_frequency"
## [19] "oxycontin_use"           "oxycontin_frequency"    
## [21] "tranquilizer_use"        "tranquilizer_frequency" 
## [23] "stimulant_use"           "stimulant_frequency"    
## [25] "meth_use"                "meth_frequency"         
## [27] "sedative_use"            "sedative_frequency"
# Define a function to create scatterplots with linear regression lines
create_scatterplot <- function(data, drug_name) {
  ggplot(data, aes(x = age, y = !!sym(drug_name))) +
    geom_point(color = "darkblue") +
    geom_smooth(method = "lm", se = FALSE, color = "red") +
    labs(title = paste("Relationship between Age and", drug_name),
         x = "Age",
         y = paste("Percentage of", drug_name, "Users"))
}

# Create scatterplots for each drug
scatterplots <- lapply(names(drug_data)[3:ncol(drug_data)], function(drug_name) {
  create_scatterplot(drug_data, drug_name)
})

# Output the scatterplots
scatterplots
## [[1]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[2]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[3]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[4]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[5]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[6]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[7]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[8]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[9]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[10]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[11]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[12]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[13]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[14]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[15]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[16]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[17]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[18]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[19]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[20]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[21]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[22]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[23]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[24]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[25]]
## `geom_smooth()` using formula = 'y ~ x'

## 
## [[26]]
## `geom_smooth()` using formula = 'y ~ x'

# Bar plot of average marijuana use by age group
bar_plot <- ggplot(drug_data, aes(x = as.factor(age), y = marijuana_use)) +
  stat_summary(fun = mean, geom = "bar", fill = "skyblue") +
  labs(title = "Average Marijuana Use by Age Group (Bar Plot)",
       x = "Age Group",
       y = "Average Percentage of Users")

# Output the plots
list(bar_plot)
## [[1]]

Research question

How does drug use vary by age, and are there any significant relationships, trends or patterns based on certain types of drugs?

Cases

What are the cases, and how many are there?

The cases are different age groups, ranging from 12 to 65+. There are 17 cases in total.

Data collection

Describe the method of data collection.

The data was collected through surveys conducted by FiveThirtyEight. Source: National Survey on Drug Use and Health from the Substance Abuse and Mental Health Data Archive.

Type of study

What type of study is this (observational/experiment)?

This study is observational.

Data Source

If you collected the data, state self-collected. If not, provide a citation/link.

The data was collected by FiveThirtyEight and is publicly available on their GitHub repository. Source: National Survey on Drug Use and Health from the Substance Abuse and Mental Health Data Archive.

Dependent Variable

What is the response variable? Is it quantitative or qualitative?

The response variable is the percentage of people who use each drug within different age groups. It is quantitative.

Independent Variable(s)

Age is the independent variable in this study.

Relevant summary statistics

Provide summary statistics for each the variables. Also include appropriate visualizations related to your research question (e.g. scatter plot, boxplots, etc). This step requires the use of R, hence a code chunk is provided below. Insert more code chunks as needed.

# Summary statistics
summary(drug_data)
##      age                  n         alcohol_use    alcohol_frequency
##  Length:17          Min.   :2223   Min.   : 3.90   Min.   : 3.00    
##  Class :character   1st Qu.:2469   1st Qu.:40.10   1st Qu.:10.00    
##  Mode  :character   Median :2798   Median :64.60   Median :48.00    
##                     Mean   :3251   Mean   :55.43   Mean   :33.35    
##                     3rd Qu.:3058   3rd Qu.:77.50   3rd Qu.:52.00    
##                     Max.   :7391   Max.   :84.20   Max.   :52.00    
##  marijuana_use   marijuana_frequency  cocaine_use    cocaine_frequency 
##  Min.   : 1.10   Min.   : 4.00       Min.   :0.000   Length:17         
##  1st Qu.: 8.70   1st Qu.:30.00       1st Qu.:0.500   Class :character  
##  Median :20.80   Median :52.00       Median :2.000   Mode  :character  
##  Mean   :18.92   Mean   :42.94       Mean   :2.176                     
##  3rd Qu.:28.40   3rd Qu.:52.00       3rd Qu.:4.000                     
##  Max.   :34.00   Max.   :72.00       Max.   :4.900                     
##    crack_use      crack_frequency      heroin_use     heroin_frequency  
##  Min.   :0.0000   Length:17          Min.   :0.0000   Length:17         
##  1st Qu.:0.0000   Class :character   1st Qu.:0.1000   Class :character  
##  Median :0.4000   Mode  :character   Median :0.2000   Mode  :character  
##  Mean   :0.2941                      Mean   :0.3529                     
##  3rd Qu.:0.5000                      3rd Qu.:0.6000                     
##  Max.   :0.6000                      Max.   :1.1000                     
##  hallucinogen_use hallucinogen_frequency  inhalant_use   inhalant_frequency
##  Min.   :0.100    Min.   : 2.000         Min.   :0.000   Length:17         
##  1st Qu.:0.600    1st Qu.: 3.000         1st Qu.:0.600   Class :character  
##  Median :3.200    Median : 3.000         Median :1.400   Mode  :character  
##  Mean   :3.394    Mean   : 8.412         Mean   :1.388                     
##  3rd Qu.:5.200    3rd Qu.: 4.000         3rd Qu.:2.000                     
##  Max.   :8.600    Max.   :52.000         Max.   :3.000                     
##  pain_releiver_use pain_releiver_frequency oxycontin_use    oxycontin_frequency
##  Min.   : 0.600    Min.   : 7.00           Min.   :0.0000   Length:17          
##  1st Qu.: 3.900    1st Qu.:12.00           1st Qu.:0.4000   Class :character   
##  Median : 6.200    Median :12.00           Median :1.1000   Mode  :character   
##  Mean   : 6.271    Mean   :14.71           Mean   :0.9353                      
##  3rd Qu.: 9.000    3rd Qu.:15.00           3rd Qu.:1.4000                      
##  Max.   :10.000    Max.   :36.00           Max.   :1.7000                      
##  tranquilizer_use tranquilizer_frequency stimulant_use   stimulant_frequency
##  Min.   :0.200    Min.   : 4.50          Min.   :0.000   Min.   :  2.00     
##  1st Qu.:1.400    1st Qu.: 6.00          1st Qu.:0.600   1st Qu.:  7.00     
##  Median :3.500    Median :10.00          Median :1.800   Median : 10.00     
##  Mean   :2.806    Mean   :11.74          Mean   :1.918   Mean   : 31.15     
##  3rd Qu.:4.200    3rd Qu.:11.00          3rd Qu.:3.000   3rd Qu.: 12.00     
##  Max.   :5.400    Max.   :52.00          Max.   :4.100   Max.   :364.00     
##     meth_use      meth_frequency      sedative_use    sedative_frequency
##  Min.   :0.0000   Length:17          Min.   :0.0000   Min.   :  3.00    
##  1st Qu.:0.2000   Class :character   1st Qu.:0.2000   1st Qu.:  6.50    
##  Median :0.4000   Mode  :character   Median :0.3000   Median : 10.00    
##  Mean   :0.3824                      Mean   :0.2824   Mean   : 19.38    
##  3rd Qu.:0.6000                      3rd Qu.:0.4000   3rd Qu.: 17.50    
##  Max.   :0.9000                      Max.   :0.5000   Max.   :104.00
# Visualizations
# Boxplot of marijuana use by age
ggplot(drug_data, aes(x = as.factor(age), y = marijuana_use)) +
  geom_point() +
  labs(title = "Marijuana Use by Age Group",
       x = "Age Group",
       y = "Percentage of Users")

# Scatterplot of cocaine use by age
ggplot(drug_data, aes(x = age, y = cocaine_use)) +
  geom_point() +
  labs(title = "Cocaine Use by Age",
       x = "Age",
       y = "Percentage of Users")