1

我想测试两个两个对比本身是否存在显着差异。

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 

...我想看看这些差异本身是否彼此不同

4

0 回答 0