residuals.grm <-
function (object, resp.patterns = NULL, order = TRUE, ...) {
if (!inherits(object, "grm"))
stop("Use only with 'grm' objects.\n")
if (any(is.na(object$X)))
warning("residuals are not meaningful for patterns with missing data.\n")
fits <- fitted(object, resp.patterns = resp.patterns)
X <- fits[, -ncol(fits), drop = FALSE]
Exp <- fits[, "Exp"]
betas <- object$coefficients
Obs <- observedFreqs(object, X)
out <- cbind(X, Obs = Obs, Exp = Exp, Resid = (Obs - Exp) / sqrt(Exp))
if (order)
out <- out[order(out[, "Resid"]), ]
out
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.