Summarising Anscombe’s quartet

Tmisc::quartet |>
  group_by(set) |>
  summarise(
    mean_x = mean(x), 
    mean_y = mean(y),
    sd_x = sd(x),
    sd_y = sd(y),
    r = cor(x, y)
  )|> 
  gt::gt() # print as table

Visualizing Anscombe’s quartet

ggplot2 \(\in\) tidyverse

• ggplot2 is tidyverse’s data visualization package.

• gg in “ggplot2” stands for Grammar of Graphics.

• Inspired by the book Grammar of Graphics by Leland Wilkinson.

Grammar of Graphics

A grammar of graphics is a tool that enables us to concisely describe the components of a graphic.

Source: BloggoType

Hello ggplot2!

ggplot2: Elegant Graphics for Data Analysis

Simple example

ggplot(data = mpg,                               # Use the mpg data
       mapping = aes(x = displ, y = hwy)) +      # map (explore) the variables displ and hwy
  geom_point()                                   # visualize the observations with points

colour, size, shape and other aesthetic attributes

ggplot(mpg, aes(displ, hwy, colour = drv))+ 
  geom_point()

ggplot(mpg, aes(displ, hwy, shape = drv))+ 
  geom_point()

ggplot(mpg, aes(displ, hwy, size = drv))+ 
  geom_point()

Faceting

facet_wrap()

ggplot(mpg, aes(displ, hwy)) + 
  geom_point() + 
  facet_wrap(~class)

Other Plot Geoms - geom_smooth()

An important argument to geom_smooth() is the method, which allows you to choose which type of model is used to fit the smooth curve.
Default is method = "loess"

ggplot(mpg, aes(displ, hwy)) + 
  geom_point() + 
  geom_smooth(method = "loess")

ggplot(mpg, aes(displ, hwy)) + 
  geom_point() + 
  geom_smooth(method = "lm") #linear model

Other Plot Geoms

ggplot(mpg, aes(drv, hwy))+
  geom_jitter()

ggplot(mpg, aes(drv, hwy))+
  geom_boxplot()

ggplot(mpg, aes(drv, hwy))+
  geom_violin()

Other Plot Geoms

geom_histogram() geom_freqploy()

ggplot(mpg, aes(hwy))+
  geom_histogram()

ggplot(mpg, aes(hwy))+
  geom_freqpoly()

Other Plot Geoms

geom_histogram() geom_freqploy()

Both geoms bin the data. You should experiment with the binwidth or bin arguments when exploring your data.

ggplot(mpg, aes(hwy)) + 
  geom_freqpoly(binwidth = 10)

ggplot(mpg, aes(hwy)) + 
  geom_freqpoly(binwidth = 1)

Other Plot Geoms geom_bar()

ggplot(mpg, aes(x = manufacturer)) + 
  geom_bar()

ggplot(mpg, aes(y =manufacturer)) + 
  geom_bar()

Other Plot Geoms

geom_bar() + fct_infreq()

mpg |> 
ggplot(aes(y = fct_infreq(manufacturer))) + 
  geom_bar()

Other Plot Geoms

geom_bar() + fct_infreq() and aes(fill = [VAR] )

mpg |> 
ggplot(aes(y = fct_infreq(manufacturer), 
           fill = class)) + 
  geom_bar()

Other Plot Geoms

geom_bar() and position = "fill"

position = "fill" works like stacking, but makes each set of stacked bars the same height. This makes it easier to compare proportions across groups.

mpg |> 
ggplot(aes(y = fct_infreq(manufacturer), 
           fill = class)) + 
  geom_bar(position = "fill")

Other Plot Geoms

geom_bar() and position = "dodge"

position = "dodge" places overlapping objects directly beside one another. This makes it easier to compare individual values.

mpg |> 
ggplot(aes(y = fct_infreq(manufacturer), 
           fill = class)) + 
  geom_bar(position = "dodge")

labs()- Title, Subtitle, Legends, and more

labs() can do it all

mpg |> 
  ggplot(aes(x = class, y = cty, fill = class))+
  geom_boxplot()+
  labs(
    title = "Milage in cities of various car classes",
    subtitle = "Car classes are colored with nice colors using the fill argument!!",
    x = "Class of car",
    y = "Milage in cities",
    fill = "Class",
    caption = "The data comes from the mpg dataset"
    )

Themes

ggplot comes prepackaged with several different themes

p <- mpg |> 
  ggplot(aes(x = class, y = cty, fill = class))+
  geom_boxplot()

p + theme_bw()

p + theme_classic()

Themes

The package ggthemes have many different themes.

p + ggthemes::theme_stata()
p + ggthemes::theme_economist()
p + ggthemes::theme_excel_new()
p + ggthemes::theme_fivethirtyeight()

color and fill scales

Notice fill is mapped to a categorical variable

p <- mpg |> 
  ggplot(aes(x = cty, y = hwy, fill=drv))+
  geom_point(shape = 21, size = 3)

p

color and fill scales

scale_fill_brewer(type = )

p + scale_fill_brewer(type = "seq")  # sequential (default)

p + scale_fill_brewer(type = "qual") # qualitative

p + scale_fill_brewer(type = "div")  # diverging

color and fill scales

scale_fill_brewer(palette = )

p + scale_fill_brewer(type = "seq", # sequential (default)
                      palette = 2)  

p + scale_fill_brewer(type = "qual", # qualitative
                      palette = 2) 

p + scale_fill_brewer(type = "div", # diverging
                      palette = 2)  

color and fill scales

Notice color is mapped to a continous variable

p <- mpg |> 
  ggplot(aes(x = cty, y = hwy, color=displ))+
  geom_point(size = 3)

p

color and fill scales

scale_colour_continuous(type = )

p + 
  scale_colour_continuous(type = "gradient")

# The virdidis color scale
p + 
  scale_colour_continuous(type = "viridis")

p + 
  scale_colour_continuous(type = "viridis", 
                          option = "magma")

The patchwork package

Arranging plots with the patchwork package

p1 <- ggplot(mpg) + 
  geom_point(aes(x = displ, y = hwy))+
  ggtitle("p1")

p2 <- ggplot(mpg) + 
  geom_bar(aes(x = as.character(year), fill = drv), position = "dodge") + 
  labs(x = "year")+
  ggtitle("p2")

p3 <- ggplot(mpg) + 
  geom_density(aes(x = hwy, fill = drv), colour = NA) + 
  facet_grid(rows = vars(drv))+
  ggtitle("p3")

p4 <- ggplot(mpg) + 
  stat_summary(aes(x = drv, y = hwy, fill = drv), geom = "col", fun.data = mean_se) +
  stat_summary(aes(x = drv, y = hwy), geom = "errorbar", fun.data = mean_se, width = 0.5)+
  ggtitle("p4")

Arranging plots with the patchwork package

Arranging plots with the patchwork package

• + does not specify any specific layout, only that the plots should be displayed together.

• patchwork automagically decides the layout.

• Adding 3 plots together will create a 1x3 grid.

• Adding 4 plots together will create a 2x2 grid.

p1 + p2

Arranging plots with the patchwork package

• + does not specify any specific layout, only that the plots should be displayed together.

• patchwork automagically decides the layout.

• Adding 3 plots together will create a 1x3 grid.

• Adding 4 plots together will create a 2x2 grid.

p1 + p2 + p3 +p4

Arranging plots with the patchwork package

plot_layout()

You have a very high degree of control via the plot_layout() function

p1 + p2 + p3 + p4 + plot_layout(nrow = 1)

Arranging plots with the patchwork package

/ and |

/ forces a single column.
| forces a single row.

p1/p2

p1|p2 # sames as p1 + p2

Arranging plots with the patchwork package

nesting plots

p3 | (p2 /(p1 | p4))

Arranging plots with the patchwork package

guides = "collect"

p1 + p2 + p3 + plot_layout(ncol = 2, guides = "collect")

Arranging plots with the patchwork package

guide_area()

p1 + p2 + p3 + plot_layout(ncol = 2, guides = "collect") + guide_area() 

Arranging plots with the patchwork package

patchwork comes with the powerful & operator that allows you change the settings on all the individual plots.

p1 + p2 + p3 + 
  plot_layout(ncol = 2, guides = "collect") + 
  guide_area() +
  scale_fill_brewer(type = "qual",
                    palette = 3) +
  theme_classic()

Can you spot the &s in the code below?

p1 + p2 + p3 + 
  plot_layout(ncol = 2, guides = "collect") + 
  guide_area() &   
  scale_fill_brewer(type = "qual",
                    palette = 3) &
  theme_classic()