Variation in Term Predictions

code.R

predEffects <- function(obj, digits=3){
    ## obj is a linear model object
    preds <- predict(obj, type="terms")
    rg <- apply(preds, 2, range)
    diff.rg <- apply(rg, 2, diff)
    out <- data.frame(
        min = rg[1,], 
        max = rg[2,], 
        diff = diff.rg, 
        iqr = apply(preds, 2, function(x)diff(quantile(x, c(.25,.75)))),
        sd = apply(preds, 2, sd)
    )
    print.data.frame(out, digits=digits)
}