R: 阿尔法α多样性计算和箱图制作,以及差异分析

R 复制代码
# install.packages("vegan")
library(vegan)
library(ggplot2)
library(ggpubr)

setwd("xxx")
# 使用read.table()函数读取数据
df <- read.table("xxx", header = TRUE, row.names = 1)

# 转置数据框
df <- t(df)

# 计算每个样品的香农多样性指数
shannon <- apply(df, 1, diversity)

# 创建一个数据框来存储样品名和对应的香农多样性指数
data <- data.frame(sample = names(shannon), shannon = shannon)

# 读取元数据
metadata <- read.table("mother_metadata.txt", header = TRUE, row.names = 1)

# 将元数据和样品数据合并到一个数据框中
data <- merge(data, metadata, by = "row.names", all.x = TRUE)

# 设置每个组的颜色
group_colors <- c("IC" = "#BD3C29", "RC" = "#0172B6", "OS" = "#78D3AC", "VS" = "#E18727")
# 首先设置比较的列表
compare_list <- list(
  c("RC","VS"))

# 使用ggplot2创建箱线图
p <- ggplot(data, aes(x = group, y = shannon, fill = group, colour = group)) +
  geom_boxplot(width = 0.5, alpha = 0.6, lwd = 1.15, outlier.shape = NA) +  # 调整箱的大小
  geom_jitter(width = 0.3, size = 3, alpha = 0.75) +  # 添加散点
  labs(y = "Shannon Diversity") +
  scale_fill_manual(values = group_colors) +  # 设置颜色
  scale_color_manual(values = group_colors) +
  theme_minimal() +
  labs(x = NULL) +
  theme(panel.border = element_rect(colour = "black", fill=NA, linewidth = 1.1)) +
  stat_compare_means(
    comparisons = compare_list,
    method = "wilcox.test",
    label = "p.signif",
    hide.ns = TRUE)
  # 添加检验结果

# 显示图形
#print(p)

# 保存图形为PDF
ggsave("barplot.pdf", p, height = 5, width = 5)
相关推荐
czliutz1 天前
R绘制股票日波动线图 中国海油600938
开发语言·r语言·r语言绘图
Biomamba生信基地2 天前
R语言绘图神器| ggplot2与其基本用法介绍
开发语言·信息可视化·r语言·生信·医药
天桥下的卖艺者3 天前
R语言基于shiny开发随机森林预测模型交互式 Web 应用程序(应用程序)
前端·随机森林·r语言·shiny
小火柴1233 天前
利用R语言绘制直方图
开发语言·r语言
zhangfeng11333 天前
R语言 读取tsv的三种方法 ,带有注释的tsv文件
开发语言·r语言·生物信息
Tiger Z4 天前
《R for Data Science (2e)》免费中文翻译 (第9章) --- Layers(2)
r语言·数据科学·中文翻译
海大哥4 天前
深入解析Java中HashMap的工作原理与性能优化策略
r语言
万粉变现经纪人5 天前
如何解决 pip install -r requirements.txt 约束文件 constraints.txt 仅允许固定版本(未锁定报错)问题
开发语言·python·r语言·django·beautifulsoup·pandas·pip
万粉变现经纪人5 天前
如何解决 pip install -r requirements.txt 无效可编辑项 ‘e .‘(-e 拼写错误)问题
开发语言·python·r语言·beautifulsoup·pandas·pip·scipy