av.frame.lm: lm method for av.frame

In gmonette/spida: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

av.frame.lm(model, variable, ...)

Arguments

model
variable
...

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (model, variable, ...) 
{
    mod.mat <- model.matrix(model)
    var.names <- colnames(mod.mat)
    omit <- grep(variable, var.names)
    if (0 == length(omit)) 
        stop(paste(variable, "is not matched among columns of the model matrix."))
    cat("x.var =", var.names[omit[1]], "\n", "omitted vars =", 
        var.names[omit[-1]], "\n")
    response <- response(model)
    x.var <- mod.mat[, omit[1]]
    Xpred <- mod.mat[, -omit]
    preds <- predict(update(model, na.action = na.exclude))
    responseName <- responseName(model)
    if (is.null(weights(model))) 
        wt <- rep(1, length(response))
    else wt <- weights(model)
    res <- lsfit(mod.mat[, -omit], cbind(mod.mat[, omit[1]], 
        response), wt = wt, intercept = FALSE)$residuals
    ret <- matrix(NA, nrow = length(preds), ncol = 2)
    ret[!is.na(preds), ] <- res
    data.frame(x.res = ret[, 1], y.res = ret[, 2])
  }

gmonette/spida documentation built on May 17, 2019, 7:25 a.m.