我想测试两个两个对比本身是否存在显着差异。
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
...我想看看这些差异本身是否彼此不同