当你想比较(测试)多个组(例如,在做anova时),就会遇到多重比较的问题。如果我们想绘制比较,这同样适用。绘制多重比较?
因此,我的问题是,你知道哪些工具(R)可以反映多重比较的绘图?
目前,我知道只有两个(虽然我相信有更多):
- TukeyHSD()与图()
- 方式箱线图选择了 “缺口”
当你想比较(测试)多个组(例如,在做anova时),就会遇到多重比较的问题。如果我们想绘制比较,这同样适用。绘制多重比较?
因此,我的问题是,你知道哪些工具(R)可以反映多重比较的绘图?
目前,我知道只有两个(虽然我相信有更多):
有一些方法围绕在GLMS多重比较
http://www.r-bloggers.com/multiple-comparisons-for-glmms-using-glmer-glht/
有一个关于从统计分析(网站)的R-项目手册同时推断文章...
http://cran.r-project.org/web/packages/HSAUR2/vignettes/Ch_simultaneous_inference.pdf
来自gplot包的plotmeans()。这包括置信区间。
然后是“psych”包的error.bars.by()函数。从数据框中分组绘制手段和SD。
一些使用密度图进行可视化。
# Compare MPG distributions for cars with
# 4,6, or 8 cylinders
library(sm)
attach(mtcars)
# create value labels
cyl.f <- factor(cyl, levels= c(4,6,8),
labels = c("4 cylinder", "6 cylinder", "8 cylinder"))
# plot densities
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")
title(main="MPG Distribution by Car Cylinders")
# add legend via mouse click
colfill<-c(2:(2+length(levels(cyl.f))))
legend(locator(1), levels(cyl.f), fill=colfill)
包multcomp有eg plot.cld()
- 你可以尝试
library(multcomp)
example(plot.cld)
此外,在http://rseek.org快速“多重比较暗算”搜索发现多了一些包和任务视图。
谢谢德克。 我会给他们看看。 – 2010-02-18 21:57:23