R/AIC.glmssn.R

In jayverhoef/SSN: Spatial Modeling on Stream Networks

AIC.glmssn <- function(object,...,k=2) {
    if(class(object)[[1]] != "glmssn") return("Not a glmssn object")
    if(!missing(k)) cat("This argument has no effect\n")
    if(object$args$family == "poisson") return(NA)
    if(object$args$family =="binomial") return(NA)
    cp <- covparms(object)[,3]
    nparmstheta <- length(cp)
    rankX <- object$sampinfo$rankX
    if(object$args$EstMeth == "REML")
        nparms <- nparmstheta
    if(object$args$EstMeth == "ML")
        nparms <- rankX + nparmstheta
    object$estimates$m2LL + k*nparms
}