我一直在使用包glmulti来运行所有子集模型选择过程。我有一个大型数据集,我的全局模型包含7个候选预测变量。问题是我想考虑两个主要预测因素之间的单一交互。当我在全局模型函数中包含交互项时,glmulti在使用' level = 1'时忽略它。仅指定考虑的主要效果。有没有人有办法让我可以考虑主要效果和单一互动的程序?
这是一个小样本:
GlobalModel
rttdfit< - glm(ROCK~AREA + DRATIO + ELEV + SECC + BLDGD + MAXD * CLIM + I(CLIM ^ 2),data = rttd,family = binomial(link =" logit") )
Allsubsets选择
rttdASMS< -glmulti(rttdfit,confsetsize = 5,crit =" aic",level = 1,family = binomial(" logit"))
weightable(rttdASMS)
答案 0 :(得分:0)
我一直在寻找有关此问题的信息,并在此处找到了一个可能对将来的读者有所帮助的链接: