library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.2.1     ✔ readr     2.2.0
## ✔ forcats   1.0.1     ✔ stringr   1.6.0
## ✔ ggplot2   4.0.3     ✔ tibble    3.3.1
## ✔ lubridate 1.9.5     ✔ tidyr     1.3.2
## ✔ purrr     1.2.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ 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
library(readxl)

load("NSDUH_2023.Rdata")
summary(data$BNGDRKMON)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##  0.0000  0.0000  0.0000  0.2154  0.0000  1.0000
summary(data$IRAMDEYR)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.     NAs 
##  0.0000  0.0000  0.0000  0.1149  0.0000  1.0000   11572
hist(data$BNGDRKMON)

hist(data$IRAMDEYR)

plot(data$BNGDRKMON,data$IRAMDEYR)

cor(data$BNGDRKMON,data$IRAMDEYR,use = "complete.obs")
## [1] 0.03124424