#' @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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.