https://github.com/cran/MuMIn
Tip revision: a0f20e50f08f1cb37bc636092dfc8f46eb3f0187 authored by Kamil BartoĊ on 14 June 2009, 00:00:00 UTC
version 0.12.2
version 0.12.2
Tip revision: a0f20e5
getAICc.R
`getAICc` <-
function(model, k = 2) {
if (any(inherits(model, c("lmer", "glmer")))) {
mLogLik <- logLik(model, model@status["REML"])
N <- NROW(model@frame)
} else {
mLogLik <- logLik(model)
N <- length(resid(model))
}
mK <- attr(mLogLik, "df")
mAIC <- -2 * c(mLogLik) + k * mK
ret <- mAIC + 2 * mK * (mK + 1)/(N - mK - 1)
attr(ret, "df") <- mK
attr(ret, "AIC") <- mAIC
return (ret)
}