https://github.com/cran/MuMIn
Raw File
Tip revision: a0f20e50f08f1cb37bc636092dfc8f46eb3f0187 authored by Kamil BartoĊ„ on 14 June 2009, 00:00:00 UTC
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)
}

back to top