为了显示字母的显着差异,multcomp包中有一个绘图功能:
library(multcomp)
tuk <- glht(model, linfct = mcp(effect = "Tukey"))
plot(cld(tuk))
如何计算方框(即上下四分位数和胡须)?我在文档中找不到任何相关信息。
答案 0 :(得分:1)
这基本上只是一个普通的boxplot
。例如,使用warpbreaks
数据:
# load data
data(warpbreaks)
# model
amod <- aov(breaks ~ tension, data = warpbreaks)
tuk <- glht(amod, linfct = mcp(tension = "Tukey"))
# plot
plot(cld(tuk))
这与
几乎相同boxplot(breaks ~ tension, data = warpbreaks)
# OR
boxplot(split(warpbreaks$breaks, warpbreaks$tension))
boxplot
使用boxplot.stats
,stats.fivenum
使用quantile
来计算必要的统计信息。请注意,这与probs=c(0,.25, .5, .75, 1)
函数与na.rm=TRUE
和NA
略有不同。最简单的形式(没有fivenum
&#39; s,长度> 0),fivenum <- function (x)
{
x <- sort(x)
n <- length(x)
n4 <- floor((n + 3)/2)/2
d <- c(1, n4, (n + 1)/2, n + 1 - n4, n)
0.5 * (x[floor(d)] + x[ceiling(d)])
}
看起来像这样:
{{1}}
答案 1 :(得分:0)
看看代码:
plot.cld()
来电boxplot()
,因此您可以在?boxplot.stats
找到您要查找的信息。