Можно ли сделать, чтобы glmulti учитывал только выборочные влияния во всех моделях подмножеств?

Question

Я использовал пакет glmulti для запуска процедуры выбора модели всех подмножеств. У меня большой набор данных, и моя глобальная модель содержит 7 кандидатов-предикторов. Проблема в том, что я хотел бы рассмотреть только одно взаимодействие между двумя моими основными предикторами. Когда я включаю термин взаимодействия в функцию глобальной модели, glmulti игнорирует его, когда я использую 'level = 1', обозначая только основные рассматриваемые эффекты. У кого-нибудь есть способ, которым я могу заставить процедуру рассмотреть основные эффекты и единственное взаимодействие?

Вот небольшой образец:

GlobalModel

rttdfit<- glm (ROCK ~ ОБЛАСТЬ + 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)

Answer 1

Я искал по этому вопросу и нашел здесь ссылку, которая может помочь будущим читателям:

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