tests/rlm.R

In MuMIn: Multi-Model Inference

if (MuMIn:::testStart("MASS")) {

    data(Cement, package = "MuMIn")

    nseq <- function(x, len = length(x)) 
        seq(min(x, na.rm = TRUE), max(x, na.rm = TRUE), length = len)

    fm1 <- rlm(y ~X1 + X2 * X3 + X4, data = Cement)
    dd <- dredge(fm1, trace = TRUE)
    gm <- get.models(dd, subset = 1:10)
    ma <- model.avg(gm)
    stopifnot(all(predict(ma) == predict(ma, Cement)))
    predict(ma, lapply(Cement, nseq, len = 30), se.fit = TRUE)
    vcov(ma)
}