Easy Beta distribution plots in base R - auto scaling axes

beta_plot.R

beta_plot <- function(shape_1 = 1, shape_2 = 1, 
                      x_min = qbeta(0.0001, shape_1, shape_2), 
                      x_max = qbeta(0.9999, shape_1, shape_2)){
  x <- seq(x_min, x_max, length.out = 1000)
  y <- dbeta(x, shape_1, shape_2)
  plot(x, y, 
       type = "l", 
       ylab = "density", 
       xlab = "probability", 
       main = paste0("Beta(", shape_1, ", ", shape_2, ")"))
}