如何从vegan包adonis结果提取R2值并批量导出至文件？

绝对可以实现！手动复制R²值这种繁琐操作完全可以避免，不管是单次分析还是批量迭代运行adonis，都能通过代码自动提取结果并导出到文件。我来给你一步步演示：

首先，vegan包的adonis（或更新的adonis2）函数返回的结果是一个列表对象，其中核心的方差分析表存放在aov.tab元素中，包含了R²、p值、自由度等所有关键统计量。

示例代码：

# 加载vegan包并准备示例数据
library(vegan)
data(dune)  # 植被群落数据
data(dune.env)  # 环境因子数据

# 运行adonis2分析（推荐用adonis2，结果结构更清晰）
adonis_res <- adonis2(dune ~ Management, data = dune.env)

# 提取完整的结果表格
result_df <- adonis_res$aov.tab

# 如果只需要提取R²值，直接调用列名即可
r2_values <- result_df$R2
# 提取p值同理
p_values <- result_df$`Pr(>F)`

当需要批量处理多个因子或多组数据时，用列表存储每次的adonis结果是最高效的方式，之后可以把列表中的结果合并成一个统一的数据框：

# 定义需要测试的环境因子列表
target_factors <- c("Management", "A1", "Moisture")

# 创建空列表用于存储每次的adonis结果
adonis_result_list <- list()

# 循环运行adonis分析
for (factor in target_factors) {
  # 动态构建公式
  analysis_formula <- as.formula(paste("dune ~", factor))
  # 运行分析并存储结果表格
  adonis_result_list[[factor]] <- adonis2(analysis_formula, data = dune.env)$aov.tab
}

# 合并列表中的所有结果为一个数据框（保留因子名称作为标识）
# 方法1：用tidyverse工具
library(dplyr)
library(tibble)

combined_results <- bind_rows(adonis_result_list, .id = "Tested_Factor")

# 方法2：用基础R实现（无需额外包）
combined_results <- do.call(rbind, lapply(names(adonis_result_list), function(factor_name) {
  cbind(Tested_Factor = factor_name, adonis_result_list[[factor_name]])
}))

最后用write.table或write.csv把整理好的结果保存到本地文件：

# 用write.table保存为制表符分隔的文本文件
write.table(combined_results, 
            file = "adonis_batch_results.txt", 
            sep = "\t", 
            row.names = FALSE, 
            quote = FALSE)

# 或者用write.csv保存为CSV文件（更适合Excel打开）
write.csv(combined_results, 
          file = "adonis_batch_results.csv", 
          row.names = FALSE)

这样你就可以直接在Excel或其他工具中打开导出的文件，完全不需要手动复制粘贴啦！

内容的提问来源于stack exchange，提问作者user7886302