多处理组植物种群Shannon指数计算的技术求助

先处理数据清洗：修正处理组名称 #

你的示例数据中Treatment列包含序号前缀（比如"1 b"），先提取纯处理标识，避免分组混乱：

# 提取Treatment列的字母部分，生成干净的分组列
d$Treatment_clean <- gsub("^\\d+ ", "", d$Treatment)

方法一：用vegan包diversity函数批量计算 #

1. 逐行计算单样本Shannon指数 #

先筛选出所有植物物种列，再逐行计算：

library(vegan)
# 选择所有植物频度列（排除原Treatment和新生成的Treatment_clean列）
species_data <- d[, !names(d) %in% c("Treatment", "Treatment_clean")]
# 为每行添加Shannon指数新变量
d$Shannon <- apply(species_data, 1, function(row) diversity(row, index = "shannon"))

2. 按处理组汇总计算组水平Shannon指数 #

如果需要统计每个处理组的整体种群多样性：

# 按处理组分组，计算每组各物种的频度总和
grouped_species <- aggregate(. ~ Treatment_clean, data = species_data, sum)
# 为每组添加汇总后的Shannon指数
grouped_species$Shannon_group <- apply(grouped_species[, -1], 1, diversity, index = "shannon")

方法二：手动实现Shannon指数计算 #

不依赖vegan包，自定义函数完成计算（自动规避log(0)报错）：

# 定义Shannon指数计算函数，过滤0频度物种
calc_shannon <- function(pop) {
  pop <- pop[pop > 0] # 移除频度为0的物种
  total <- sum(pop)
  prop <- pop / total
  -sum(prop * log(prop))
}

# 逐行计算单样本指数
d$Shannon_manual <- apply(species_data, 1, calc_shannon)

# 按处理组汇总计算
grouped_manual <- aggregate(. ~ Treatment_clean, data = species_data, sum)
grouped_manual$Shannon_group_manual <- apply(grouped_manual[, -1], 1, calc_shannon)

问题排查说明 #

• 之前使用diversity失败，大概率是没正确筛选物种列，或者Treatment列的序号前缀导致分组逻辑混乱，先清洗分组列是核心前提。
• 手动计算时必须过滤0频度物种，否则会因log(0)出现NaN结果。

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