library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.6.1
## Warning: package 'ggplot2' was built under R version 4.6.1
## Warning: package 'tibble' was built under R version 4.6.1
## Warning: package 'tidyr' was built under R version 4.6.1
## Warning: package 'readr' was built under R version 4.6.1
## Warning: package 'purrr' was built under R version 4.6.1
## Warning: package 'dplyr' was built under R version 4.6.1
## Warning: package 'stringr' was built under R version 4.6.1
## Warning: package 'forcats' was built under R version 4.6.1
## Warning: package 'lubridate' was built under R version 4.6.1
## ── 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)
## Warning: package 'readxl' was built under R version 4.6.1
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
which variable as missing values? DPFPASPEP has missing values.
data2<-data1 %>% na.omit()
1,207 observations in data1. 1,202 observations are left in data2.
cor(data2$DPFPASPEP,data2$DPETSPEP)
## [1] 0.3700234
my interpretation: positive correlation, but not a strong correlation. a positive weak correlation between spending on special education and the percent of students in special education.