R/many_normal_plots.R

In TeachBayes: Teaching Bayesian Inference

many_normal_plots <- function(list_normal_par){
  N <- length(list_normal_par)
  for(j in 1:N)
    list_normal_par[[j]] <- unlist(list_normal_par[[j]])
  Means <- sapply(list_normal_par, function(y) y[1])
  SDs <- sapply(list_normal_par, function(y) y[2])
  Means <- round(Means, 1)
  SDs <- round(SDs, 1)
  gmin <- min(Means - 4 * SDs)
  gmax <- max(Means + 4 * SDs)
  labels <- paste("N(", Means, ", ", SDs, ")", sep="")
  x <- seq(gmin, gmax, length.out=200)
  df <- NULL
  for (j in 1:N){
    y <- x
    f <- dnorm(x, Means[j], SDs[j])
    Model <- labels[j]
    df <- rbind(df, data.frame(Model, y, f))
  }
  p <- ggplot(df, aes(y, f, group=Model, color=Model)) +
    geom_line(size=1.5)
  if(N == 2){
    p <- p +
      scale_colour_manual(values = c("blue", "red"))
  }
  p
}