Контраст в контрастах в emmeans - PullRequest
1 голос
/ 06 января 2020

Я хотел бы проверить, являются ли два двух контраста сами по себе существенно различающимися.

library(tidyverse)
library(emmeans)

data(mtcars)

lm1 <- lm(
  mpg ~ factor(vs) * factor(am),
  data = mtcars
  )

lm1 %>% 
  emmeans(
    pairwise ~ vs | am
  )

возвращает разницу в mpg для каждого значения vs

$contrasts
am = 0:
 contrast estimate   SE df t.ratio p.value
 0 - 1       -5.69 1.65 28 -3.448  0.0018 

am = 1:
 contrast estimate   SE df t.ratio p.value
 0 - 1       -8.62 1.93 28 -4.464  0.0001 

... Я хочу посмотреть, отличаются ли эти различия друг от друга

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...