Exploratory Analysis II

Code for Quiz 8.

1. Load the R package we will use.

2. Quiz questions

Replace all the ???s. These are answers on your moodle quiz.

Run all the individual code chunks to make sure the answers in this file correspond with your quiz answers

After you check all your code chunks run then you can knit it. It won’t knit until the ??? are replaced

The quiz assumes that you have watched the videos, downloaded (to your examples folder) and worked through the exercises in exercises_slides-50-61.Rmd

  1. Pick one of your plots to save as your preview plot. Use the ggsave command at the end of the chunk of the plot that you want to preview.

Question: Modify slide 51

Create a plot with the mpg dataset

ggplot(data = mpg) + 
   geom_point(aes(x = displ, y = hwy)) +
   facet_wrap(facets = vars(manufacturer))

Question: Modify facet-ex-2

ggplot(mpg) + 
  geom_bar(aes(y = manufacturer)) + 
  facet_grid(vars(class), scales = "free_y", space = "free_y")

Question: spend_time

To help you complete this question use: the patchwork slides and

the vignette: https://patchwork.data-imaginist.com/articles/patchwork.html

Download the file spend_time.csv from moodle

spend_time contains 10 years of data on how many hours Americans spend each day on 5 activities

read it into spend_time

spend_time <- read_csv("https://estanny.com/static/week8/spend_time.csv")
p1  <- spend_time %>% filter(year == "2018")  %>% 
ggplot() + 
  geom_col(aes(x = activity, y = avg_hours, fill = activity)) +
  scale_y_continuous(breaks = seq(0, 6, by = 1)) +
  labs(subtitle = "Avg hours per day: 2018", x = NULL, y = NULL)
p1

Avg hours per day: 2028

Start with spend_time

p2  <- spend_time  %>% 
ggplot() + 
  geom_col(aes(x = year, y = avg_hours, fill = activity)) +
  labs(subtitle  = "Avg hours per day: 2010-2019", x = NULL, y = NULL)
p2

Use patchwork to display p1 on top of p2

p_all  <-  p1 / p2 

p_all

Start with p_all

p_all_no_legend  <- p_all & theme(legend.position = 'none')
p_all_no_legend

Start with p_all_no_legend

p_all_no_legend  +
 plot_annotation(title = "How much time Americans spent on selected activities", 
                  caption = "Source: American Time of Use Survey, https://data.bls.gov/cgi-bin/surveymost?tu")

Patchwork 2

use spend_time from last question patchwork slides Start with spend_time

p4  <- 
spend_time %>% filter(activity == "leisure/sports")  %>% 
ggplot() + 
  geom_point(aes(x = year, y = avg_hours)) +
  geom_smooth(aes(x = year, y = avg_hours)) +
  scale_x_continuous(breaks = seq(2010, 2019, by = 1)) +
  labs(subtitle = "Avg hours per day: leisure/sports", x = NULL, y = NULL)
p4

Start with p4

p5 <-  p4 + coord_cartesian(ylim = c(0, 6))
p5

Start with spend_time

p6   <- 
 spend_time  %>% 
ggplot() + 
  geom_point(aes(x = year, y = avg_hours, color = activity, group = activity)) +
  geom_smooth(aes(x = year, y = avg_hours, color = activity, group = activity)) +
  scale_x_continuous(breaks = seq(2010, 2019, by = 1)) +
  coord_cartesian(ylim = c(0, 6)) + 
  labs(x = NULL, y = NULL) 
p6

Use patchwork to display p4 and p5 on top of p6

( p4 | p5 ) / p6