R语言计算含0值的相对差异：解决Inf值问题

首先明确你需要的计算逻辑：Case相对于Control的倍数百分比，即(Case / Control) * 100，对应第一例AB的结果是400%。针对Control为0导致的Inf问题，有两种常用解决方案：

方案1：标记无意义结果（推荐） #

当Control为0时，相对差异的计算逻辑不成立，直接标记为N/A或NA，避免误导。

代码实现（Base R） #

# 加载数据
dat <- data.frame(Var = c("AB", "TG", "GF", "ZH", "ED", "MS"),
                  Case = c(0.4, 3.2, 0.23, 14, 7.1, 5.5),
                  Control = c(0.1, 6.3, 0.18, 14, 0, 5.6))

# 新增相对差异列，Control为0时标记为N/A，其余转为百分比格式
dat$RelDiff <- ifelse(dat$Control == 0, 
                      "N/A", 
                      sprintf("%.0f%%", (dat$Case / dat$Control) * 100))

# 查看结果
print(dat)

运行结果 #

Var  Case Control RelDiff
1  AB  0.40    0.10    400%
2  TG  3.20    6.30     51%
3  GF  0.23    0.18    128%
4  ZH 14.00   14.00    100%
5  ED  7.10    0.00     N/A
6  MS  5.50    5.60     98%

如果需要数值型结果而非字符型百分比，可改为：

dat$RelDiff_Num <- ifelse(dat$Control == 0, 
                          NA_real_, 
                          (dat$Case / dat$Control) * 100)

方案2：替换0为极小值 #

若必须得到数值型结果，可将Control列的0替换为一个极小值（如1e-6），避免除以0。但需注意，这种情况下Control为0对应的结果会是极端大值，后续分析时需单独处理。

代码实现 #

dat$RelDiff_Num_Alt <- (dat$Case / pmax(dat$Control, 1e-6)) * 100

说明 #

pmax(dat$Control, 1e-6)会将Control中小于等于1e-6的值替换为1e-6，确保分母不为0。ED行的结果会是7.1 / 1e-6 * 100 = 710000000%，需在报告中明确说明该值是替换0后的近似结果。

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