0

我正在使用 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]])  

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

在此先感谢,

4

1 回答 1

2

get.models()返回list与您的 相同类的模型对象,因此通过或循环对每个项目global.model使用 eg或任何相关函数。confintlapplysapplyfor

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

可重现的例子:

library(glmmTMB)
library(MuMIn)

# from example(glmmTMB)
m2 <- glmmTMB(count ~ spp + mined + (1|site), family=nbinom2, data=Salamanders))
models <- get.models(dredge(m2), TRUE)

# list of CI for each model's parameters:
lapply(models, confint)
于 2019-12-04T12:11:48.090 回答