misc/lmbeta.R

In Rkabacoff/qacReg: Multiple Regression and Logistic Regression

# from lm.beta package
lm.beta <- function (object)
{
  if (!"lm" %in% attr(object, "class"))
    stop("object has to be of class lm")
  object$standardized.coefficients <- coef(object) * apply(as.matrix(model.matrix(object)),
                                                           2, function(x) sqrt(sum((x - mean(x, na.rm = T) * attr(attr(object$model,
                                                                                                                       "terms"), "intercept"))^2, na.rm = T)))/apply(as.matrix(model.frame(object)[,
                                                                                                                                                                                                   1]), 2, function(x) sqrt(sum((x - mean(x, na.rm = T) *
                                                                                                                                                                                                                                   attr(attr(object$model, "terms"), "intercept"))^2,
                                                                                                                                                                                                                                na.rm = T)))
  attr(object, "class") <- c("lm.beta", "lm")
  return(object)
}