R/binom.prior.R

In StudyPrior: Distributions for Historical Information

#' Prior distribution for binomial data
#'
#' @param x Vector of historical events
#' @param n Vector of historical sample sizes
#' @param type Which prior to provide
#' @param ... Parameters to provide to prior
#' @param verbose Print messages
#'
#' @return A density function
#' @export
#'

binom.prior <- function(type = c("MAP.FB", "MAP.EB","PP.FB","PP.EB","PP.EB.Beta","Beta.EB","PP.EB.Sep", "PP.Fix","PP.Cor"),
                        x, n,
                        verbose=FALSE,
                        ...){

  type <- match.arg(type)

  f.to.call <- switch(type,
                      "MAP.FB" = binom.MAP.FB,
                      "MAP.EB" = binom.MAP.EB,
                      "PP.FB" = binom.PP.FB.MC.BE,
                      "PP.EB" = binom.PP.EB,
                      "PP.EB.Beta" = binom.PP.EB.Beta,
                      "PP.EB.Sep" = binom.PP.EB.Sep,
                      "PP.Fix" = binom.PP.FIX,
                      "PP.Cor" = binom.PP.FB.MC.COR,
                      "Beta.EB" = binom.Beta.EB
                      )
v <- verbose

  do.call(f.to.call, args=alist(x=x,n=n, verbose=v, ...))

}