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

机器算法验证 r 广义线性模型
2022-03-08 07:04:47

我在 R 中有两个逻辑回归模型,使用glm(). 它们都使用相同的变量,但使用矩阵的不同子集制作。有没有一种简单的方法来获得一个给出系数平均值的平均模型,然后将它与 predict() 函数一起使用?

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

谢谢

1个回答

你想取预测概率的平均值,还是系数的平均值?它们会给出不同的结果,因为逻辑回归涉及线性预测变量的非线性变换。

做任何一个的功能都是这样的。设置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)
}