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

Question

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

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

您必须 select species:life:exposure 在 tukey 中列出，然后 select 只需要比较，像这样：

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

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

How to extract a few levels from a large statistical result in R?

r

extract