如何用R语言将结果与不确定度以科学计数法±形式合并并自动调整指数？

要实现将结果值和不确定度合并为类似0.20±1.52E-17的格式，核心是分别格式化结果值与不确定度后拼接；若需要统一两者的指数，可通过提取指数、调整数值来实现。

方案1：结果值用普通小数，不确定度自动调整为标准科学计数法 #

这种方式直接匹配你给出的示例格式，不确定度会自动调整指数，确保数值部分在1~10之间：

merge_result_uncert <- function(value, uncert, value_decimals = 2, uncert_decimals = 2) {
  # 格式化结果值为指定小数位数的普通数字
  formatted_value <- sprintf(paste0("%.", value_decimals, "f"), value)
  # 格式化不确定度为标准科学计数法，自动调整指数
  formatted_uncert <- sprintf(paste0("%.", uncert_decimals, "E"), uncert)
  # 拼接成最终字符串
  paste(formatted_value, formatted_uncert, sep = "±")
}

# 测试示例
merge_result_uncert(0.20, 1.52e-17)
# 输出："0.20±1.52E-17"

# 测试需要调整指数的情况：比如不确定度为1.52e-18
merge_result_uncert(0.20, 1.52e-18)
# 输出："0.20±1.52E-18"
# 若想将其调整为E-17（即数值部分变为0.15），可手动调整不确定度后再传入：
merge_result_uncert(0.20, 1.52e-18 * 10)
# 输出："0.20±1.52E-17"

方案2：统一结果值与不确定度的指数 #

如果需要让两者使用同一个指数（即实际值为a×10^exp ± b×10^exp，格式为a±bEexp），可按以下步骤实现：

merge_with_common_exp <- function(value, uncert, decimals = 2) {
  # 将结果值转为科学计数法，提取指数作为共同指数
  value_sci <- format(value, scientific = TRUE, digits = decimals)
  common_exp <- as.integer(sub(".*E([+-]\\d+)", "\\1", value_sci))
  
  # 提取结果值的数值部分
  value_num <- as.numeric(sub("E.*", "", value_sci))
  
  # 调整不确定度的数值部分，匹配共同指数
  uncert_sci <- format(uncert, scientific = TRUE, digits = decimals)
  uncert_exp <- as.integer(sub(".*E([+-]\\d+)", "\\1", uncert_sci))
  uncert_num <- as.numeric(sub("E.*", "", uncert_sci)) * 10^(uncert_exp - common_exp)
  
  # 拼接成最终格式
  sprintf(paste0("%.", decimals, "f±%.", decimals, "fE%d"), value_num, uncert_num, common_exp)
}

# 测试：结果值0.20（科学计数法为2.00E-01），不确定度1.52e-17
merge_with_common_exp(0.20, 1.52e-17)
# 输出："2.00±0.00E-01"

# 若想对齐到不确定度的指数，修改共同指数的来源即可：
merge_with_common_exp_uncert <- function(value, uncert, decimals = 2) {
  uncert_sci <- format(uncert, scientific = TRUE, digits = decimals)
  common_exp <- as.integer(sub(".*E([+-]\\d+)", "\\1", uncert_sci))
  
  uncert_num <- as.numeric(sub("E.*", "", uncert_sci))
  
  value_sci <- format(value, scientific = TRUE, digits = decimals)
  value_exp <- as.integer(sub(".*E([+-]\\d+)", "\\1", value_sci))
  value_num <- as.numeric(sub("E.*", "", value_sci)) * 10^(value_exp - common_exp)
  
  sprintf(paste0("%.", decimals, "f±%.", decimals, "fE%d"), value_num, uncert_num, common_exp)
}

# 测试
merge_with_common_exp_uncert(0.20, 1.52e-17)
# 输出："20000000000000000.00±1.52E-17"

说明 #

• 方案1适合结果值与不确定度数量级差异较大的场景，既保持结果值的可读性，又让不确定度以标准科学计数法呈现。
• 方案2适合需要严格统一指数的场景，可根据需求选择对齐结果值或不确定度的指数。

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