data(Cement)
options(na.action = "na.fail") # change the default "na.omit" to prevent models
# from being fitted to different datasets in
# case of missing values.
fm1 <- lm(y ~ ., data = Cement)
ms1 <- dredge(fm1)
plot(ms1)
model.avg(ms1, subset = delta < 4)
confset.95p <- get.models(ms1, cumsum(weight) <= .95)
avgmod.95p <- model.avg(confset.95p)
summary(avgmod.95p)
confint(avgmod.95p)Run the code above in your browser using DataLab