R语言Dataframe转换：构建物种-站点存在/缺失矩阵求助

我明白你要做的是把这种按站点列、每行存放物种的长格式数据，转换成以站点为行、所有不重复物种为列的0/1存在矩阵对吧？这是生态数据处理里很常见的格式转换需求，我给你两种实现方式，一种用tidyverse（更简洁直观），一种用base R（不需要额外包），你可以按需选择。

首先先修正一下你提供的模拟数据代码（原来的cbind会生成矩阵，改成c更贴合后续的dataframe操作）：

# 模拟原始数据
A1 <- c("sp1","sp2","sp3", "sp4", "sp7", "sp8") 
A2 <- c("sp1","sp3", "sp4", "sp7", "sp9") 
A3 <- c("sp5","sp6","sp7", "sp10") 
A4 <- c("sp1","sp2","sp7", "sp9", "sp10") 
A5 <- c("sp3","sp4") 

max_row <- 6
# 补NA到统一长度
A1 <- c(A1, rep(NA, max_row - length(A1)))
A2 <- c(A2, rep(NA, max_row - length(A2)))
A3 <- c(A3, rep(NA, max_row - length(A3)))
A4 <- c(A4, rep(NA, max_row - length(A4)))
A5 <- c(A5, rep(NA, max_row - length(A5)))

# 生成原始Dataframe
df <- data.frame(A1, A2, A3, A4, A5, stringsAsFactors = FALSE)

方法一：用tidyverse包（推荐） #

tidyverse里的tidyr包专门处理这种宽长格式转换，代码逻辑清晰，容易理解：

library(tidyverse)

# 1. 把宽格式转成长格式，提取所有站点-物种的有效组合
long_df <- df %>%
  pivot_longer(cols = everything(), names_to = "站点", values_to = "物种") %>%
  filter(!is.na(物种)) # 去掉填充的NA值

# 2. 转成目标的宽格式，生成0/1存在矩阵
result_df <- long_df %>%
  pivot_wider(
    names_from = 物种, 
    values_from = 物种, 
    values_fn = ~1,  # 存在的物种标记为1
    values_fill = 0  # 缺失的物种标记为0
  )

# 查看最终结果
print(result_df)

代码解释： #

• pivot_longer：把原来的A1-A5这些列名转换成“站点”列的取值，每个站点下的物种都拆成单独的行，这样就得到了所有真实存在的站点-物种配对。
• filter(!is.na(物种))：过滤掉之前为了统一行数填充的NA，这些不是真实的物种数据。
• pivot_wider：把物种名转成列，只要站点存在该物种就填1，不存在就填0，完美匹配你的需求。

方法二：用Base R（无需额外包） #

如果你不想安装tidyverse，用Base R也能实现：

# 1. 提取所有不重复的物种（去掉NA）
all_species <- unique(unlist(df))
all_species <- all_species[!is.na(all_species)]

# 2. 遍历每个站点，生成对应物种的0/1向量
site_species_matrix <- lapply(df, function(site_species) {
  sapply(all_species, function(sp) as.integer(sp %in% site_species))
})

# 3. 转换为Dataframe并添加站点列
result_df_base <- data.frame(
  站点 = names(site_species_matrix),
  do.call(rbind, site_species_matrix),
  row.names = NULL
)

# 查看结果
print(result_df_base)

代码解释： #

• 先收集所有出现过的不重复物种，作为最终矩阵的列名。
• 对每个站点列，逐个检查每个物种是否存在，用as.integer把逻辑值转成1/0。
• 最后把所有站点的结果合并成Dataframe，加上站点名称列。

两种方法得到的结果是完全一致的，你可以根据自己的习惯选择。如果你的原始数据还有其他特殊情况（比如物种名有重复、站点名特殊等），可以再调整代码细节。

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