R/BIC.drcOrdinal.R

In DoseResponse/bmd: Benchmark dose estimation for dose-response data

BIC.drcOrdinal <- function(object, epsilon = 10^(-16)){
  n.parameters <- sum(sapply(object$drmList, function(mod) length(mod$coefficients)))
  n.obs <- sum(object$data[[object$weights]])
  n.parameters * log(n.obs) - 2 * logLik(object, epsilon) 
}