beta_prior_post <- function(prior_shapes, post_shapes,
label_1="Prior",
label_2="Posterior"){
TH <- theme(
plot.title = element_text(
colour = "blue",
size = 18,
hjust = 0.5,
vjust = 0.8,
angle = 0
)
)
x <- NULL
ggplot(data.frame(x=c(0, 1)), aes(x)) +
stat_function(fun=dbeta, geom="line",
aes(linetype=label_1), size=1.5,
args=list(shape1=prior_shapes[1],
shape2=prior_shapes[2])) +
stat_function(fun=dbeta, geom="line",
aes(linetype=label_2), size=1.5,
args=list(shape1=post_shapes[1],
shape2=post_shapes[2])) +
scale_colour_manual(values=c("red", "blue")) +
labs(colour = "Type") +
xlab("P") + ylab("Density")
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.