解释堵塞逻辑回归的估计

机器算法验证 物流 回归系数
2022-01-19 19:37:21

有人可以建议我如何使用 cloglog 链接解释逻辑回归的估计值吗?

我已经安装了以下模型lme4

glm(cbind(dead, live) ~ time + factor(temp) * biomass,
    data=mussel, family=binomial(link=cloglog))

例如,时间的估计值为 0.015。说每单位时间的死亡率乘以 exp(0.015) = 1.015113(每单位时间增加约 1.5%)是否正确。
换句话说,在 logit 逻辑回归中是否以 log 几率表示在 cloglog 中获得的估计值?

1个回答

使用互补-对数-对数链接功能,它不是逻辑回归——术语“逻辑”意味着一个逻辑链接。当然,它仍然是二项式回归。

时间估计为 0.015。说每单位时间的死亡率乘以 exp(0.015) = 1.015113 是否正确(每单位时间增加约 1.5%)

不,因为它没有根据对数赔率建模。这就是您使用 logit 链接所拥有的;如果您想要一个在对数赔率方面有效的模型,请使用 logit-link。

补充日志链接功能说

η(x)=log(log(1πx))=xβ

其中πx=P(Y=1|X=x)

所以不是优势比;确实exp(η)exp(η)=log(1πx)

因此因此,如果您需要某个特定的优势比,您可以计算一个,但就对数优势的贡献而言,参数没有直接的简单解释。exp(exp(η))=(1πx)1exp(exp(η))=πxx

相反(不出所料)一个参数显示(对于的单位变化)对互补对数对数的贡献。x


正如本在评论中的问题中轻轻暗示的那样:

说单位时间内死亡的概率(即危险)增加了 1.5% 是真的吗?

互补对数模型中的参数在风险比方面确实有一个简洁的解释。我们有:

eη(x)=log(1πx)=log(Sx),其中是生存函数。S

(因此,在示例中,每单位时间的对数生存率将下降约 1.5%。)

现在危险,,所以确实在例子中问题中给出的,每单位时间的死亡概率*增加了约 1.5%h(x)=ddxlog(Sx)=ddxeη(x)

*(或更一般地,对于具有 cloglog 链接的二项式模型,P(Y=1)