Я работаю с набором данных, достаточно большим, чтобы использовать более продвинутые алгоритмы, чем stats glm, для вычисления моделей биномиальной регрессии. glm вычисление занимает больше недели.
Большая биномиальная регрессия с фиксированными эффектами в R предполагает rxGlm и speedglm. Я попробовал оба и мне удалось добиться значительного улучшения скорости (до менее чем 2 часов). Однако я столкнулся со следующей проблемой:
В обоих случаях я не могу использовать функцию effects plot(effect()) для создания графиков взаимодействия, подобных следующим, потому что effects не может работать с такими объектами и/или отсутствуют необходимые данные.
Я уже разместил вопрос с вознаграждением на rxGlm (Как отображать эффекты взаимодействия из очень больших наборов данных (особенно из выходных данных rxGlm)), но безрезультатно. Так что я думаю, что нет способа сделать это с помощью rxGlm.
Теперь мой вопрос: могу ли я отобразить такие термины эффекта взаимодействия, используя speedglm, и если да, то как?
