SOC 4930 & SOC 505 - Week 08 - Tips for Paired Data

week-08-tips.R

# reprex of issues related to working with paired data
## dependencies
library(stlData)
library(dplyr)
library(ggplot2)
library(tidyr)

## data preparation
income <- stlIncome

income %>%
  select(geoID, mi10_inflate, mi15) %>%
  gather(key = year, value = medianInc, mi10_inflate, mi15) %>% 
  mutate(year = ifelse(year == "mi10_inflate", 2010, 2015)) %>% 
  arrange(geoID, year) -> incomeLong

incomeWide <- spread(incomeLong, key = year, value = medianInc)

## plot 1
ggplot() +
  geom_boxplot(data = incomeLong, mapping = aes(x = year, y = medianInc))

## plot 2
ggplot() +
  geom_boxplot(data = incomeLong, mapping = aes(x = year, y = medianInc, group = year))

## plot 3
ggplot() +
  geom_boxplot(data = incomeLong, mapping = aes(x = as.factor(year), y = medianInc, group = year))

## using variables with numeric variable names
mean(incomeWide$`2010`)