cohens_d_r <- function(x, object) UseMethod("cohens_d_r", x)

cohens_d_r.lm <- function(x, object) {
  x / sd(resid(object))
}

cohens_d_r.mlm <- function(x, object) {
  x / apply(resid(object$lm), 2, sd)
}