自写R代码整合了单因素ANOVA(整体的多组间比较)、Tukey HSD多重比较(事后两两比较)以及ggplot2作图,并且用不同的显著性标记方法标记,,便于自动完成一系列的需要。
时间: 2024-05-29 10:10:44 浏览: 19
以下是自写的R代码:
#导入必要的包
library(ggplot2)
library(multcomp)
#设置数据
data <- data.frame(group = rep(c("A", "B", "C"), each = 10),
value = c(rnorm(10, 20, 3), rnorm(10, 25, 3), rnorm(10, 30, 3)))
#进行单因素ANOVA
model <- aov(value ~ group, data = data)
summary(model)
#进行Tukey HSD多重比较
tukey <- TukeyHSD(model)
tukey
#利用ggplot2进行作图
plot <- ggplot(data, aes(x = group, y = value)) +
geom_boxplot(fill = "gray") +
labs(title = "Comparison of Groups A, B, and C",
x = "Group", y = "Value") +
theme_classic()
#添加显著性标记
plot +
geom_signif(comparisons = list(c("A", "B"), c("A", "C"), c("B", "C")),
textsize = 6, tip_length = 0.02) +
geom_text(data = tukey$`group`, aes(x = group1, y = group2, label = as.character(p.adj)),
color = "red", size = 4, position = position_nudge(x = -0.2, y = -0.2))