#' Formats a model as an equation (string), for
#' description purposes. It shows the linear expectation of the mean, and assumes you have not manipulated
#' default assumptions (e.g., via offsets and link functions).
#'
#' @param object The model.
#' @importFrom stats coef formula
#' @export
Equation <- function(object)
{
coefs <- coef(object)
parameter.names <- names(coefs)
dependent.name <- all.vars(formula(object$call))[1]
parameter.names[1] <- "" #Intercept.
signs <- sign(coefs)
operator <- rep(" + ", length(signs))
if (sum(signs == -1) > 0)
operator[signs == -1] <- " - "
operator[1] <- ifelse(signs[1] == 1, "", " -")
coefs <- FormatAsReal(abs(coefs))
equation <- paste0(dependent.name, " = ",
paste0(operator, coefs, parameter.names, collapse = ""))
if (object$type == "Poisson" || object$type == "Quasi-Poisson" )
equation <- paste0("exp(", equation, ")")
strwrap(equation)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.