Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.

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)

district<-read_excel("district.xls")
data1<-district %>% select(DISTNAME,DPETSPEP,DPFPASPEP)
summary(data1)
##       DISTNAME       DPETSPEP       DPFPASPEP     
##  Length   :1207   Min.   : 0.00   Min.   : 0.000  
##  N.unique :1196   1st Qu.: 9.90   1st Qu.: 5.800  
##  N.blank  :   0   Median :12.10   Median : 8.900  
##  Min.nchar:   7   Mean   :12.27   Mean   : 9.711  
##  Max.nchar:  50   3rd Qu.:14.20   3rd Qu.:12.500  
##                   Max.   :51.70   Max.   :49.000  
##                                   NAs    :5

The missing values are in DPFPASPEP.

data2<-data1 %>% na.omit()

There are 1,202 observations remaining in data2.

cor(data2$DPFPASPEP,data2$DPETSPEP)
## [1] 0.3700234

There’s a weak correlation between special ed spending and students enrolled in special ed.