泊松回归斜率估计的解释

机器算法验证 r 广义线性模型 泊松回归
2022-03-28 12:21:04

我看过其他关于解释泊松回归斜率估计的帖子,基于此,这是我对泊松回归的解释。如果我能写下下面的语句,谁能告诉我:

mdl<-glm(y ~ year, family="quasipoisson")
summary(mdl) 

Deviance Residuals: 
Min       1Q   Median       3Q      Max  
-2.8258  -1.4108  -0.5760   0.8562   3.2575  

Coefficients:
         Estimate Std. Error t value  Pr(>|t|)    
(Intercept) 68.918820  14.838432   4.645 0.0000216 ***
year        -0.034351   0.007531  -4.561 0.0000289 ***

slope<-exp(mdl$coefficients[2])
# 0.9662321

slope - 1
# -0.03376786

我可以做如下陈述:增加一个单位year导致y减少一个0.033单位?

谢谢

1个回答

没有。有两个问题,一个是算术到英文的翻译,一个是哲学问题。

短语“按0.033单位减少”应解释为“减去0.033单位y”,这是不正确的。

更好的是

增加一个单位对应year于乘以y0.966

或者

增加一个单位year对应于3.3%减少y

百分号非常重要,它携带的信息是减少是乘法的,而不是加法的。

另一个是你对原因这个词的使用。这是一个沉重的词,你不应该在没有认真考虑的情况下使用它。当然,回归不能单独证明因果关系,它必须与随机分配的实验或一些科学推理相结合才能相信因果关系存在(在这种情况下,回归只是估计其影响大小)。