Question

我想从glht对象获取p值和相关因子（最后的示例代码）

summary(mcptu)$test$pvalues适用于p值，但rownames(summary(mcptu)$linfct)很棘手。

使用rownames(summary(mcptu)$linfct)我得到“2 - 1”“3 - 1”“4 - 1”“5 - 1”“6 - 1”“3 - 2”我可以拆分等等。但是可能有一种优雅的方式来获得被比较的因素？

library(mratios); library(multcomp)

data(Penicillin)
Penicillin$strain <- as.factor(Penicillin$strain)

linmod <- lm(diameter ~ strain, Penicillin)
mcptu  <- glht(linmod, mcp(strain="Tukey"))
summary(mcptu)

THX 克里斯托夫

Answer 1

对只是（2,1），（3,1），...（n，1），然后是（3,2），...（n，2），......（n， N）。使用seq和rep变体创建它们。

get_level_pairs <- function(n)
{
  data.frame(
    x = unlist(lapply(seq_len(n)[-1], seq.int, to = n)),
    y = rep.int(seq_len(n - 1), rev(seq_len(n - 1)))
  )
}

get_level_pairs(nlevels(Penicillin$strain))

或拆分这些rownames。

library(stringr)
str_split_fixed(rownames(mcptu$linfct), " - ", 2)

从glht对象获取详细信息

1 个答案: