Я создал обобщенную линейную модель с отрицательным биномиальным отображением, используя функцию glm.nb:
glm_mod <- glm.nb(obs ~ offset(log(time_interval)) + pred1*pred2, data = df)
Я хотел бы смоделировать скорость (наблюдения, проводимые в единичный интервал времени) на основе pred1, pred2 и их взаимодействия. Эта часть работает нормально.
Далее я бы хотел представить эффекты взаимодействия между двумя предикторами pred1:pred2 с использованием пакета effects. Это - то, когда я сталкиваюсь с проблемой, потому что кажется, что в настоящее время этот пакет не может обрабатывать переменные смещения, определенные в формуле, и функция glm.nb не принимает и смещает переменную как аргумент, как стандартная функция glm.
Есть ли способ обойти эту несовместимость между effects и glm.nb?
Есть ли другие соответствующие пакеты, которые я могу использовать вместо этого?