我一直在使用包glmulti来运行所有子集模型选择过程 . 我有一个大型数据集，我的全局模型包含7个候选预测变量 . 问题是我想考虑两个主要预测因子之间的单一交互 . 当我在全局模型函数中包含交互项时，当我使用'level = 1'仅指定考虑的主效应时，glmulti会忽略它 . 有没有人有办法让我可以考虑主要效果和单一互动的程序？
这是一个小样本：
GlobalModel rttdfit < - glm（ROCK~AREA DRATIO ELEV SECC BLDGD MAXD * CLIM I（CLIM ^ 2），data = rttd，family = binomial（link =“logit”））Allsubsets selection rttdASMS <-glmulti（rttdfit，confsetsize = 5，crit =“aic”，level = 1，family = binomial（“logit”））weightable（rttdASMS）