.vlarge[How to make better spaghetti (plots) 🍝 🍝 ]

class: left, middle, inverse, title-slide

# .vlarge[How to make better spaghetti (plots) 🍝 🍝 ]
### .vlarge[Nicholas Tierney, Monash University]
### <div class="medium">
rstudioconf2020, January 29 <br><br> <svg viewBox='0 0 512 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M326.612 185.391c59.747 59.809 58.927 155.698.36 214.59-.11.12-.24.25-.36.37l-67.2 67.2c-59.27 59.27-155.699 59.262-214.96 0-59.27-59.26-59.27-155.7 0-214.96l37.106-37.106c9.84-9.84 26.786-3.3 27.294 10.606.648 17.722 3.826 35.527 9.69 52.721 1.986 5.822.567 12.262-3.783 16.612l-13.087 13.087c-28.026 28.026-28.905 73.66-1.155 101.96 28.024 28.579 74.086 28.749 102.325.51l67.2-67.19c28.191-28.191 28.073-73.757 0-101.83-3.701-3.694-7.429-6.564-10.341-8.569a16.037 16.037 0 0 1-6.947-12.606c-.396-10.567 3.348-21.456 11.698-29.806l21.054-21.055c5.521-5.521 14.182-6.199 20.584-1.731a152.482 152.482 0 0 1 20.522 17.197zM467.547 44.449c-59.261-59.262-155.69-59.27-214.96 0l-67.2 67.2c-.12.12-.25.25-.36.37-58.566 58.892-59.387 154.781.36 214.59a152.454 152.454 0 0 0 20.521 17.196c6.402 4.468 15.064 3.789 20.584-1.731l21.054-21.055c8.35-8.35 12.094-19.239 11.698-29.806a16.037 16.037 0 0 0-6.947-12.606c-2.912-2.005-6.64-4.875-10.341-8.569-28.073-28.073-28.191-73.639 0-101.83l67.2-67.19c28.239-28.239 74.3-28.069 102.325.51 27.75 28.3 26.872 73.934-1.155 101.96l-13.087 13.087c-4.35 4.35-5.769 10.79-3.783 16.612 5.864 17.194 9.042 34.999 9.69 52.721.509 13.906 17.454 20.446 27.294 10.606l37.106-37.106c59.271-59.259 59.271-155.699.001-214.959z'></path>
</svg> <strong>bit.ly/njt-rstudio</strong> <br><br> <svg viewBox='0 0 512 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M459.37 151.716c.325 4.548.325 9.097.325 13.645 0 138.72-105.583 298.558-298.558 298.558-59.452 0-114.68-17.219-161.137-47.106 8.447.974 16.568 1.299 25.34 1.299 49.055 0 94.213-16.568 130.274-44.832-46.132-.975-84.792-31.188-98.112-72.772 6.498.974 12.995 1.624 19.818 1.624 9.421 0 18.843-1.3 27.614-3.573-48.081-9.747-84.143-51.98-84.143-102.985v-1.299c13.969 7.797 30.214 12.67 47.431 13.319-28.264-18.843-46.781-51.005-46.781-87.391 0-19.492 5.197-37.36 14.294-52.954 51.655 63.675 129.3 105.258 216.365 109.807-1.624-7.797-2.599-15.918-2.599-24.04 0-57.828 46.782-104.934 104.934-104.934 30.213 0 57.502 12.67 76.67 33.137 23.715-4.548 46.456-13.32 66.599-25.34-7.798 24.366-24.366 44.833-46.132 57.827 21.117-2.273 41.584-8.122 60.426-16.243-14.292 20.791-32.161 39.308-52.628 54.253z'></path>
</svg> <strong>nj_tierney</strong>
</div>

---

layout: true
<div class="my-footer"><span>bit.ly/njt-rstudio • @nj_tierney</span></div>

---
class: middle, center, inverse

.huge[
How to make better spaghetti 🍝
]

---
class: middle, center, inverse

.huge[
~~How to make better spaghetti 🍝~~
]

.huge[
How to explore longitudinal data effectively
]

---
class: inverse, middle,

# What even **is** longitudinal data?

.huge[
> Individuals repeatedly measured through time
]

---

# What even **is** longitudinal data?

---

# What even **is** longitudinal data?

---

# What even **is** longitudinal data?

.large[
<table>
 <thead>
  <tr>
   <th style="text-align:left;"> country </th>
   <th style="text-align:right;"> year </th>
   <th style="text-align:right;"> height_cm </th>
  </tr>
 </thead>
<tbody>
  <tr>
   <td style="text-align:left;"> Australia </td>
   <td style="text-align:right;"> 1910 </td>
   <td style="text-align:right;"> 173 </td>
  </tr>
  <tr>
   <td style="text-align:left;"> Australia </td>
   <td style="text-align:right;"> 1920 </td>
   <td style="text-align:right;"> 173 </td>
  </tr>
  <tr>
   <td style="text-align:left;"> Australia </td>
   <td style="text-align:right;"> 1960 </td>
   <td style="text-align:right;"> 176 </td>
  </tr>
</tbody>
</table>
]

---

# What even **is** longitudinal data?

.large[
🎉 Individuals repeatedly measured through time
]

---

---

# All of Australia

---

# ...And New Zealand

---
# And the rest?

---
class: animated, heartBeat
# And the rest?

---
class: bg-black
background-image: url("gifs/noodle-explode.gif")
background-size: contain
background-position: 50% 50%
class: center, bottom, white

---

# Problems:

.large.pull-left[
<svg viewBox='0 0 512 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M12.41 148.02l232.94 105.67c6.8 3.09 14.49 3.09 21.29 0l232.94-105.67c16.55-7.51 16.55-32.52 0-40.03L266.65 2.31a25.607 25.607 0 0 0-21.29 0L12.41 107.98c-16.55 7.51-16.55 32.53 0 40.04zm487.18 88.28l-58.09-26.33-161.64 73.27c-7.56 3.43-15.59 5.17-23.86 5.17s-16.29-1.74-23.86-5.17L70.51 209.97l-58.1 26.33c-16.55 7.5-16.55 32.5 0 40l232.94 105.59c6.8 3.08 14.49 3.08 21.29 0L499.59 276.3c16.55-7.5 16.55-32.5 0-40zm0 127.8l-57.87-26.23-161.86 73.37c-7.56 3.43-15.59 5.17-23.86 5.17s-16.29-1.74-23.86-5.17L70.29 337.87 12.41 364.1c-16.55 7.5-16.55 32.5 0 40l232.94 105.59c6.8 3.08 14.49 3.08 21.29 0L499.59 404.1c16.55-7.5 16.55-32.5 0-40z'></path>
</svg> Overplotting

🙈 We don't see the individuals

🤷 Looking at **144** plots doesn't really help
]

.pull-right.animated.shake[
<img src="figures/gg-heights-heights-again-1.png" width="90%" style="display: block; margin: auto;" />
]

---

# Answers: Transparency?

---

# Answers: Transparency + a model?

???

- This helps reduce the overplotting
- It's not that this is wrong, it is useful - but we lose the individuals
- We only get the overall average. We dont get the rest of the information
- How do we even get started?
---

# But we forget about the **individuals**

???

- The model might make some good overall predictions
- But it can be really _ill suited_ for some individual
- Exploring this is somewhat clumsy - we need another way to explore

---
class: inverse, middle

# Three problems in exploring longitudinal data

.large[
Problem #1: How do I look at **some** of the data?
]

.large[
Problem #2: How do I find **interesting** observations?
]

.large[
Problem #3: How do I **understand** a model?
]

---

# Introducing `brolgar`: [brolgar.njtierney.com](https://brolgar.njtierney.com)

.pull-left.mlarge[
**br**owsing <br>
**o**ver <br>
**l**ongitudinal data <br>
**g**raphically, and <br>
**a**nalytically, in <br>
**r**
]

.pull-right[
<img src="imgs/brolga-bird.jpg" width="90%" style="display: block; margin: auto;" />
]

???

* It's a crane, it fishes, and it's a native Australian bird

---
class: inverse, middle

# What is longitudinal data?

.huge[
> Individuals repeatedly measured through time
]

---
class: inverse, middle,

# What is longitudinal data?

.huge[
> Individuals **repeatedly** measured through time
]

---
class: inverse, middle,

# What is longitudinal data?

.huge[
> Individuals **repeatedly** measured through **time**
]

---

class: inverse, middle

# 🤔 longitudinal data as a **time series**?

.vlarge[
> Anything that is observed sequentially over time **is a time series**
]

.large[
[-- Rob Hyndman and George Athanasopolous,
Forecasting: Principles and Practice](https://otexts.com/fpp2/data-methods.html)
]

---

# Longitudinal data as a time series <img src="https://tsibble.tidyverts.org/reference/figures/logo.png" align="right" height=140/>

```r
heights <- as_tsibble(heights,
                      index = year,
                      key = country,
*                     regular = FALSE)
```

.medium[
1. **index**: Your time variable
2. **key**: Variable(s) defining individual groups (or series)

`1. +  2.` determine distinct rows in a tsibble.

(From Dr. Earo Wang's talk: [Melt the clock](https://slides.earo.me/rstudioconf19/#8))
]

---

.large[

```
## # A tsibble: 1,490 x 3 [!]
## # Key:       country [144]
##   country      year height_cm
##   <chr>       <dbl>     <dbl>
## 1 Afghanistan  1870      168.
## 2 Afghanistan  1880      166.
## 3 Afghanistan  1930      167.
## 4 Afghanistan  1990      167.
## 5 Afghanistan  2000      161.
## 6 Albania      1880      170.
## # … with 1,484 more rows
```
]

---

# Longitudinal data as a time series <img src="https://tsibble.tidyverts.org/reference/figures/logo.png" align="right" height=140/>

.large[
We add information about **index** + **key**:

📐 Index = Year
  
  🔑   Key = Country
  ]
  
--

.large[
Record important time series information **once**

Use it **many** times in other places

]

---
class: inverse, center, middle
# Problem #1: How do I look at **some** of the data?

# How do you **eat** spaghetti?

---

.pull-left[
<img src="gifs/full-house-spag.gif" width="90%" style="display: block; margin: auto;" />
]

.pull-right[
<img src="gifs/garfield-spaghetti.gif" width="90%" style="display: block; margin: auto;" />
]

.animated.rollIn.huge.center[
data indigestion
]

---

.pull-left[

<img src="gifs/cat-single-spag.gif" width="100%" style="display: block; margin: auto;" />
]

.pull-right[
<img src="gifs/rat-spaghetti.gif" width="100%" style="display: block; margin: auto;" />
]

---

# Portion out your spaghetti! 🍝 🍝 🍝 🍝

---

# Look at **one** set of subsamples 🍝

---

# Look at **many** subsamples 🍝 🍝 🍝 🍝 🍝 🍝

---
class: inverse

# **How** do I look at many subsamples? 🤔

.vlarge[
How many keys are there?
]
--

.vlarge[
How many facets do I want?
]

.vlarge[
How many keys per facet?
]

---
class: inverse

# **How** do I look at many subsamples? 🤔

.vlarge[
How to keep the same number of keys per plot?
]

.vlarge[
What is `rep`, `rep.int`, and `rep_len`?
]

.vlarge[
Do I want `length.out` or `times`?
]

---
background-image: url("imgs/what-year-is-it.jpg")
background-size: contain
background-position: 50% 50%
class: center, bottom, white

---

# Distraction threshold ⏲ 🐰 🕳

.large[

(Something I made up)
]

.large[
> If solving a problem requires solving **3+ smaller problems**

> Your focus shifts from the current goal to something else.

> You are distracted.

]

???

- Task one

- Task one being overshadowed slightly by minor task 1
- Task one being overshadowed slightly by minor task 2
- Task one being overshadowed slightly by minor task 3

---

# Avoiding the rabbit hole

.large[
We can blame ourselves when we are distracted for **not being better**.
]

.large[
It's not that we should be better, rather **with better tools we could be more efficient**. 
]

.large[
We need to make things **as easy as reasonable**, with the least amount of distraction.
]

---

# `facet_sample()`: See more individuals 👀

.large[

```r
ggplot(heights, aes(x = year, 
                    y = height_cm, 
                    group = country)) + 
  geom_line() 
```
]

---

# `facet_sample()`: See more individuals 👀

# `facet_sample()`: See more individuals 👀
---

# `facet_sample()`: See more individuals 👀

.large[

```r
ggplot(heights, aes(x = year,
                    y = height_cm,
                    group = country)) + 
  geom_line() + 
* facet_sample()
```
]

---

# `facet_sample()`: See more individuals 👀

---

# Remove distraction. Ask relevant questions

.large[
> How many keys per facet?

> How many plots do I want to look at?
]

.large[

```r
gg_heights +  
  facet_sample(
*   n_per_facet = 3,
*   n_facets = 9
  )
```
]

---

---

# How to see **all** individuals?

.huge[

```r
gg_heights +  
* facet_strata()
```
]

---

# `facet_strata()`: See **all** individuals

---

# 🤔 ... can we re-order these facets?

???

In asking these questions we can solve something else interesting

---

# We can re-order these facets?! 😄

---

# See all individuals **along** some variable

.huge[

```r
gg_heights + 
  facet_strata(
*   along = -year
    ) 
```
]

---

# See all individuals **along** some variable

---

# Magic facets: Focus on relevant questions instead of minutiae:

.pull-left[

**`facet_sample()`**

"How many lines per facet"

"How many facets?"

```r
gg_heights +  
  facet_sample(
*   n_per_facet = 10,
*   n_facets = 12
    ) 
```

]

.pull-right[

**`facet_strata()`**

"How many facets / strata?"

"**What** to arrange plots **along**?"

```r
gg_heights +  
  facet_strata(
*   n_strata = 10,
*   along = -year
    ) 
```

]

---

# `facet_strata()` & `facet_sample()`

.pull-left[
<img src="imgs/car-hood-small.jpg" width="100%" style="display: block; margin: auto;" />
]

.pull-right.large[
<svg viewBox='0 0 320 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M296 160H180.6l42.6-129.8C227.2 15 215.7 0 200 0H56C44 0 33.8 8.9 32.2 20.8l-32 240C-1.7 275.2 9.5 288 24 288h118.7L96.6 482.5c-3.6 15.2 8 29.5 23.3 29.5 8.4 0 16.4-4.4 20.8-12l176-304c9.3-15.9-2.2-36-20.7-36z'></path>
</svg> with `sample_n_keys()` & `stratify_keys()`

You can still get at data and do manipulations
]

---

# Problem #1: How do I look at some of the data?

.left-code.large[

`as_tsibble()`

`facet_sample()`

`facet_strata()`

]

.right-plot.large[

Store useful information

View **many** subsamples

View **all** subsamples

]

---

# ~~Problem #1: How do I look at some of the data?~~

.left-code.large[

`as_tsibble()`

`facet_sample()`

`facet_strata()`

]

.right-plot.large[

Store useful information

View **many** subsamples

View **all** subsamples

]

---

## Problem #2: How do I find **interesting** observations?

---
class: inverse, center, middle

.huge[
A workflow
]

.huge[
**Define** what is interesting:
]

.huge.animated.bounceInLeft[
_maximum height_
]

---

## Identify features: one observation per **key**

---

## Identify features: one observation per **key**

---

## Identify features: one observation per **key**

---

## Identify **important features** and **decide** how to filter

---

## Identify **important features** and **decide** how to filter

---

## Join this feature back to the data

---

## Join this feature back to the data

---

## 🎉 Countries with smallest and largest max height

---
class: inverse, middle

.vhuge[
Let's see that **one more time**, but with the data
]

---

## Identify features: one observation per **key**

.large[

---

## Identify features: one observation per **key**

.large[

```
## # A tibble: 144 x 2
##   country       max
##   <chr>       <dbl>
## 1 Afghanistan  168.
## 2 Albania      170.
## 3 Algeria      171.
## 4 Angola       169.
## 5 Argentina    174.
## 6 Armenia      172.
## # … with 138 more rows
```
]

---

## Identify important features and decide how to filter

```r
heights_five %>% 
  filter(max == max(max) | max == min(max))
```

```
## # A tibble: 2 x 2
##   country            max
##   <chr>            <dbl>
## 1 Denmark           183.
## 2 Papua New Guinea  161.
```

---

## Join summaries back to data

```r
heights_five %>% 
  filter(max == max(max) | max == min(max)) %>% 
  left_join(heights, by = "country")
```

```
## # A tibble: 21 x 4
##   country   max  year height_cm
##   <chr>   <dbl> <dbl>     <dbl>
## 1 Denmark  183.  1820      167.
## 2 Denmark  183.  1830      165.
## 3 Denmark  183.  1850      167.
## 4 Denmark  183.  1860      168.
## 5 Denmark  183.  1870      168.
## 6 Denmark  183.  1880      170.
## # … with 15 more rows
```

---

background-image: url("gifs/dog-solve-problem.gif")
background-size: contain
background-position: 50% 50%
class: center, bottom, white

---

.vhuge[
But Nick, how did you create those summaries?
]

---

## Identify features: one per **key** <img src="https://feasts.tidyverts.org/reference/figures/logo.png" align="right" height=140/>

```r
heights %>%
* features(height_cm, feat_five_num)
```

```
## # A tibble: 144 x 6
##   country       min   q25   med   q75   max
##   <chr>       <dbl> <dbl> <dbl> <dbl> <dbl>
## 1 Afghanistan  161.  164.  167.  168.  168.
## 2 Albania      168.  168.  170.  170.  170.
## 3 Algeria      166.  168.  169   170.  171.
## 4 Angola       159.  160.  167.  168.  169.
## 5 Argentina    167.  168.  168.  170.  174.
## 6 Armenia      164.  166.  169.  172.  172.
## # … with 138 more rows
```

---

## What is the range of the data? `feat_ranges`

```r
heights %>%
  features(height_cm, feat_ranges)
```

```
## # A tibble: 144 x 5
##   country       min   max range_diff   iqr
##   <chr>       <dbl> <dbl>      <dbl> <dbl>
## 1 Afghanistan  161.  168.       7     3.27
## 2 Albania      168.  170.       2.20  1.53
## 3 Algeria      166.  171.       5.06  2.15
## 4 Angola       159.  169.      10.5   7.87
## 5 Argentina    167.  174.       7     2.21
## 6 Armenia      164.  172.       8.82  5.30
## # … with 138 more rows
```

---

## Does it only increase or decrease? `feat_monotonic`

```r
heights %>%
  features(height_cm, feat_monotonic)
```

```
## # A tibble: 144 x 5
##   country     increase decrease unvary monotonic
##   <chr>       <lgl>    <lgl>    <lgl>  <lgl>    
## 1 Afghanistan FALSE    FALSE    FALSE  FALSE    
## 2 Albania     FALSE    TRUE     FALSE  TRUE     
## 3 Algeria     FALSE    FALSE    FALSE  FALSE    
## 4 Angola      FALSE    FALSE    FALSE  FALSE    
## 5 Argentina   FALSE    FALSE    FALSE  FALSE    
## 6 Armenia     FALSE    FALSE    FALSE  FALSE    
## # … with 138 more rows
```

---

## What is the spread of my data? `feat_spread`

```r
heights %>%
  features(height_cm, feat_spread) 
```

```
## # A tibble: 144 x 5
##   country        var    sd   mad   iqr
##   <chr>        <dbl> <dbl> <dbl> <dbl>
## 1 Afghanistan  7.20  2.68  1.65   3.27
## 2 Albania      0.950 0.975 0.667  1.53
## 3 Algeria      3.30  1.82  0.741  2.15
## 4 Angola      16.9   4.12  3.11   7.87
## 5 Argentina    2.89  1.70  1.36   2.21
## 6 Armenia     10.6   3.26  3.60   5.30
## # … with 138 more rows
```

---

## features: MANY more features in `feasts` <img src="https://feasts.tidyverts.org/reference/figures/logo.png" align="right" height=140/>

.mlarge[

`feat_acf`: autocorrelation-based features

`feat_stl`: STL (Seasonal, Trend, and Remainder by LOESS) decomposition

[Create your own features](http://brolgar.njtierney.com/articles/finding-features.html#creating-your-own-features)

]

---

# ~~Problem #1: How do I look at **some** of the data?~~

# Problem #2: How do I find **interesting** observations?

---

# ~~Problem #1: How do I look at **some** of the data?~~

# ~~Problem #2: How do I find **interesting** observations?~~

# Problem #3: How do I **understand** a model?

.vhuge.center[
🤔
]

---

# Problem #3: How do I **understand** a model?

.medium[
Let's fit a mixed effects model.

Fixed effect of year + Random intercept for country
]

```r
heights_fit <- lmer(height_cm ~ year + (1|country), heights)

heights_aug <- heights %>%
  add_predictions(heights_fit, var = "pred") %>%
  add_residuals(heights_fit, var = "res")
```

---

# Problem #3: How do I **understand** a model?

.large[

```
## # A tsibble: 1,490 x 5 [!]
## # Key:       country [144]
##   country      year height_cm  pred    res
##   <chr>       <dbl>     <dbl> <dbl>  <dbl>
## 1 Afghanistan  1870      168.  164.  4.59 
## 2 Afghanistan  1880      166.  164.  1.52 
## 3 Afghanistan  1930      167.  166.  0.823
## 4 Afghanistan  1990      167.  168. -1.04 
## 5 Afghanistan  2000      161.  169. -7.10 
## 6 Albania      1880      170.  168.  2.39 
## # … with 1,484 more rows
```
]

---

# Problem #3: How do I **understand** a model?

---

# Look at **many** subsamples? `facet_sample()`

```r
gg_heights_fit + facet_sample()
```

---

# Look at **all** subsamples? `facet_strata()`

```r
gg_heights_fit + facet_strata()
```

---

# Look at **all** subsamples **along** residuals?

```r
gg_heights_fit + facet_strata(along = -res)
```

---

# Look at the predictions with the data?

```r
set.seed(2020-01-21)
heights_sample <- 
heights_aug %>%
* sample_n_keys(size = 9) %>%
  ggplot(aes(x = year,
             y = pred,
             group = country)) + 
  geom_line() +
  facet_wrap(~country)

heights_sample
```

---

# Look at the predictions with the data?

---

# Look at the predictions with the data?

```r
heights_sample + geom_point(aes(y = height_cm))
```

---

# Take homes

## Problem #1: How do I look at **some** of the data?

.mlarge[
1. Longitudinal data as a time series 💹
2. Specify structure, get a free lunch. 🥪
3. Look at as much of the raw data as possible  🍣
4. Use `facet_sample()` / `facet_strata()`
]

---

# Take homes

## Problem #2: How do I find **interesting** observations?

.mlarge[
1. Decide what features are interesting
2. Summarise down to one observation per key
3. Decide how to filter
4. Join this feature back to the data
]

---

# Take homes

## Problem #3: How do I **understand** a model?

.mlarge[
1. Look at (one, more or all!) subsamples
1. Arrange subsamples
1. (actually use similar approaches to earlier!)
]

---

# Future Directions

.large[
<svg viewBox='0 0 512 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M501.54 60.91c17.22-17.22 12.51-46.25-9.27-57.14a35.696 35.696 0 0 0-37.37 3.37L251.09 160h151.37l99.08-99.09zM496 192H16c-8.84 0-16 7.16-16 16v32c0 8.84 7.16 16 16 16h16c0 80.98 50.2 150.11 121.13 178.32-12.76 16.87-21.72 36.8-24.95 58.69-1.46 9.92 6.04 18.98 16.07 18.98h223.5c10.03 0 17.53-9.06 16.07-18.98-3.22-21.89-12.18-41.82-24.95-58.69C429.8 406.11 480 336.98 480 256h16c8.84 0 16-7.16 16-16v-32c0-8.84-7.16-16-16-16z'></path>
</svg> More features (summaries)

---

# Thanks

.large.pull-left[
- Di Cook
- Tania Prvan
- Stuart Lee
- Mitchell O'Hara Wild
- Earo Wang
]

.large.pull-right[
- Rob Hyndman
- Miles McBain
- Hadley Wickham
- Garrick Aden-Buie
- Monash University
- ACEMS
]

---

# Resources

.large[
<svg viewBox='0 0 640 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M425.7 256c-16.9 0-32.8-9-41.4-23.4L320 126l-64.2 106.6c-8.7 14.5-24.6 23.5-41.5 23.5-4.5 0-9-.6-13.3-1.9L64 215v178c0 14.7 10 27.5 24.2 31l216.2 54.1c10.2 2.5 20.9 2.5 31 0L551.8 424c14.2-3.6 24.2-16.4 24.2-31V215l-137 39.1c-4.3 1.3-8.8 1.9-13.3 1.9zm212.6-112.2L586.8 41c-3.1-6.2-9.8-9.8-16.7-8.9L320 64l91.7 152.1c3.8 6.3 11.4 9.3 18.5 7.3l197.9-56.5c9.9-2.9 14.7-13.9 10.2-23.1zM53.2 41L1.7 143.8c-4.6 9.2.3 20.2 10.1 23l197.9 56.5c7.1 2 14.7-1 18.5-7.3L320 64 69.8 32.1c-6.9-.8-13.5 2.7-16.6 8.9z'></path>
</svg> [feasts](http://feasts.tidyverts.org/)

---

# Colophon

.medium[
<svg viewBox='0 0 512 512' xmlns='http://www.w3.org/2000/svg' style='height:1em;fill:currentColor;position:relative;display:inline-block;top:.1em;'>
  <path d='M12.41 148.02l232.94 105.67c6.8 3.09 14.49 3.09 21.29 0l232.94-105.67c16.55-7.51 16.55-32.52 0-40.03L266.65 2.31a25.607 25.607 0 0 0-21.29 0L12.41 107.98c-16.55 7.51-16.55 32.53 0 40.04zm487.18 88.28l-58.09-26.33-161.64 73.27c-7.56 3.43-15.59 5.17-23.86 5.17s-16.29-1.74-23.86-5.17L70.51 209.97l-58.1 26.33c-16.55 7.5-16.55 32.5 0 40l232.94 105.59c6.8 3.08 14.49 3.08 21.29 0L499.59 276.3c16.55-7.5 16.55-32.5 0-40zm0 127.8l-57.87-26.23-161.86 73.37c-7.56 3.43-15.59 5.17-23.86 5.17s-16.29-1.74-23.86-5.17L70.29 337.87 12.41 364.1c-16.55 7.5-16.55 32.5 0 40l232.94 105.59c6.8 3.08 14.49 3.08 21.29 0L499.59 404.1c16.55-7.5 16.55-32.5 0-40z'></path>
</svg> Slides made using [xaringan](https://github.com/yihui/xaringan)

---

# Learning more

]

---
class: animated, heartBeat

.vhuge[
BONUS ROUND
]

---

---

---

# What about interactive graphics?

.vlarge[
> if you make a tooltip or rollover, **assume no one will ever see it**"
-- Archie Tse, NYT
]

---

# What if we grabbed a sample of those who have the best, middle, and worst residuals?

.medium[

```r
summary(heights_aug$res)
```

```
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## -8.1707 -1.6202 -0.1558  0.0000  1.3545 12.1729
```
]

.mlarge[
Which countries are **nearest** to these statistics?
]

---

# use `keys_near()`

.pull-left[

```r
keys_near(heights_aug, 
          var = res)
```

```
## # A tibble: 6 x 5
##   country       res stat  stat_value stat_diff
##   <chr>       <dbl> <fct>      <dbl>     <dbl>
## 1 Ireland    -8.17  min       -8.17   0       
## 2 Azerbaijan -1.62  q_25      -1.62   0.000269
## 3 Laos       -0.157 med       -0.156  0.00125 
## 4 Mongolia   -0.155 med       -0.156  0.00125 
## 5 Egypt       1.35  q_75       1.35   0.000302
## 6 Poland     12.2   max       12.2    0
```
]

.pull-right[

> 🔑 🔑 🔑 that best match the 5 number summary.

]

---

# Join data back and explore

---

.vhuge.left-code[
**End.**
]

.right-plot[
<img src="imgs/explorer.jpg" width="100%" style="display: block; margin: auto;" />

]