R/calc_bic.R

In tece-lab/DAISIEutils: Utility Functions for the DAISIE Package

#' Compute Bayesian Information Criterion of a model
#'
#' @inheritParams default_params_doc
#'
#' @return Numeric value with BIC
#' @export
#'
#' @author Joshua W. Lambert, Pedro Santos Neves
calc_bic <- function(results, daisie_data) {
  k <- results$df
  if (is.null(daisie_data[[1]]$not_present)) {
    m <- daisie_data[[1]]$not_present_type1 +
      daisie_data[[1]]$not_present_type2 + length(daisie_data) - 1
  } else {
    m <- daisie_data[[1]]$not_present + length(daisie_data) - 1
  }
  bic <- k * (log(m) + log(2 * pi)) - 2 * results$loglik
  testit::assert(is.numeric(bic))
  testit::assert(identical(length(bic), 1L))
  bic
}