为什么二项式 glm 会给出负面预测?

机器算法验证 r 广义线性模型 二项分布
2022-03-24 16:08:40

我以非常简单的方式使用计数数据,但我无法理解二项式 glm 如何返回负面预测

示例代码,其中成功计数随着响应变量的增加而增加:

    suc=c(1:10)
    fail=c(10:1)
    predict(glm(cbind(suc,fail)~c(1:10),family=binomial))

这导致:

    -1.9974174 -1.5535469 -1.1096763 -0.6658058 -0.2219353  0.2219353  0.6658058 1.1096763  1.5535469  1.9974174 

我不明白这一点:二项式模型如何给出这些预测?它应该是整数正预测,不是吗?

1个回答

假设您使用的predict.glm()stats包中的。

手册中的引用,在解释类型参数的条目下:

因此,对于默认二项式模型,默认预测是对数赔率(logit 标度上的概率),并且“type = "response"' 给出了预测概率。

因此,请尝试以下操作:

predict(glm(cbind(suc,fail)~c(1:10),family=binomial), type="response")