如何使用R语言Vegan包计算全数据集的Gamma多样性Shannon H指数

解决方案 #

方法1：直接计算Gamma多样性，无需生成汇总数据集 #

你不需要额外生成单站点汇总表，直接对所有物种的汇总丰度调用diversity()函数即可，一行代码就能完成：

# 计算基于物种总丰度的Gamma Shannon指数（主流通用的计算逻辑）
vegan::diversity(colSums(df[, -1]), index = "shannon")

# 如果你的分析场景确实需要按平均丰度计算，替换为colMeans即可
vegan::diversity(colMeans(df[, -1]), index = "shannon")

colSums/colMeans会直接返回所有物种列的汇总结果，属于一维向量，vegan::diversity()可以直接识别该输入为单样本的丰度向量，输出的结果就是你需要的单一Gamma多样性值。

方法2：高效生成单站点汇总数据集 #

如果你确实需要生成示例中的单站点汇总表，用base R或者dplyr都可以快速完成，不需要手动统计：

base R版本

# 生成平均丰度的汇总表
gamma_df <- data.frame(site = 1, t(colMeans(df[, -1])))

# 生成总丰度的汇总表把colMeans换成colSums即可

dplyr版本

library(dplyr)
gamma_df <- df %>%
  summarise(across(starts_with("Spec"), mean)) %>% # 自动匹配所有物种列，不用手动写列名
  mutate(site = 1, .before = 1)

注意点 #

常规Gamma多样性计算通常使用物种的总丰度而非平均丰度，避免站点数量差异对结果造成干扰，你可以根据自己的分析场景选择对应的汇总逻辑。

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