Displaying 1 result from an estimated 1 matches for "extractaicc".
Did you mean:
extractaic
2009 May 05
0
stepAICc function (based on MASS:::stepAIC.default)
...if (!length(scope))
stop("no terms in scope for adding to object")
ns <- length(scope)
ans <- matrix(nrow = ns + 1L, ncol = 2L, dimnames = list(c("<none>",
scope), c("df", "AICc")))
n=length(residuals(object))
edf=extractAICc(object, scale, k = k, ...)[1L]
ans[1, ] <- rbind(
edf,
extractAICc(object, scale, k = k, ...)[2L]+((2*edf*(edf+1))/(n-edf-1)) #calculate AICc
)
n0 <- length(object$residuals)
env <- environment(formula(object))
for (i in seq(ns)) {
tt <- scope[i]...