在R语言中统计物种-站点矩阵中各物种的出现次数

我完全懂你的需求——你要统计的是每个物种在站点中的出现频次（只要在某个站点里存在，不管丰度多少都算1次），而不是直接用colSums统计的丰度总和对吧？这个问题其实很容易解决，核心思路是先把原矩阵转换成「存在/不存在」的二进制矩阵，再用求和函数统计。

具体步骤示例 #

假设你的物种-站点矩阵是这样的（和你描述的结构匹配）：

# 构造示例矩阵
species_matrix <- matrix(
  c(5, 0, 2, 1,
    3, 4, 0, 0,
    1, 0, 0, 5,
    2, 3, 0, 4),
  nrow = 4, ncol = 4,
  dimnames = list(c("W", "X", "Y", "Z"), c("a", "b", "c", "d"))
)

接下来只需要两步：

1. 把矩阵中所有大于0的元素替换成1（代表物种存在），0保持不变（代表物种未出现）
2. 用colSums统计每列的和，就是每个物种的出现次数

代码实现：

# 转换成存在/不存在二进制矩阵
presence_absence <- ifelse(species_matrix > 0, 1, 0)

# 统计每个物种的出现次数
species_occurrence <- colSums(presence_absence)

# 查看结果
print(species_occurrence)

运行后你会得到和预期一致的结果：

a b c d 
4 2 1 2 

简化写法 #

如果你不想单独保存中间矩阵，也可以一行代码搞定：

species_occurrence <- colSums(ifelse(species_matrix > 0, 1, 0))

这个方法不管你的矩阵是基础的matrix类型还是data.frame类型都适用，非常灵活～

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