search

tokami / TropFishR

Tropical Fisheries Analysis with R
24 stars 19 forks source link

QAICc for select_Millar #39

Open tokami opened 3 years ago

tokami commented 3 years ago

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
}