Я использовал пакет 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))
Выбор всех подмножеств
rttdASMS ‹-glmulti (rttdfit, confsetsize = 5, crit = aic, level = 1, family = binomial (logit))
взвешиваемый (rttdASMS)