如何在 glm Tweedie 回归模型中使用 emmeans？

Question

我和家人Tweedie调整了一个glm。现在我需要获得响应量表的调整均值和标准误差。到目前为止：

mtcars$am <- as.factor(mtcars$am)
m1 <- glm(hp ~ am, data = mtcars,
          family = statmod::tweedie(var.power = 2.2, 
                                    link.power = 0 #link log
                                    ))
summary(emmeans(m1, "am"), type = "response")

 am response      SE  df asymp.LCL asymp.UCL
 0    5.0768 0.11516 Inf    4.8511    5.3025
 1    4.8430 0.13600 Inf    4.5764    5.1095

Unknown transformation "mu^0": no transformation done 
Confidence level used: 0.95

如何告知我使用了 link log()？如何处理："Unknown transformation "mu^0": no transformation done"?

Answer 1

这样做：

emm <- update(emmeans(m1, "am"), tran = "log")
summary(emm, type = "response")

有关详细信息，请参阅转换中的插图。

如何在 glm Tweedie 回归模型中使用 emmeans？

How to use emmeans in a glm Tweedie regression model?

r

glm

emmeans