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]]
How does drug use vary by age, and are there any significant relationships, trends or patterns based on certain types of drugs?
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.
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.
What type of study is this (observational/experiment)?
This study is observational.
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.
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.
Age is the independent variable in this study.
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")