如何用R将地球化学数据的含量与误差列合并为“值±误差”列？

需求可行性说明 #

完全可行，不管是单元素列合并，还是批量处理200+变量的大量数据，都可以通过R代码高效实现，且能完美适配你使用xlsx包生成正式报告的需求。

具体实现步骤 #

1. 单元素列合并（示例：镁元素） #

如果只需要处理单个元素的列对，直接使用paste()函数将含量、±符号、误差拼接成字符串，之后删除原误差列即可：

# 合并镁的含量与误差列
PPI10 <- PPI10 %>%
  mutate(`Magnesium Abundance [ppm]` = paste(`Magnesium Abundance [ppm]`, "±", `Magnesium error [ppm]`, sep = " ")) %>%
  select(-`Magnesium error [ppm]`) # 删除多余的误差列

2. 批量处理所有元素（适配200+变量场景） #

针对大量元素的批量处理，可通过匹配列名规则自动配对“含量列”与“误差列”，无需逐个手动处理：

# 筛选所有含"Abundance"的含量列
abundance_cols <- grep("Abundance", colnames(PPI10), value = TRUE)

# 循环处理每一对含量-误差列
for (col in abundance_cols) {
  # 提取元素名称，构造对应的误差列名
  element_name <- sub(" Abundance \\[ppm\\]", "", col)
  error_col <- paste(element_name, "error [ppm]", sep = " ")
  
  # 合并列并删除原误差列（使用tidyverse的非标准求值确保列名正确识别）
  PPI10 <- PPI10 %>%
    mutate(!!col := paste(!!sym(col), "±", !!sym(error_col), sep = " ")) %>%
    select(-!!sym(error_col))
}

3. 整合到现有代码流程 #

将上述合并代码插入到数据清理完成后，写入Excel之前的位置，也就是现有代码中PPI10 <- PPI10.pXRF[PPI10.pXRF$Method Name != "Cal Check", ]语句之后，addDataFrame(PPI10, sheet, ...)语句之前。

4. 验证处理结果 #

可以通过以下代码查看处理后的列内容，确认格式正确：

# 查看镁元素列的处理结果示例
head(PPI10$`Magnesium Abundance [ppm]`)

关键注意事项 #

• 列名匹配规则：确保误差列命名严格遵循“元素名 error [ppm]”格式，与含量列的“元素名 Abundance [ppm]”对应；如果你的列名有其他规则，需调整sub()函数的字符串匹配逻辑。
• ±符号显示：直接使用中文输入法输入的±或Unicode编码"\u00B1"，都能在Excel中正常显示。
• 数据类型：合并后的列会转为字符型，完全满足报告展示需求，不影响Excel中的格式渲染。

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