为何R代码计算PAR达最大值0.01%时的Pressure返回NA？

问题描述 #

现有包含PAR与Pressure字段的数据集，已知最大PAR值为588.46（对应Pressure=3），需求是查找**PAR≤最大值0.01%**时的最小Pressure。数据中Pressure=54时PAR=5.72符合条件，但运行提供的R代码却返回NA。代码在多数文件中正常运行，且已设置na.rm=TRUE，需排查原因。

附原数据集及代码：

test <- structure(list(PAR = c(79, 374.99, NA, 588.46, 390.78, 351.5, 
298.24, 246.64, 240.41, NA, 189.88, 167.27, 163.61, 152.5, 131.08, 
156.68, 139.7, 121.61, 107.22, 96.73, 89.56, 83.34, 77.06, 70.82, 
65.65, 67.04, 54.24, 51.22, 49.62, 46.63, 44.26, 42.66, 40.75, 
39.72, 29.53, 31.6, 30.68, 28.2, 25.73, 23.77, 21.84, 19.94, 
18.08, 16.39, 14.79, 13.44, 12.31, 11.3, 10.34, 9.38, 8.61, 8, 
7.41, 6.54, 5.72, 5.37, 5.02, 4.69, 4.36, 4.04, 3.75, 3.47, 3.22, 
3, 2.81, 2.64, 2.48, 2.35, 2.22, 2.1, 1.98, 1.87, 1.76, 1.66, 
1.57, 1.48, 1.39, 1.32, 1.24, 1.23, 1.11, 1.07, 1.01, 0.95, 0.89, 
0.85, 0.8, 0.76, 0.72, 0.7, 0.66, 0.64, 0.61, 0.57, 0.55, 0.52, 
0.49, 0.47, 0.44, 0.4), Pressure = c(0, 1, 2, 3, 4, 5, 6, 7, 
8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 
24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 
40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 
56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 
72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 
88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99)), row.names = c(NA, 
100L), class = "data.frame")

miss <- function(x) ifelse(is.finite(x),x,NA)
d.depths <- test %>% 
summarize(
euphotic=miss(min(Pressure[PAR <= 0.01*max(PAR)], Inf, na.rm=TRUE)))

原因分析 #

• 核心问题：max(PAR)未处理NA值：数据中存在NA的PAR值，原代码里max(PAR)没有设置na.rm=TRUE，导致max(PAR)返回NA，进而0.01*max(PAR)也为NA。此时PAR <= NA的判断结果全为NA，min()函数处理全NA的向量时，即便设置了na.rm=TRUE，也会返回NA。
• 冗余函数miss()：min()配合na.rm=TRUE已经能返回有效有限值，额外的miss()函数完全多余，甚至可能干扰结果。

修正方案 #

修改代码，在max(PAR)中添加na.rm=TRUE，并移除冗余的miss()函数：

d.depths <- test %>% 
  summarize(
    euphotic = min(Pressure[PAR <= 0.01 * max(PAR, na.rm = TRUE)], na.rm = TRUE)
  )

修正说明 #

1. max(PAR, na.rm = TRUE)会正确计算出PAR的最大值588.46，0.01%的阈值即为5.8846，此时PAR=5.72（对应Pressure=54）满足条件。
2. 移除miss()函数后，min()直接返回符合条件的最小Pressure值54，结果正确。

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