要在R中执行方差分析,我通常遵循两个步骤:
1(我用函数aov计算方差分析摘要2(我重新组织数据聚合主题和条件以可视化情节
我想知道是否总是有必要对数据进行重组以查看结果,或者是否存在快速绘制结果的 f(x(。
感谢您的建议
G.
我想你的意思是用数字来说明你的测试结果?方差分析通常用箱线图来说明。
set.seed(1234)
data <- data.frame(group = c(rep("group_1",25),rep("group_2",25)), scores = c(runif(25,0,1),runif(25,1.5,2.5)))
mod1<-aov(scores~group,data=data)
summary(mod1)
您可以使用实现的函数图或箱线图制作箱线图
boxplot(scores~group,data=data)
plot(scores~group,data=data)
或者用 ggplot
require(ggplot2)
require(ggsignif)
ggplot(data, aes(x = group, y = scores)) +
geom_boxplot(fill = "grey80", colour = "blue") +
scale_x_discrete() + xlab("Group") +
ylab("Scores") +
geom_signif(comparisons = list(c("group_1", "group_2")),
map_signif_level=TRUE)
希望这有帮助