3d NMDS多样性分析图 R语言

CAAS_IFR_zp

已于 2024-10-16 10:19:54 修改

阅读量586

点赞数 4

文章标签： r语言开发语言

于 2024-10-16 10:17:55 首次发布

本文链接：https://blog.csdn.net/m0_53945548/article/details/142972844

版权

# 安装并加载必要的包
if (!require("vegan")) install.packages("vegan")
if (!require("ggplot2")) install.packages("ggplot2")
if (!require("plotly")) install.packages("plotly")
if (!require("reticulate")) install.packages("reticulate")

# 加载包
library(vegan)
library(ggplot2)
library(plotly)
library(reticulate)

# 设置工作目录
setwd("C:/Users/fordata/Desktop/研究生/第二个想法(16s肠型＋宏基因组功能)/第二篇病毒组/result/β多样性")

# 载入数据
tpm <- read.table("fetal_2.0_rpkm.txt", header = TRUE, row.names = 1)
tpm <- t(tpm)
metadata <- read.table("fetal_2.0_metadata.txt", header = TRUE, row.names = 1)
metadata$group <- as.factor(metadata$group)

# 计算 Bray-Curtis 距离矩阵
dist_matrix <- vegdist(tpm, method = "bray")

# 执行 NMDS 分析
nmds_result <- metaMDS(dist_matrix, k = 3, trymax = 200)  # 增加 trymax 的值

# 提取NMDS的坐标得分（提取样本点的坐标）
scores_nmds <- as.data.frame(nmds_result$points)
colnames(scores_nmds) <- c("NMDS1", "NMDS2", "NMDS3")

# 将元数据添加到结果中
df_nmds <- cbind(scores_nmds, metadata)

# 查看列名 (确保有 NMDS1, NMDS2, NMDS3)
colnames(df_nmds)

# 绘制3D NMDS图
fig <- plot_ly(df_nmds, 
               x = ~NMDS1, 
               y = ~NMDS2, 
               z = ~NMDS3, 
               color = ~group, 
               colors = c("AF" = "#BD3C29", "M" = "#0172B6", "P" = "#78D3AC", "UB"="#E18727", "NC"="black"),
               type = 'scatter3d', 
               mode = 'markers', 
               marker = list(size = 13)) %>%
  layout(scene = list(
    xaxis = list(title = "NMDS1"),
    yaxis = list(title = "NMDS2"),
    zaxis = list(title = "NMDS3")
  ),
  legend = list(x = 1, y = 0.5, font = list(size = 16)))  # 调整图例位置

# 显示图形
fig

# 计算组间差异的PERMANOVA分析
adonis_result <- adonis2(dist_matrix ~ group, data = metadata, permutations = 999)
cat("PERMANOVA结果:\n")
print(adonis_result)

# 计算组间差异的ANOSIM分析
anosim_result <- anosim(dist_matrix, metadata$group, permutations = 999)
cat("\nANOSIM结果:\n")
print(anosim_result)

# 保存为高分辨率 PDF
# 使用外部工具或 `orca` 保存为 PDF（例如需要安装 `orca` 工具）
# fig %>% save_image("fetal_nmds_3d_plot.pdf", width = 800, height = 600, scale = 2)

保存图片是个问题，我试了一下，export出html文件，然后到网页截图会比较清楚

PERMANOVA 结果（通过 adonis2() 计算）：你将会得到组间差异的 R² 值、F 值以及 p 值。p 值越低，组间差异越显著。
ANOSIM 结果（通过 anosim() 计算）：你将会得到组间相似性的 R 值和相应的 p 值。R 值越接近 1，组间差异越大；R 值越接近 0，说明组内和组间的差异相似。