多核包中的plot.cld如何计算框?

时间:2014-08-18 11:40:23

标签: r boxplot mixed-models

为了显示字母的显着差异,multcomp包中有一个绘图功能:

library(multcomp)
tuk <- glht(model, linfct = mcp(effect = "Tukey"))
plot(cld(tuk))

如何计算方框(即上下四分位数和胡须)?我在文档中找不到任何相关信息。

2 个答案:

答案 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.statsstats.fivenum使用quantile来计算必要的统计信息。请注意,这与probs=c(0,.25, .5, .75, 1)函数与na.rm=TRUENA略有不同。最简单的形式(没有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找到您要查找的信息。