R Under development (unstable) (2024-06-18 r86781 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > 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) + } Loading required package: MuMIn Fixed term is "(Intercept)" 0 : rlm(formula = y ~ 1, data = Cement) 1 : rlm(formula = y ~ X1 + 1, data = Cement) 2 : rlm(formula = y ~ X2 + 1, data = Cement) 3 : rlm(formula = y ~ X1 + X2 + 1, data = Cement) 4 : rlm(formula = y ~ X3 + 1, data = Cement) 5 : rlm(formula = y ~ X1 + X3 + 1, data = Cement) 6 : rlm(formula = y ~ X2 + X3 + 1, data = Cement) 7 : rlm(formula = y ~ X1 + X2 + X3 + 1, data = Cement) 8 : rlm(formula = y ~ X4 + 1, data = Cement) 9 : rlm(formula = y ~ X1 + X4 + 1, data = Cement) 10 : rlm(formula = y ~ X2 + X4 + 1, data = Cement) 11 : rlm(formula = y ~ X1 + X2 + X4 + 1, data = Cement) 12 : rlm(formula = y ~ X3 + X4 + 1, data = Cement) 13 : rlm(formula = y ~ X1 + X3 + X4 + 1, data = Cement) 14 : rlm(formula = y ~ X2 + X3 + X4 + 1, data = Cement) 15 : rlm(formula = y ~ X1 + X2 + X3 + X4 + 1, data = Cement) 22 : rlm(formula = y ~ X2 + X3 + X2:X3 + 1, data = Cement) 23 : rlm(formula = y ~ X1 + X2 + X3 + X2:X3 + 1, data = Cement) 30 : rlm(formula = y ~ X2 + X3 + X4 + X2:X3 + 1, data = Cement) 31 : rlm(formula = y ~ X1 + X2 + X3 + X4 + X2:X3 + 1, data = Cement) (Intercept) X1 X2 X3 X4 (Intercept) 657.3291410 -2.7536391499 -3.5010654693 -17.745284127 -7.683626226 X1 -2.7536391 0.0488568894 0.0034839126 0.150515200 0.024212557 X2 -3.5010655 0.0034839126 0.0403462579 0.148368689 0.097570988 X3 -17.7452841 0.1505151999 0.1483686893 0.227058309 0.188165359 X4 -7.6836262 0.0242125572 0.0975709883 0.188165359 0.082900719 X2:X3 0.0148248 0.0006357005 -0.0009424436 -0.003691552 0.003716187 X2:X3 (Intercept) 0.0148247985 X1 0.0006357005 X2 -0.0009424436 X3 -0.0036915519 X4 0.0037161868 X2:X3 0.0001039456 > > proc.time() user system elapsed 1.65 0.25 1.89