Извлечь вероятности из парного контраста в emmeans

Question

Я пытаюсь выполнить парное сравнение после glm с биномиальным dv, а emmeans сообщает соотношение шансов, тогда как мне требуется разница в вероятностях.

library(magrittr)
library(emmeans)

glm(
  am ~  wt * factor(vs), 
  family = binomial(),
  data = mtcars
) %>% 
  emmeans(
~ vs | wt, at = list(wt = seq(2.6, 3.6, length.out = 10))
  ) %>% 
  pairs(type = "response")

возвращает

wt = 2.6:
 contrast odds.ratio         SE  df z.ratio p.value
 0 / 1     668.84887 3388.29452 Inf   1.284  0.1991

wt = 2.71111111111111:
 contrast odds.ratio         SE  df z.ratio p.value
 0 / 1     453.70452 2028.70137 Inf   1.368  0.1713

Answer 1

Добавить %>% regrid() до pairs().Вызов regrid() преобразует все результаты в новый набор EMM, как будто никогда не было преобразования.

См. Виньетку о трансформациях.