R/NAME_bayesbox.R

In sebastiansauer/prada: Data and R functions for practical data analysis

#' Compute Bayes Box
#'
#' Similar to Allen Downey's concept of a Bayes Box, such a table is computed by this function.
#' It can be used for simple Bayes models of discrete (non-continuous) variables.
#'
#' @param hyps Vector of hypotheses (all types allowed), serves as ID
#' @param priors Vector of prior probabilites (numeric)
#' @param liks Vector of likelihoods (numeric)
#' @param round number of decimal places, NULL for not rounding (integer)
#'
#' @return data.frame with Bayes Box
#' @export
#'
#' @examples
#'
#' hyps <- seq(0,1, by = .1)
#' bayesbox(hyps = hyps, priors = 1, liks = dbinom(x = 5, size = 7, prob = hyps, round = 3))
bayesbox <- function(hyps, priors, liks, round = 2L) {

  out <-
    data.frame(
      hyps = as.factor(hyps),
      priors = priors,
      liks = liks)

  out$post_unstand <- priors*liks
  out$post_std = out$post_unstand / sum(out$post_unstand)

  return(out)
}