+
+
+
+
+
+
+
+
library(tidyverse)
+
## ── Attaching packages ─────────────────────────────────────────────────────────── tidyverse 1.3.0 ──
+
## ✓ ggplot2 3.3.0 ✓ purrr 0.3.3
+## ✓ tibble 2.1.3 ✓ dplyr 0.8.3
+## ✓ tidyr 1.0.2 ✓ stringr 1.4.0
+## ✓ readr 1.3.1 ✓ forcats 0.4.0
+
## ── Conflicts ────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
+## x dplyr::filter() masks stats::filter()
+## x dplyr::lag() masks stats::lag()
+
+
For the writeup, see CS506 Final Project - Deliverable 1.pdf
+
+
Cleaning Data
+
data <- read.csv(file = '../data/vaccinationrate_and_exemption.csv')
+spending_data_headers <- read.csv(file = '../data/healthExpendituresState.csv', skip = 1, nrows = 1, as.is = TRUE)
+spending_data = read.csv(file = '../data/healthExpendituresState.csv', skip = 3, nrows = 51, as.is = TRUE)
+colnames(spending_data) = spending_data_headers
+spending_data$`Total Health Spending` = gsub("\\$", "", spending_data$`Total Health Spending`)
+spending_data$`Total Health Spending` = as.numeric(gsub("\\,", "", spending_data$`Total Health Spending`))
+spending_data <- spending_data %>% rename(
+ States = Location
+)
+population_data <- read.csv(file = '../data/population_by_state.csv')
+population_data <- population_data[order(population_data$States),]
+population_data$Population = as.numeric(gsub("\\,", "", population_data$Population))
+all_data <- merge(data, spending_data, by="States")
+all_data <- merge(all_data, population_data, by="States")
+all_data <- all_data %>% rename(VaxRate = `X.`)
+all_data <- all_data %>% mutate(spending_per_capita = (`Total Health Spending`/Population)*1000000)
+
+
+
Exemptions
+
ggplot(all_data) +
+ #geom_point(aes(x = Exemption, y=X.)) +
+ geom_boxplot(aes(x = Exemption, y=VaxRate)) +
+ ylab("Vaccination Rate") +
+ ggtitle("Vaccination Rates in States with Various Exemptions")
+
+
+
+
Spending per Capita by State
+
model <- lm(VaxRate~spending_per_capita, data = all_data)
+summary(model)
+
##
+## Call:
+## lm(formula = VaxRate ~ spending_per_capita, data = all_data)
+##
+## Residuals:
+## Min 1Q Median 3Q Max
+## -5.7266 -1.5730 0.4823 1.4642 6.4950
+##
+## Coefficients:
+## Estimate Std. Error t value Pr(>|t|)
+## (Intercept) 8.628e+01 2.099e+00 41.096 <2e-16 ***
+## spending_per_capita 6.555e-04 2.549e-04 2.572 0.0132 *
+## ---
+## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
+##
+## Residual standard error: 2.402 on 49 degrees of freedom
+## Multiple R-squared: 0.1189, Adjusted R-squared: 0.1009
+## F-statistic: 6.613 on 1 and 49 DF, p-value: 0.01321
+
+
+
+
+
+
+