library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.1     ✔ tibble    3.2.1
## ✔ lubridate 1.9.4     ✔ tidyr     1.3.1
## ✔ purrr     1.0.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)

district<-read_excel("district.xls")
bilingual<-district%>%select(DPETBILP,DPSTBIFP)
summary(district$DPETBILP)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    0.00    2.90    7.30   12.58   16.80  100.00
summary(district$DPSTBIFP)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's 
##   0.000   0.000   0.000   2.311   2.100  94.300       3
bilingual_clean<-bilingual%>%na.omit(.)
ggplot(bilingual_clean,aes(DPETBILP,DPSTBIFP))+geom_point()

cor(bilingual_clean$DPETBILP,bilingual_clean$DPSTBIFP)
## [1] 0.5593534