有没有一种简单的方法可以在R中组合两个glm模型？

我在用制成的R中有两个逻辑回归模型glm()。它们都使用相同的变量，但使用不同的矩阵子集。是否有一种简单的方法来获取给出系数均值的平均模型，然后将其与predict（）函数一起使用？

[抱歉，如果此类问题应发布在编程网站上，请告诉我，我会在此处发布]

谢谢

您是否要获取预测概率的平均值或系数的平均值？它们将给出不同的结果，因为逻辑回归涉及线性预测变量的非线性变换。

一个可以做的功能就是这样。设置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)
}