R/resid_raw.R

In goodekat/redres: Residuals and Diagnostic Plots for Mixed Models

# rawres: This is a function that is used to compute the conditional or marginal raw
# residuals from the model. Default is the conditional residuals.

#' @importFrom stats predict

rawres <- function(model, cond = TRUE) {

  # Compute the conditional raw residuals
  if (cond == TRUE){
    res <- as.vector(model@resp$y - stats::predict(model))
  # Compute the marginal raw residuals
  } else {
    res <- as.vector(model@resp$y - (model@pp$X %*% matrix(model@beta, ncol = 1)))
  }

  # Return the computed residuals
  return(res)

}