R/bic.R

In tharris0924/btld: Triangular and Bimodal Triangular Linked Distributions

#' Bayesian Information Criterion (BIC)
#'
#' Calculate the BIC for a model.
#'
#' @param loglik log-likelihood
#' @param k number of paramters
#' @param n sample size
#' @return BIC
#' @import mixtools
#' @export
#' @examples
#' gmm<-mixtools::normalmixEM(runif(1000))
#' bic(gmm$loglik,3,1000)



bic<-function(loglik,k,n){
  return((-2*loglik)+(k*log(n)))
}