library(tidyverse)
tibble(x = seq(-3, 3, by = 0.1)) %>%
ggplot(aes(x)) +
stat_function(fun = pnorm, colour = "steelblue3",
args = list(mean = 0, sd = 1),
size = 1.5) +
labs(y = "Probability", x = "Z Score") +
scale_y_continuous(limits = c(0, 1)) +
scale_x_continuous(limits = c(-3.5, 3.5), breaks = -3:3) +
ggtitle("Cumulative distribution function / pnorm()") +
theme_linedraw()