Я хотел бы проверить, являются ли два двух контраста сами по себе существенно различающимися.
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
... Я хочу посмотреть, отличаются ли эти различия друг от друга