1 Load packages 2 Simulate data 3 Plot 1 4 Plot 2 5 Reproducibility 1 Load packages library(tidyverse) # data wrangling 2 Simulate data low_spread <- tibble(var = rnorm(n = 100), id = 1:100, type = "low spread") high_spread <- tibble(var= rnorm(n = 100, sd = 10), id = 1:100, type = "high spread") d <- low_spread %>% bind_rows(high_spread) 3 Plot 1 ggplot(d) + aes(x = id, y = var) + facet_wrap(~ type) + geom_hline(yintercept = 0, color = "grey40") + geom_point() + theme_minimal() 4 Plot 2 ggplot(d) + aes(x = type, y = var) + geom_boxplot() 5 …