R/devianceGLM.r

#' @export
#' @title devianceGLM
#' @description couldn't accurately describe
#' @family statistical
#' @param \code{x} the glm model object
#' @param \code{data} 
#' @param \code{test} defaults to FALSE
#' @param \code{digits} defaults to 2
#' @return \code{result}
#' @author  unknown, \email{<unknown>@@dfo-mpo.gc.ca}
#' @export
devianceGLM <- function (x, data=NULL, test = "F", digits = 2) {
        if(is.null(data)) data = x$data
    #deviance explained by glm model 
    nexpl <- x$deviance
    tot <- x$null.deviance
    expl <- tot - nexpl
    ratio <- round((expl/tot * 100), digits = digits)
    expl <- round(expl, digits = digits)
    tot <- round(tot, digits = digits)
    cat("Deviance explained: ", expl, "/", tot, "(", ratio, "percent)\n\n")
    formula0 <- paste(names(x$model)[1], "~ 1")
    model0 <- glm(as.formula(formula0), family = x$family, data = na.omit(data))
    result <- anova(model0, x, test = test)
    return(result)
    }
AtlanticR/bio.utilities documentation built on June 21, 2020, 7:43 p.m.