我在 R 中创建了一个简单的线性回归模型,其中我对响应进行了对数转换。自变量是有 80 个级别的分类变量。然后我使用 emmeans 进行成对测试。如何将 emmeans 反向转换为与响应具有相同的比例(而不是对数比例)?我在 R 中使用 emmeans 包。
最佳答案
我认为这只是
summary(emmeans(...), type = "response")
?有一个vignette on transformations and link functions
例如
library(emmeans)
m1 <- lm(log(mpg) ~ factor(cyl), mtcars)
summary(emmeans(m1, ~cyl), type = "response")
cyl response SE df lower.CL upper.CL
4 26.3 1.268 29 23.9 29.0
6 19.7 1.190 29 17.4 22.3
8 14.9 0.636 29 13.6 16.2
Confidence level used: 0.95
Intervals are back-transformed from the log scale
请注意,如果您使用pairs(emmeans(m1, ~cyl), type = "response")
,它将为您提供成对比率,但请注意测试是在对数刻度上执行的(正如他们应该的那样)
关于r - R 中的反向变换 lsmeans,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/76323429/