2ksamp.R

n_ <- 2000
m_ <- 10000

g <- ggplot(roll(m = m_, n = n_),
            mapping = aes(x = V1)) +
  geom_vline(xintercept = 3.5, colour = "tomato3") +
  labs(
    subtitle = str_interp('Density of means of dice 
                          rolls for ${n_} dice over ${m_} rolls.'),
    x = 'Mean Score',
    y = 'Proportion of the time the value was observed'
  ) +
  geom_bar(aes(y = ..prop..), alpha = 0.4) +
  scale_x_continuous(
    breaks = 1:6,
    lim =  c(0, 7)
  ) +
  geom_density()

g