R:箱线图
两组对比,t.test计算差异
rm (list = ls ()) #清除所有变量
library(ggplot2)
library(ggpubr)
setwd("C:\\Users\\Administrator\\Desktop\\新建文件夹") #设置工作目录
index <- read.table('group.txt', header = TRUE, row.names = 1)
my_comparisons = list( c("B73", "Mo17") )
# 绘制箱线图并美化
p <- ggboxplot(index, x = "Group", y = "Shannon", color = "Group", add = "jitter") +
stat_compare_means(comparisons = my_comparisons, method = "t.test", label = "p.signif") +
theme_minimal(base_size = 14) + # 使用简洁主题,调整基础字体大小
labs(title = "Comparison of Shannon Diversity",
x = "Group",
y = "Shannon Index") + # 添加标题和坐标轴标签
theme(legend.position = "top", # 将图例放置在顶部
panel.grid.major = element_blank(), # 去掉主要网格线
panel.grid.minor = element_blank(), # 去掉次要网格线
panel.border = element_blank(), # 去掉面板边框
axis.line = element_line(colour = "black")) + # 显示坐标轴的线
scale_color_manual(values = c("B73" = "#8FC9E2", "Mo17" = "#ECC97F")) # 自定义颜色
# 显示图形
print(p)
# 保存图形为PNG文件
ggsave("boxplot_shannon_diversity.png", plot = p, width = 10, height = 8, dpi = 600, bg = "white")
多组对比,AVNOA计算差异
rm (list = ls ()) #清除所有变量
library(ggplot2)
library(ggpubr)
setwd("C:\\Users\\Administrator\\Desktop\\新建文件夹") #设置工作目录
index <- read.table('group.txt', header = TRUE, row.names = 1)
# 进行ANOVA分析并进行Tukey多重比较
anova_result <- aov(Simpson ~ Group, data = index)
tukey_result <- TukeyHSD(anova_result)
print(summary(anova_result))
print(tukey_result)
# 手动设置基于Tukey多重比较结果的标签
# 注意:需要根据实际的Tukey测试结果来手动设置这些标签
group_labels <- data.frame(
Group = c("Mo17_DAS28", "Mo17_DAS42", "Mo17_DAS56", "Mo17_DAS70"),
label = c("a", "a", "a", "b") # 所有组别都标记为"a"
)
# 绘制箱线图并美化
p <- ggboxplot(index, x = "Group", y = "Simpson", color = "Group", add = "jitter") +
theme_minimal(base_size = 14) +
labs(title = "Comparison of Simpson Diversity", x = "Group", y = "Simpson Index") +
theme(legend.position = "top",
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.border = element_blank(),
axis.line = element_line(colour = "black")) +
scale_color_manual(values = c("Mo17_DAS28" = "#9BC985", "Mo17_DAS42" = "#F7D58B", "Mo17_DAS56" = "#B595BF", "Mo17_DAS70" = "#797BB7")) +
geom_text(data = group_labels, aes(x = Group, y = max(index$Simpson) * 1.05, label = label), vjust = 0) # 添加Tukey标签
# 显示图形
print(p)
# 保存图形为PNG文件
ggsave("boxplot_Simpson_diversity.png", plot = p, width = 10, height = 8, dpi = 600, bg = "white")