如何从R中的大量统计结果中提取几个级别？

Question

我不确定问题标题中是否使用了正确的词，所以请让我解释一下。 可以说，我有一个包含三个因素的数据集：物种，生命和暴露以及一个响应变量：mean_v。 总共有108个级别。 当我在Tukey的事后检验中应用三向方差分析时，它会产生大量的统计比较。

life <- rep(c("5d", "15d", "45d"), 2, each = 18)
species <- rep(c("SP1", "SP2"), each = 54)
exposure <- rep(c("c1", "c2", "c3"), 6, each = 6)
mean_v <- runif(108, 4, 80)
data1 <- data.frame(species, life, exposure, mean_v)
model1 <- aov(mean_v~species*life*exposure, data1)
tukey <- TukeyHSD(model1)
tukey

但是，我只对几个感兴趣，例如SP2：5d：c3-SP1：5d：c3之间的比较。 有没有办法从整个分析中提取特定的比较？ 提前致谢

Answer 1

您必须选择tukey内的species:life:exposure list，然后才选择所需的比较，例如：

tukey$`species:life:exposure`[rownames(tukey$`species:life:exposure`)=="SP2:5d:c3-SP1:5d:c3",]
      diff        lwr        upr      p adj 
 -2.123986 -47.461693  43.213722   1.000000