1

我一直在使用glmulti包来运行所有子集模型选择过程。我有一个大型数据集,我的全局模型包含 7 个候选预测变量。问题是我只想考虑两个主要预测变量之间的单一交互。当我在全局模型函数中包含交互项时,glmulti 在我使用 'level = 1' 指定仅考虑的主要影响时会忽略它。有没有人可以让我获得考虑主要影响和单一交互的程序?

这是一个小样本:

全球模型

rttdfit<- glm(ROCK~AREA+DRATIO+ELEV+SECC+BLDGD+MAXD*CLIM+I(CLIM^2),data=rttd,family=binomial(link="logit"))

所有子集选择

rttdASMS<-glmulti(rttdfit, confsetsize = 5, crit = "aic", level = 1, family = binomial("logit"))

可加权(rttdASMS)

4

1 回答 1

0

我正在搜索这个问题,并在这里找到了一个可能对未来读者有所帮助的链接:

https://vcalcagnoresearch.wordpress.com/package-glmulti/

于 2020-02-11T15:24:45.130 回答