我一直在使用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)