---
title: "Bivariate Viz (Notes)"
format:
  html:
    toc: true
    toc-depth: 2
    embed-resources: true
---

```{r include = FALSE}
# This chunk just sets up some styling (eg: default size of our images)
knitr::opts_chunk$set(
  collapse = TRUE, 
  warning = FALSE,
  message = FALSE,
  fig.height = 2.75, 
  fig.width = 4.25,
  fig.align = 'center')
```



> **TODAY'S GOALS**
> 
> - Review univariate viz ideas and code.
> - Explore how to visualize *relationships* among 2 variables.


\
\
\
\


> **FOR MORE ON THESE TOPICS**
> 
> Watch:
> 
> - [ggplot demo](https://www.youtube.com/watch?v=1SYzVMH62yw&list=PLyEH7o09I467e8zck95awweg_bGuLzqjz&index=8) (Lisa Lendway)
> - [common ggplot mistakes](https://www.youtube.com/watch?v=j7F-IdJAOJs&list=PLyEH7o09I467e8zck95awweg_bGuLzqjz&index=8)(Lisa Lendway)
> 
> Read:     
>
> - [Visualizing amounts](https://clauswilke.com/dataviz/visualizing-amounts.html) (Wilke)
> - [Visualizing many distributions at once](https://clauswilke.com/dataviz/boxplots-violins.html) (Wilke)





\
\
\
\



## Review


Let's review some *univariate* concepts and code using our class survey data.
If the answers aren't at the top of your mind, don't fret!
We've barely started speaking this new language, and learned a ton of vocab last week, so you naturally won't remember it all.

```{r}
# Import data
survey <- read.csv("https://ajohns24.github.io/data/112/about_us_2024.csv")

# How many students have now filled out the survey?


# What type of variables do we have?

```


\
\
\
\



**EXAMPLE 1: Hangout preferences**

Students were asked, in that moment, where they'd most like to spend time outside.
How did they answer?
Was there a lot of agreement or a lot of variability in answers?
Build and interpret a plot that helps address these questions while reviewing:

- "code as communication"
- connecting with the components of a plot:       
    - set up a **frame**
    - add a **layer** / geometric element
    - change the **theme**, e.g. axis labels, color, fill

```{r}
# Attach a package needed to use the ggplot function


# Make a ggplot

```

```{r echo=FALSE, message=FALSE}
library(tidyverse)
```




\
\
\
\



**EXAMPLE 2: Temperature preferences**

Students were asked about their ideal outdoor temperature, in degrees Celsius.
How did they answer?
What was the typical response?
What was the range in responses?
Were there any outliers?
Build and interpret 2 plots that help address these questions.

```{r}

```





\
\
\
\


> **Subtleties**
> 
> Bar charts & histograms can appear pretty similar, but they do different things.
> 
> - Bar charts count up the number of observations of *each* outcome of a variable.
>     They're good for categorical variables, or quantitative variables with only a handful of possible outcomes.
> - Histograms count up the number of observations that fall into different numerical *ranges* of variable.
>     They're good for quantitative variables, especially those with many different observed outcomes.
> 






\
\
\
\




## New stuff


Thus far, we've been studying one variable at a time, using *univariate* plots.
But once we get a sense of how individual variables behave on their own, our questions often turn to *relationships* among variables.
For example, in our `hikes` data:

- How much `time` does it take to complete a hike? ------> How is `time` related to a hike's `elevation`? What about its `length`?
- How does difficult `rating` vary from hike to hike? -------> How is difficulty `rating` related to a hike's `ascent`?


\
\
\
\



**Exploring relationships**

Exploring univariate patterns often sparks follow-up questions about **relationships** between 2+ variables.
Often, but not always, variables take on specific roles:

- **response variable**: the variable whose variability we would like to explain (`time` to complete a hike)
- **predictors**:  variables that might explain some of the variability in the response (a hike's `elevation` or `length`)

Visualizations can help explore:     

- relationship *trends* (direction and form)   
- relationship *strength* (degree of variability from the trend)    
- *outliers* in the relationship






\
\
\
\



**EXAMPLE 3**

For each pair of variables below, sketch on paper a visualization of their relationship.
Focus on general viz process, don't worry about the exact details.
The data here are totally made up.

a. 3pm temperature (response) vs 9am temperature (predictor)       

```{r}
data.frame(temp_3pm = c(24, 26, 20, 15, 15, 15), temp_9am = c(14, 18, 15, 13, 11, 11))
```
    
    

b. 3pm temperature (response) vs location (predictor)        

```{r}
weather <- data.frame(temp_3pm = c(24, 26, 20, 15, 15, 0, 40, 60, 57, 44, 51, 75),
                          location = rep(c("A", "B"), each = 6))
    weather
```
    
Think: How might we modify the below density plot of `temp_3pm` to distinguish between locations?
    
```{r}
ggplot(weather, aes(x = temp_3pm)) +
      geom_density()
```    


c. `rain_today` (the response) and `location` (the predictor)

```{r}
weather <- data.frame(rain_today = c("no", "no", "no", "no", "yes", "no", "yes", "no", "yes", "yes", "no", "yes"),
                          location = c(rep("A", 7), rep("B", 5)))
    weather
```
    
Think: How might we modify the below bar plot of `location` to distinguish between days on which it did or didn't rain?
    
```{r}
ggplot(weather, aes(x = location)) +
      geom_bar()
```    






\
\
\
\





**General guidance for building bivariate plots**

As with univariate plots, an appropriate visualization for the relationship between 2 variables depends upon whether the variables are quantitative or categorical.
In *general*:

- Each **quantitative variable** requires a new *axis* (or a quantitative *scale* if we run out of axes).
- Each **categorical variable** requires a new way to "group" the graphic (eg: using colors, shapes, separate facets, etc)    
- For visualizations in which **overlap** in glyphs or plots obscures the patterns, try faceting or transparency.







\
\
\
\








## Exercises (required)


[Github user Tony McGovern](https://github.com/tonmcg/US_County_Level_Election_Results_08-20) has compiled and made available 2020/2016/2012 presidential election results for most of 3000+ U.S. counties, except Alaska. (Image: [Wikimedia Commons](https://commons.wikimedia.org/wiki/File:Map_of_USA_with_county_outlines_(black_%26_white).png))

![](https://upload.wikimedia.org/wikipedia/commons/6/68/Map_of_USA_with_county_outlines_%28black_%26_white%29.png)

A wrangled version of this data, is imported below, after being combined with:

- 2013 county-level demographics from the `df_county_demographics` data set within the `choroplethr` R package
- historical voting trends in the state in which the county falls. red = consistently Republican, blue = consistently Democratic, purple = something in between (from <https://www.270towin.com/content/blue-and-red-states>)

```{r}
# Load data
elections <- read.csv("https://mac-stat.github.io/data/election_2020_county.csv")

# Check it out
head(elections)
```


We'll use this data to explore voting outcomes within the U.S.'s 2-party system.
Here's a list of candidates by year:


year    Republican candidate    Democratic candidate
------- ----------------------- -------------------------
2020    Donald Trump            Joe Biden
2016    Donald Trump            Hillary Clinton
2012    Mitt Romney             Barack Obama







\
\
\
\




**Goals**

Use this data to explore how to build and interpret visualizations of **bivariate**, i.e. two variable, relationships.



\
\




**Directions**

- Go slow! Don't worry about memorizing the code, but *do* slow down to examine it, *even if* you're not asked explicitly about the code. Think about:        
    - How the data variables are used in the plot.
    - What `geom` the plot uses.
    - How color and fill work for the `geom` you're using.

- There are some coding "intuition checks" along the way. The goal here is to start connecting to and expanding upon the structure of our `ggplot()` code. Don't spend more than a few minutes trying these out. If you don't get it, no problem! You can move on and learn the code in the next exercise.

- The goal is for every student to finish the required exercises before we meet again.
    After that, choose your own adventure!        
    - If you're feeling like the required exercises gave you plenty to chew on, simply stop the activity. Go back and review the material you learned. *You never need to think about the optional exercises.*
    - If you're feeling like your brain still has some space after completing the required exercises, try the optional ones. These exercises are either "bonus" visualizations or intuition checks about some material coming up.
    



\
\
\
\



### Exercise 0: Review {-}

#### Part a {-}

How many, or roughly what percent, of the 3000+ counties did the Republican candidate win in 2020?

- Take a guess.
- Then make a plot of the `winner` variable.
- Then discuss what follow-up questions you might have (and that our data might help us answer).

```{r}

```





#### Part b {-}

The `repub_pct_20` variable provides more detail about the Republican support in each county.
Construct a plot of `repub_pct_20`.

Notice that the distribution of Republican support from county to county is slightly **left skewed** or **negatively skewed**.

What follow-up questions do you have?

```{r}


```








\
\
\
\




### Exercise 1: quantitative vs quantitative intuition check {-}

(Don't spend more than 3 minutes on this!)
Below is a **scatterplot** of the Republican support in 2020 vs 2016.
Notice that:

- both variables are quantitative, and get their own axes
- the *response* variable is on the y-axis, demonstrating how `repub_pct_20` might be predicted by `repub_pct_16`, not vice versa

Try to replicate this using `ggplot()`.
THINK:

- What info do you need to set up the canvas?
- What geometric layer (`geom_???`) might add these dots / points for each county? We haven't learned this yet, just take some guesses.


![](https://mac-stat.github.io/images/112/scatterplot_demo.png){width=50%}


```{r}

```





\
\
\
\





### Exercise 2: 2 quantitiative variables {-}

Run each chunk below to build up a a scatterplot of `repub_pct_20` vs `repub_pct_16` with different *glyphs* representing each county.
Address or think about any prompts in the comments (`#`).

```{r eval = FALSE}
# Set up the plotting frame
# How does this differ than the frame for our histogram of repub_pct_20 alone?
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16))
```

```{r eval = FALSE}
# Add a layer of points for each county
# Take note of the geom!
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_point()
```

```{r eval = FALSE}
# Change the shape of the points
# What happens if you change the shape to another number?
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_point(shape = 3)
```

```{r}
# YOU TRY: Modify the code to make the points "orange"
# NOTE: Try to anticipate if "color" or "fill" will be useful here. Then try both.
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_point()
```

```{r eval = FALSE}
# Add a layer that represents each county by the state it's in
# Take note of the geom and the info it needs to run!
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_text(aes(label = state_abbr))
```    


\
\
\
\




### Exercise 3: Reflect {-}

Summarize the relationship between the Republican support in 2020 and 2016.
Be sure to comment on:    

- the strength of the relationship (weak/moderate/strong)    
- the direction of the relationship (positive/negative)    
- outliers (in what state do counties deviate from the national trend? Any ideas why this might be the case?)    








\
\
\
\



### Exercise 4: Visualizing trend {-}

The trend of the relationship between `repub_pct_20` and `repub_pct_16` is clearly positive and (mostly) linear.
We can highlight this trend by adding a model "smooth" to the plot:    

```{r eval = FALSE}
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_point() +
  geom_smooth()
```


#### Part a {-}

Construct a new plot that contains the model smooth but does *not* include the individual point glyphs.

```{r}

```


#### Part b {-}

By default, `geom_smooth()` adds a smooth, localized model line.
To examine the "best" *linear model*, we can specify `method = "lm"`.
It's pretty similar in this example!
  
```{r eval = FALSE}
ggplot(elections, aes(y = repub_pct_20, x = repub_pct_16)) +
  geom_point() +
  geom_smooth(method = "lm")
```    










\
\
\
\


### Exercise 5: Your turn {-}

 
To examine how the 2020 results are related to some county demographics, construct scatterplots of `repub_pct_20` vs `median_rent`, and `repub_pct_20` vs `median_age`.
Summarize the relationship between these two variables and comment on which is the better predictor of `repub_pct_20`, `median_rent` or `median_age`.

```{r}
# Scatterplot of repub_pct_20 vs median_rent



# Scatterplot of repub_pct_20 vs median_age



```







\
\
\
\





### Exercise 6: A sad scatterplot {-}

Next, let's explore the relationship between a county's 2020 Republican support `repub_pct_20` and the `historical` political trends in its state.
In this case `repub_pct_20` is **quantitative**, but `historical` is **categorical**.
Explain why a scatterplot might *not* be an effective visualization for exploring this relationship.
(What questions does / doesn't it help answer?)

```{r}
ggplot(elections, aes(y = repub_pct_20, x = historical)) +
  geom_point()
```







\
\
\
\






### Exercise 7: quantitative vs categorical -- violins & boxes {-}

Though the above scatterplot *did* group the counties by `historical` category, it's nearly impossible to pick out meaningful patterns in 2020 Republican support in each category.
Let's try adding 2 different `geom` layers to the frame:

```{r eval = FALSE}
# Side-by-side violin plots
ggplot(elections, aes(y = repub_pct_20, x = historical)) +
  geom_violin()
```

```{r eval = FALSE}
# Side-by-side boxplots (defined below)
ggplot(elections, aes(y = repub_pct_20, x = historical)) +
  geom_boxplot()
```    


Box plots are constructed from five numbers - the minimum, 25th percentile, median, 75th percentile, and maximum value of a quantitative variable: 

![](https://mac-stat.github.io/images/112/boxplot.png){width=50%}


**REFLECT:**

Summarize what you've learned about the 2020 Republican county-level support within and between red/purple/blue states.  




\
\
\
\



### Exercise 8: quantitative vs categorical -- intuition check {-}

(Don't spend more than 3 minutes on this!)
We can also visualize the relationship between `repub_pct_20` and `historical` using our familiar density plots.
In the plot below, notice that we simply created a separate density plot for each `historical` category.
(The plot itself is "bad" but we'll fix it below.)
Try to adjust the code chunk below, which starts with a density plot of `repub_pct_20` alone, to re-create this image.


![](https://mac-stat.github.io/images/112/density_demo.png){width=50%}


```{r}
ggplot(elections, aes(x = repub_pct_20)) +
  geom_density()
```








\
\
\
\



### Exercise 9: quantitative vs categorical -- density plots {-}

Work through the chunks below and address the comments therein.

```{r eval = FALSE}
# Name two "bad" things about this plot
ggplot(elections, aes(x = repub_pct_20, fill = historical)) +
  geom_density()
```    

```{r eval = FALSE}
# What does scale_fill_manual do?
ggplot(elections, aes(x = repub_pct_20, fill = historical)) +
  geom_density() +
  scale_fill_manual(values = c("blue", "purple", "red"))
```    

```{r eval = FALSE}
# What does alpha = 0.5 do?
# Play around with different values of alpha, between 0 and 1
ggplot(elections, aes(x = repub_pct_20, fill = historical)) +
  geom_density(alpha = 0.5) +
  scale_fill_manual(values = c("blue", "purple", "red"))
```    

```{r eval = FALSE}
# What does facet_wrap do?!
ggplot(elections, aes(x = repub_pct_20, fill = historical)) +
  geom_density() +
  scale_fill_manual(values = c("blue", "purple", "red")) +
  facet_wrap(~ historical)
```

```{r eval = FALSE}
# Let's try a similar grouping strategy with a histogram instead of density plot.
# Why is this terrible?
ggplot(elections, aes(x = repub_pct_20, fill = historical)) +
  geom_histogram(color = "white") +
  scale_fill_manual(values = c("blue", "purple", "red"))
```





\
\
\
\







### Exercise 10 {-}

We've now learned 3 (of many) ways to visualize the relationship between a quantitative and categorical variable: side-by-side violins, boxplots, and density plots.

- Which do you like best?

- Name one pro of density plots relative to boxplots.

- Name one con of density plots relative to boxplots.








\
\
\
\







### Exercise 11: categorical vs categorical intuition check {-}

Finally, let's simply explore who *won* each county in 2020 (`winner_20`) and how this breaks down by `historical` voting trends in the state.
That is, let's explore the relationship between 2 categorical variables!
Following the same themes as above, we can utilize grouping features such as fill/color or facets to distinguish between different categories of `winner_20` and `historical`.

Spend at most 5 minutes on the following intuition check: adjust the code below to recreate the following two plots.


![](https://mac-stat.github.io/images/112/barplot_demo.png){width=60%}


```{r}
# Plot 1: adjust this to recreate the top plot
ggplot(elections, aes(x = historical)) +
  geom_bar()
```

```{r}
# Plot 2: adjust this to recreate the bottom plot
ggplot(elections, aes(x = winner_20)) +
  geom_bar()
```









\
\
\
\







### Exercise 12: categorical vs categorical {-}

Construct the following 4 bar plot visualizations.    

```{r eval = FALSE}
# A stacked bar plot
# How are the "historical" and "winner_20" variables mapped to the plot, i.e. what roles do they play?
ggplot(elections, aes(x = historical, fill = winner_20)) +
  geom_bar()
```

```{r eval = FALSE}
# A faceted bar plot
ggplot(elections, aes(x = winner_20)) +
  geom_bar() +
  facet_wrap(~ historical)
```    

```{r eval = FALSE}
# A side-by-side bar plot
# Note the new argument to geom_bar
ggplot(elections, aes(x = historical, fill = winner_20)) +
  geom_bar(position = "dodge")
```

```{r eval = FALSE}
# A proportional bar plot
# Note the new argument to geom_bar
ggplot(elections, aes(x = historical, fill = winner_20)) +
  geom_bar(position = "fill")
```


#### Part a {-}

Name one pro and one con of using the "proportional bar plot" instead of one of the other three options.    


#### Part b {-}

What's your favorite bar plot from part and why?








\
\
\
\




### Exercise 13: Practice (now or later) {-}

Decide what's best for you:

- Try this extra practice now.
- Reflect on the above exercises and come back to this extra practice later (but before the next class).

Import some daily weather data from a few locations in Australia:

```{r}
weather <- read.csv("https://mac-stat.github.io/data/weather_3_locations.csv")
```

Construct plots that address the research questions in each chunk.
You might make multiple plots (there are many ways to do things!), but don't just throw spaghetti at the wall.

Reflect before doing anything.
What types of variables are these?
How might you plot just 1 of the variables, and then tweak the plot to incorporate the other?

```{r}
# How do 3pm temperatures (temp3pm) differ by location?



```


```{r}
# How might we predict the 3pm temperature (temp3pm) by the 9am temperature (temp9am)?



```


```{r}
# How do the number of rainy days (raintoday) differ by location?


```













\
\
\
\




## Exercises (optional)

The above visualizations are foundational and important.
But they're not the only way to visualize the variables in our dataset.


### Optional exercise 1: Dealing with lots of categories  {-}

Suppose we wanted to better understand how the 2020 Republican support varied from county to county within each **state**.
Since `repub_pct_20` is quantitative and `state_abbr` is categorical, we *could* make a density plot of `repub_pct_20` for each state.
Reflect and why this is bad:

```{r}
ggplot(elections, aes(x = repub_pct_20, fill = state_abbr)) + 
  geom_density(alpha = 0.5)
```

A facet wrap would also be bad!

```{r fig.width = 8, fig.height = 8}
ggplot(elections, aes(x = repub_pct_20)) + 
  geom_density(alpha = 0.5) + 
  facet_wrap(~ state_abbr)
```


When we want to compare the distribution of some quantitative outcome among *many* groups / categories, a **ridgeline plot** can be a good option.
These are also called **joy plots**, named after the album cover for "Unknown Pleasures" by Joy Division.
(Look it up!)
To make a ridgeline plot, we can use the `geom_density_ridges()` function in the `ggridges` package.

```{r fig.height = 10, fig.width = 8, eval = FALSE}
# Install ggridges package
library(ggridges)

# Make our first joy plot
# THINK: What DON'T you like about this?
ggplot(elections, aes(x = repub_pct_20, y = state_abbr)) + 
  geom_density_ridges()
```

```{r fig.height = 10, fig.width = 8, eval = FALSE}
# Let's put the states in order by Republican support, not alphabet
# How do you think fct_reorder works? We'll learn about this later in the semester.
ggplot(elections, aes(x = repub_pct_20, y = fct_reorder(state_abbr, repub_pct_20))) + 
  geom_density_ridges(alpha = 0.5)
```

```{r fig.height = 10, fig.width = 8}
# YOUR TURN: color/fill the ridges according to a state's historical voting patterns 
# and add meaningful axis labels


```



**Follow-up questions**

- Which states tend to have the most variability in outcomes from county to county? The least?
- What other interesting patterns do you notice?
- Does this plot prompt any other questions?






\
\
\
\




### Optional exercise 2: total outcomes by state  {-}

Let's import some new data and counts up the total votes (Republican and Democratic) by state, not county.
This was **wrangled** from the elections data!

```{r}
elections_by_state <- read.csv("https://mac-stat.github.io/data/election_2020_by_state.csv")
```

For example, we might make a scatterplot of the 2020 vs 2016 outcomes:

```{r}
ggplot(elections_by_state, aes(y = repub_pct_20, x = repub_pct_16)) + 
  geom_point()
```


BUT this isn't the easiest way to communicate or identify the changes from 1 year to the next.


```{r fig.height = 10, fig.width = 8}
# YOU TRY
# Start by creating a "scatterplot" of state_abbr (y-axis) by 2020 Republican support on the x-axis
# Color the points red
# Scroll to solutions below when you're ready


```

```{r fig.height = 10, fig.width = 8}
# Check it out
ggplot(elections_by_state, aes(x = repub_pct_20, y = state_abbr)) + 
  geom_point(color = "red")
```

```{r fig.height = 10, fig.width = 8}
# YOU TRY
# Reorder the states in terms of their 2020 Republican support (not alphabet)
# Scroll to solutions below when you're ready


```

```{r fig.height = 10, fig.width = 8}
# Check it out
ggplot(elections_by_state, aes(x = repub_pct_20, y = fct_reorder(state_abbr, repub_pct_20))) + 
  geom_point(color = "red")
```

```{r fig.height = 10, fig.width = 8}
# Finally, add ANOTHER layer of points for the 2016 outcomes
# What info does this new geom_point() layer need to run?
ggplot(elections_by_state, aes(x = repub_pct_20, y = fct_reorder(state_abbr, repub_pct_20))) + 
  geom_point(color = "red") + 
  geom_point(aes(x = repub_pct_16, y = state_abbr))
```


**Reflect on the following**

- What do you think this plot needs? Try it! You might need to do some digging online.
- Summarize the main takeaways from the plots. Which states changed the most from 2016 to 2020? The least? Where did the Republican support increase? Where did it decrease?
- What other questions are you left with?









\
\
\
\




