當你想比較(測試)多個組(例如,在做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