Question

我正在使用lmer模型查看环境预测变量对景观变量的影响。为此，我正在使用挖泥函数创建一个包含所有可能的预测变量组合的模型候选集。

m3 <- lmer(div~scale(log(travel.time))+scale(spinsandplain)+scale(ThreeYearRain)+scale(claylake)+scale(ThreeYearRain)*scale(log(travel.time))+(1|circleID),na.action=na.fail,
          data=data, REML=FALSE)
s <-dredge(m3, extra = list("R^2"))
s

summary(get.models(s, 1)[[1]])

我现在想从每个顶级模型中提取每个变量的置信区间。除了模型平均以外，我似乎找不到任何代码。你有密码吗？这不可能吗？

预先感谢，莉安

Answer 1

get.models()返回与您的list相同类的模型对象global.model，因此请使用confint或通过lapply，sapply或for循环在每个项目上的任何相关功能。

例如：lapply(get.models(s, 1:10), confint)

挖泥函数顶层模型（lmer）的变量置信区间

1 个答案: