--- title: "Homework 5: 2016 Elections" author: "Key" date: "`r Sys.time()`" output: html_document: preserve_yaml: true toc: true toc_float: true published: false --- # Instructions > For Homework 5, you will fill out this RMarkdown template. Throughout the document, create code (in chunks) and write text answers (outside chunks) to answer the provided questions in an analysis of King County election data from 2016. IMPORTANT: Do NOT add any additional code chunks, and do NOT modify any chunk options! This week's homework will become "part 1" for next week's homework. ```{r setup, include=FALSE} knitr::opts_chunk$set(echo = TRUE) library(readr) library(dplyr) library(tidyr) library(ggplot2) ``` # Homework 5 ## Importing the Data > Download the data from . It is a plain text file of data, about 60 MB in size. Values are separated with commas. Read the file into R. Note the `cache=TRUE` chunk option, which allows R to store the file between "knits" of the RMarkdown document and thus save memory/time. ```{r import_data, message=FALSE,warning=FALSE,cache=TRUE} king_county_elections_2016 <- read_csv("https://raw.githubusercontent.com/breonh/breonh.github.io/main/csss_508/homework/homework_5/king_county_elections_2016.txt") ``` ## Exploratory Data Analysis (EDA) > Use functions `str` and/or `summary` to look at the data. Describe the data in their current state. How many rows are there? What variables are there? What kinds of values do they take? Are the column types sensible? ```{r} dim(king_county_elections_2016) summary(king_county_elections_2016) ``` The data includes 643,163 rows and 9 columns. The variables include Precint, Race, LEG, CC, CG, CounterGroup, Party, CounterType, and SumOfCount, some of which are numeric and some of which are characters. The current column types (numeric vs. char) seem sensible. > This real-world election data is provided to you in "tidy" format! That is, each row is an observation: The number of votes given to a candidate/ballot measure/answer type in a given political race across voters in a precinct. We will ignore the variables `LEG`, `CC`, and `CG` as they are not of practical use. For the remaining variables, we will explore them graphically and attempt to figure out what they mean. Remember that in real world data work, you often have to get by with intuition or poking around online to figure out the nature of the data. > In each code chunk below, present a summary of each variable on it's own (such as a histogram, frequency table, barplot, etc.). If there are *many* categories, it's fine to print just the first few. If you create a figure, use ggplot. After, write down a one sentence interpretation for each summary. ### Precinct ```{r} table(king_county_elections_2016$Precinct) %>% head(20) ``` This variable appears to contain election precincts from King County, in slightly varying numbers. ### Race ```{r} table(king_county_elections_2016$Race) %>% head(20) ``` This variable seems to contain election races from King County in 2016. ### CounterGroup ```{r} table(king_county_elections_2016$CounterGroup) ``` Every value in this variable is "Total", which makes the variable useless! ### Party ```{r} ggplot(king_county_elections_2016,aes(Party))+ geom_bar()+theme_bw()+ylab("Count")+ggtitle("Number of Observations by Party") ``` This variable seems to contain political parties, although, "NP" (which probably stands for "No Party" or "Non-Partisan") is very prevalent. Perhaps this stands for races which are non-partisan. ### CounterType ```{r} table(king_county_elections_2016$CounterType) %>% head(10) ``` This variable seems to contain candidate names. > Notice something odd about CounterType in particular? It tells you what a given row of votes was for... but it also has `Registered Voters` and `Times Counted`. What are these values? ### SumOfCount ```{r} ggplot(king_county_elections_2016,aes(SumOfCount))+geom_histogram(bins=30)+ xlab("Number of Votes by Race and Precinct")+ylab("Count")+theme_bw() ``` This variable seems to include the number of votes by precinct and race. Many precincts have no votes for particular races, it appears. ## The quantities of interest > In this assignment (and the next), we will focus on three major races in Washington in 2016: > * "US President & Vice President" > * "Governor" > * "Lieutenant Governor" > With these races, we are interested in: > 1. **Turnout rates** for each of these races in each precinct. We will measure turnout as total number of submitted votes (including for a candidate, blank, write-in, or "over vote") divided by the number of registered voters. > 2. Differences between precincts *in Seattle* and precincts *elsewhere in King County*. > 3. Precinct-level support for the Democratic candidates in King County in 2016 for each contest. We will measure support as the percentage of votes in a precinct for the Democratic candidate out of all votes for candidates or write-ins. > You will answer Questions #1 and #2 in this homework (Question #3 will be completed in homework 6). The sections below describe steps you may want to take to answer Questions 1 and 2. I suggest loading `dplyr` and `tidyr` (in the very first code chunk of this Rmd) to start! ## Filtering down the data > For what we want to do, there are a lot of rows that are not useful. Start by filtering the dataset to only includes rows in which the race is one of: "US President & Vice President", "Governor", or "Lieutenant Governor". Save this subsetted dataset as a new object, called `king_county_elections_2016_Exec` ```{r} king_county_elections_2016_Exec <- king_county_elections_2016 %>% filter(Race %in% c("US President & Vice President", "Governor", "Lieutenant Governor")) ``` ## Seattle precincts > In our subsetted data, we want to add a "boolean" variable (TRUE or FALSE) for a precinct belonging to Seattle. The following code will create a vector of booleans. Using this code, add it to your dataset king_county_elections_2016_Exec as a new variable called "Seattle" `ifelse(substr(king_county_elections_2016_Exec$Precinct, start = 1, stop = 4) == "SEA ","Seattle","Not Seattle")` ```{r} king_county_elections_2016_Exec$Seattle <- ifelse(substr(king_county_elections_2016_Exec$Precinct, start = 1, stop = 4) == "SEA ","Seattle","Not Seattle") ``` ## Registered voters and turnout rates > We want to calculate turnout rates for each race. We define $Turnout = \frac{Total Votes}{Registered Voters}$, where total votes are listed in rows where the variable `CounterType` equals 'Times Counted', and registered votes are listed in rows where `CounterType` equals 'Registered Voters'. We can calculate turnout rates in three steps: Total votes by race/precinct, Registered votes by race/precinct, and finally turnout by race/precinct. Let's do it! > First, create a dataset called `Votes` in which you filter `king_county_elections_2016_Exec` to contain rows only where CounterType == 'Times Counted'. You should now have on row per Precint/Race. Add a variable called "TotalVotes" to the `Votes` dataset, which contains the number of total votes by race/precinct (currently in the "SumOfCount" variable in `Votes`). ```{r} Votes <- king_county_elections_2016_Exec %>% filter(CounterType == 'Times Counted') Votes$TotalVotes <- Votes$SumOfCount ``` > Second, create a dataset called `Registered` which contains rows only where CounterType == 'Registered Voters'. The "SumOfCount" variable in this dataset includes the number of registered votes by precinct/race. Add a variable called "Registered" to the `Votes` dataset, which contains the number of registered voters by race/precinct (currently in the "SumOfCount" variable in `Registered`). ```{r} Registered <- king_county_elections_2016_Exec %>% filter(CounterType == 'Registered Voters') Votes$Registered <- Registered$SumOfCount ``` > Third, create a new variable in `Votes` called "Turnout", which includes turnout calculated by dividing total votes by registered voters. Then, subset `Votes` to contain only the following variables: Precinct, Seattle, Race, TotalVotes, Registered, and Turnout. Display the first 10 rows of `Votes`, and print the number of rows/columns (there should be 7551 rows and 6 columns!). ```{r} Votes$Turnout <- Votes$TotalVotes/Votes$Registered Votes <- Votes %>% select(Precinct,Seattle,Race,TotalVotes,Registered,Turnout) Votes %>% head(10) dim(Votes) ``` ## Quick Plot! > Create ggplot histograms of turnout rates, first by "Race" and second by "Seattle". Do you notice any changes in turnout based on race or whether or not a precinct was located in Seattle? ```{r} ggplot(Votes,aes(Turnout))+theme_bw()+ geom_histogram(bins=30)+facet_wrap(~Race)+ xlab("Turnout Rate")+ylab("Count")+ ggtitle("Turnout Race by Executive Race", "2016 King County Election Results")+ xlim(c(0,1)) ggplot(Votes,aes(Turnout))+theme_bw()+ geom_histogram(bins=30)+facet_wrap(~Seattle)+ xlab("Turnout Rate")+ylab("Count")+ ggtitle("Turnout Race by Seattle vs. Non-Seattle Precincts", "2016 King County Election Results")+ xlim(c(0,1)) ``` Turnout seems quite consistent across different races; however, turnout in Seattle seems slightly higher on average than turnout in other parts of King County. > That's it for Homework 5! # Next Time > In the next homework, we'll continue this analysis!