R/cal.prior.R

In bayclumpr: Bayesian Analysis of Clumped Isotope Datasets

#' Generate a dataset reflecting the priors used to run the analyses
#'
#' @param prior Informative or not
#' @param n number of observations to simulate
#'
#' @return A \code{data.frame} with prior distributions.
#'
#' @importFrom stats rnorm
#'
#' @export

cal.prior <- function(prior, n = 1000) {
  if (prior == "Informative") {
    params <- cbind.data.frame(
      parameter = c("alpha", "beta"),
      mean = c(0.231, 0.039),
      sd = c(0.065, 0.004)
    )
    params
  } else {
    params <- cbind.data.frame(
      parameter = c("alpha", "beta"),
      mean = c(0, 0.01),
      sd = c(0, 0.01)
    )
    params
  }

  data <- cbind.data.frame(
    alpha = rnorm(n, params[1, 2], params[1, 3]),
    beta = rnorm(n, params[2, 2], params[2, 3])
  )
  attr(data, "priors") <- prior
  attr(data, "params") <- params
  data
}