QAICc for select_Millar

[tokami]

check + include

QAICc <- function(mod){
  ll <- as.numeric(mod$value) 
  K <- length(mod$estimates[,1])   
  n <- nnzero(CatchPerNet_mat) 
  c_hat = as.numeric(mod$out[4]/mod$out[6])
  qaicc = -2*ll/c_hat + 2*K + 2*K*(K + 1)/(n - K - 1)
  qaicc
}