如何为R包vegan配置α多样性分析的计数数据

数据格式适用性 #

你的数据格式完全适配vegan包的分析要求：vegan的多样性分析函数默认要求行代表样本（你已将同一年份的站点数据汇总为单个样本），列代表物种，单元格为物种计数数据，你的矩阵结构完全匹配该规则。仅需注意将Year列从计数矩阵中分离，作为样本标识（行名）即可。

完整操作步骤 #

1. 数据导入与预处理 #

先将Excel文件另存为CSV格式（规避Excel格式的导入兼容性问题），再在R中执行以下操作：

# 首次使用需安装vegan包
install.packages("vegan")
library(vegan)

# 读取CSV数据，注意提前清理Year列的多余符号（比如2022.改为2022）
data <- read.csv("你的数据文件名.csv", stringsAsFactors = FALSE)

# 将Year列设为行名，移除原Year列，保留纯物种计数矩阵
rownames(data) <- data$Year
data <- data[, -which(colnames(data) == "Year")]

# 检查数据结构，确保所有列均为数值型计数
str(data)

常见导入问题解决：

• 若Year列被识别为字符串，可手动在Excel中改为数值型，或用data$Year <- as.integer(gsub("\\.", "", data$Year))清理格式
• 若物种列存在非数值内容，用data <- data %>% mutate(across(everything(), as.numeric))（需加载dplyr包）转换类型

2. α多样性指标计算 #

物种丰富度（现存物种数）

# 计算各年份的物种丰富度
richness <- specnumber(data)
# 转为数据框方便查看
richness_df <- data.frame(Year = names(richness), Richness = richness)
print(richness_df)

Shannon-Wiener指数

# 计算Shannon指数，默认以自然对数为底，需以2为底可加参数base=2
shannon <- diversity(data, index = "shannon")
shannon_df <- data.frame(Year = names(shannon), Shannon = shannon)
print(shannon_df)

合并所有指标

alpha_div <- merge(richness_df, shannon_df, by = "Year")
print(alpha_div)

3. 丰度曲线（Rank-Abundance Curve） #

按年份绘制物种丰度排序曲线：

# 设置多图布局，按年份数量分配面板
par(mfrow = c(1, length(rownames(data))))
for (year in rownames(data)) {
  # 提取该年份非0丰度的物种并降序排序
  abundances <- sort(data[year, data[year,] > 0], decreasing = TRUE)
  # 绘制曲线，y轴取对数更易观察差异
  plot(abundances, type = "l", main = year, xlab = "物种排序", ylab = "物种丰度", log = "y")
}
par(mfrow = c(1,1)) # 重置绘图布局

也可使用vegan的radfit()函数拟合丰度分布模型：

rad_model <- radfit(data)
plot(rad_model)

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