R/classes_ineq.R

In multinomineq: Bayesian Inference for Multinomial Models with Inequality Constraints

# ' Count of Samples that Satisfy Inequality Constraints
# '
# ' The class \code{ineq_count} defines a class for the Bayes factor for inequality-constrained Bayes factors.
# '
# ' @name ineq_count-class
# ' @rdname ineq_count-class
# ' @exportClass ineq_count


as_ineq_count <- function(count) {
  class(count) <- c("ineq_count", "matrix")
  attr(count, "proportion") <- prod(count[, "count"] / count[, "M"])
  s_prop <- sampling_proportion(count = count[, "count"], M = count[, "M"], log = FALSE)
  attr(count, "se") <- sd(s_prop)
  attr(count, "ci90") <- quantile(s_prop, c(.05, .95))
  count
}

#' @export
print.ineq_count <- function(x, ...) {
  cat("Number of samples satisfying the inequality constraints:\n")
  print.table(x)
  p <- attr(x, "proportion")
  se <- attr(x, "se")
  cat("\nTo extract the proporiton of samples, use:\n   attr(count, \"proportion\") = ",
    p, " (SE = ", se, ").\n",
    sep = ""
  )
}


## #' @export
## show.ineq_count <- function(object){
##   print(object)
## }