set.seed(2023)
differences <- replicate(5000, mean(rnorm(30, 100, 15)) - mean(rnorm(30, 100, 15)))
ggplot(data.frame(differences), aes(differences)) +
geom_histogram(bins = 40, color = "white") +
geom_vline(xintercept = 0) +
labs(title = "Differences of Means of Pairs of Samples",
x = "Xbar1 - Xbar2", y = "Frequency")
