R/dic.R

In DDPstar: Density Regression via Dirichlet Process Mixtures of Normal Structured Additive Regression Models

dic <-
function(y, X, res, termsum = NULL) {
  if(is.null(termsum)) {
    termsum <- inf_criteria(y, X, res)
  }
  
  logtermsum <- log(termsum)
  deviance <- -2* mean(apply(logtermsum, 1, sum)) 
  pd <- deviance + 2 * sum(log(apply(exp(logtermsum),2,mean)))
  DIC <- deviance + pd
  
  res <- list()
  res$pD <- pd
  res$DIC <- DIC
  res
}