Question

当人们希望比较（测试）多个组时（例如，在做anova时），人们面临着多重比较的问题。如果我们希望绘制比较，同样适用。

我的问题是，您知道哪些工具（在R中）允许绘制反映多重比较的绘图？

目前，我只知道两个（虽然我相信还有更多）：

TukeyHSD（）结合plot（）
boxplot选择“notches”的方式

Answer 1

套餐multcomp有例如plot.cld() - 您可以尝试

library(multcomp)
example(plot.cld)

此外，http://rseek.org上的快速“多重比较图”搜索显示了更多的包和任务视图。

Answer 2

有一些方法可以在GLM中进行多重比较

http://www.r-bloggers.com/multiple-comparisons-for-glmms-using-glmer-glht/

有一篇关于R-Project统计分析手册（网站）同步推断的文章......

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)

绘制多重比较？

2 个答案: