Есть ли простой способ объединить две модели glm в R?

У меня есть две модели логистической регрессии в R, сделанные с glm(). Они оба используют одни и те же переменные, но были созданы с использованием разных подмножеств матрицы. Есть ли простой способ получить среднюю модель, которая дает средние значения коэффициентов, а затем использовать ее с функцией предиката ()?

[извините, если этот тип вопроса должен быть размещен на сайте программирования, дайте мне знать, и я опубликую его там]

Спасибо

Вы хотите взять среднее значение прогнозируемых вероятностей или среднее значение коэффициентов? Они будут давать разные результаты, потому что логистическая регрессия включает в себя нелинейное преобразование линейного предиктора.

Функция для выполнения любого из них будет выглядеть примерно так. Установите, avgчтобы "prob"получить первое, или что-то еще для последнего.

pred_comb <- function(mod1, mod2, dat, avg="prob", ...)
{
    xb1 <- predict(mod1, dat, type="link", ...)
    xb2 <- predict(mod2, dat, type="link", ...)
    if(avg == "prob")
        (plogis(xb1) + plogis(xb2))/2
    else plogis((xb1 + xb2)/2)
}