Introduction

The article where this data set originates from concentrates on baby boomers and their overall drug use.

Original article: https://fivethirtyeight.com/features/how-baby-boomers-get-high/

I wish to investigate the marijuana and alcohol usage of young adults (early 20s) until mid-30s.

#Retrieve Original Data

DrugUse <- "https://raw.githubusercontent.com/fivethirtyeight/data/master/drug-use-by-age/drug-use-by-age.csv"

DrugData <- read.csv(DrugUse)

#Subset of data concentrating on early 20s through mid-30s and those who use alcohol and marijuana
SubDD <- DrugData[c(9:14), c(1 ,3:6)]

#Change column names
setnames(SubDD, old = c("age", "alcohol.use", "alcohol.frequency", "marijuana.use","marijuana.frequency"), new = c("Age Groups", "Percent(ETOH Users)", "Times Used(ETOH)", "Percent(MJ Users)", "Times Used(MJ)"))

#Provide summary of subset data
summary(SubDD)
##    Age Groups Percent(ETOH Users) Times Used(ETOH) Percent(MJ Users)
##  20     :1    Min.   :69.70       Min.   :48.00    Min.   :16.40    
##  21     :1    1st Qu.:78.30       1st Qu.:52.00    1st Qu.:21.82    
##  22-23  :1    Median :81.90       Median :52.00    Median :26.65    
##  24-25  :1    Mean   :79.73       Mean   :51.33    Mean   :26.25    
##  26-29  :1    3rd Qu.:83.17       3rd Qu.:52.00    3rd Qu.:31.85    
##  30-34  :1    Max.   :84.20       Max.   :52.00    Max.   :34.00    
##  (Other):0                                                          
##  Times Used(MJ)
##  Min.   :52    
##  1st Qu.:52    
##  Median :56    
##  Mean   :58    
##  3rd Qu.:60    
##  Max.   :72    
##

Including Graphs

You can also embed plots, for example: ###Bar Graph ###Histogram Graph

hist(SubDD$`Times Used(ETOH)`)

hist(SubDD$`Times Used(MJ)`)

##Conclusion There is a significantly larger percentage of alcohol users than marijuana users. However, the frequency of marijuana use is on par or greater than alcohol. Suggesting those who do use marijuana use it and like it more than those who use alcohol.