Simple function for simulating/plotting beta distributions for Dataverse 2017 Package Development Workshop

beta_plot_basic.R

#' @import ggplot2
#' @export

beta_plot <- function(n = 10000, a = 1, b = 3) {
    # draw distributions
    sims <- rbeta(n = n, shape1 = a, shape2 = b)

    # convert to data frame for ggplot2 compatability
    sims <- data.frame(x = sims)

    # plot probability density function
    ggplot(sims, aes(x)) +
        geom_density() +
        xlab("") + ylab("Probability Density Function") +
        theme_bw()
}