是否可以让glmulti仅考虑所有子集模型选择中的选择入口?

时间:2016-06-16 15:35:56

标签: r glm consensus

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

1 个答案:

答案 0 :(得分:0)

我一直在寻找有关此问题的信息,并在此处找到了一个可能对将来的读者有所帮助的链接:

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