R中相对丰度列和不为1但Excel中正常的技术问题

核心问题 #

你计算的物种相对丰度列和与预期严重不符（R中结果偏差极大，Excel中显示为1），调整容差无效——从你提供的sums_relative_abundance结果看，部分列和仅0.5左右，显然不是浮点精度问题，而是数据逻辑或计算流程出错。

可能原因及排查步骤 #

1. 行列含义搞反（最常见） #

物种相对丰度通常是每个样本的所有物种丰度归一化到和为1：

• 若你的数据是「行=物种，列=样本」，则列和应为1；
• 若数据是「行=样本，列=物种」，则行和应为1，此时计算列和必然不符合预期。

排查代码：

# 查看数据维度，确认行列对应关系
dim(relative_abundances)
# 检查行和是否接近1
row_sums <- rowSums(relative_abundances)
head(row_sums)
all(abs(row_sums - 1) < 1e-6)

2. 相对丰度计算逻辑错误 #

确认你计算相对丰度时是否按正确维度归一化：

• 正确的列归一化（行=物种，列=样本）代码：

# 假设raw_data是原始丰度数据
relative_abundances <- apply(raw_data, 2, function(x) x / sum(x, na.rm = TRUE))

• 若误按行归一化，会导致列和完全偏离1：

# 错误示例：按行归一化，此时行和为1，列和无意义
relative_abundances <- apply(raw_data, 1, function(x) x / sum(x, na.rm = TRUE))

3. 数据导入异常 #

检查数据是否存在非数值列或缺失值：

# 查看数据结构，确认所有列都是数值型
str(relative_abundances)
# 检查是否有缺失值
any(is.na(relative_abundances))

如果存在非数值列，colSums会自动跳过或报错，导致计算结果异常；若有缺失值，需确认sum时是否忽略了NA（Excel默认将NA视为0，而R中sum遇到NA会返回NA，需用na.rm=TRUE）。

4. Excel计算的视觉误差 #

Excel默认会对数字进行四舍五入显示，你手动计算的「和为1」可能是显示值，而非实际存储值。可在Excel中查看单元格的完整小数位，或用=SUM(A:A)公式确认精确值。

验证后的修正方案 #

如果排查后确认是行列搞反，可转置数据后再验证列和：

# 转置数据，使样本为列、物种为行
transposed_data <- t(relative_abundances)
# 计算转置后的列和（即原数据的行和）
sums_corrected <- colSums(transposed_data)
all(abs(sums_corrected - 1) < 1e-6)

如果是相对丰度计算错误，重新按正确维度执行归一化即可。

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