在R中计算多PAHs总浓度的随机误差传播方法咨询

你已经完成了单样本PAHs各污染物均值与标准差的计算，现在要基于公式σx = sqrt(σa² + σb² + ... + σn²)计算总浓度的误差传播值，无需使用复杂的propagate包，直接用dplyr即可实现：

实现代码 #

# 在已有means数据框基础上，同时计算总浓度和对应误差
summary_with_error <- means %>% 
  group_by(sample_ID) %>% 
  summarise(
    sum_pollutant = sum(mean_conc),
    sum_error = sqrt(sum(sd_conc^2))  # 应用误差传播公式
  )

# 查看结果
summary_with_error

结果说明 #

运行上述代码后，summary_with_error会包含两列：

• sum_pollutant：16种PAHs的浓度均值之和（即你之前得到的21.5）
• sum_error：总浓度的误差传播值，通过对每种污染物的标准差平方求和后开根号得到

针对你的示例数据，计算出的sum_error约为0.97（实际值以代码运行结果为准）。

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