基于配对测试的紧凑型字母显示

Question

我想用我在线性混合效应模型(lmer)上做的临时测试来创建一个紧凑的字母显示器。

下面是我在做成对的t.test时想要做的一个例子

df <- read.table("https://pastebin.com/raw/Dzfh7b2f", header=T,sep="")

library(rcompanion)
library(multcompView)

PT <- pairwise.t.test(df$fit,df$treatment, method=bonferroni)
PT = PT$p.value
PT1 = fullPTable(PT)
multcompLetters(PT1,
                compare="<",
                threshold=0.05,
                Letters=letters,
                reversed = FALSE)

--这是我们伟大的，因为从pairwise.t.test中提取p值和创建我想要的表很容易。

现在假设我运行一个线性模型，做一个成对的比较，并且像上面一样创建一个表，它从提取的pvalue中为我创建一个紧凑的字母显示。

library(multcomp)
mult<- summary(glht(model, linfct = mcp(treatment = "Tukey")), test = adjusted("holm"))
mult

我可以看到p值，但在过去的2-3个小时里，我一直试图找出如何只提取这些值(就像我在上面使用pairwise.t.test时所做的那样)，然后创建一个紧凑的字母显示表。

任何帮助都是非常感谢的。万事如意

Answer 1

多亏了@roland的建议，答案很简单：

mult<- summary(glht(model, linfct = mcp(treatment = "Tukey")), test = adjusted("holm"))

letter_display <- cld(mult)

letter_display