SAS PROC SQL中Case When数值比较异常原因咨询

你遇到的这个问题本质上是二进制浮点数无法精确表示部分十进制小数导致的，这是所有使用IEEE 754标准浮点数的编程语言/工具（包括SAS）都会遇到的共性问题。

为什么0.015会引发异常？ #

十进制的0.015转换成二进制是一个无限循环的小数：

0.015₁₀ = 0.0000011110101110000101000111101011100001010001111...₂

SAS默认使用64位双精度浮点数，只能存储有限位的二进制小数，所以程序中直接写的.015其实是它的近似值，和理论上的十进制0.015存在微小误差。当你计算MOD1 - MOD0（比如0.955 - 0.94）时，得到的结果看起来是0.015，但实际上也是一个接近0.015的近似值——这两个近似值并不完全相等，所以NOT = .015的判断会返回真，触发错误标记。

为什么1 - .985能正常运行，.2 - .185却不行？ #

这是因为浮点数计算的误差抵消情况不同：

• 计算1 - .985时，两个数的存储误差刚好相互抵消，最终得到的近似值和SAS存储的.015近似值完全一致，所以相等判断能正确返回。
• 而.2 - .185的计算误差和.015的存储误差无法抵消，两个近似值仍存在细微差异，判断自然会失效。

验证这个误差的方法 #

你可以用SAS的PUT函数以高精度格式输出这些值，就能看到背后的真实存储内容：

data float_check;
  mod0 = .94;
  mod1 = .955;
  calc_diff = mod1 - mod0;
  target_val = .015;
  /* 用32位精度输出，显示浮点数的真实存储值 */
  put calc_diff= BEST32. target_val= BEST32.;
run;

运行后你会看到，calc_diff和target_val的实际值存在极细微的差别，比如可能是0.015000000000000003和0.014999999999999999，正是这点差异导致了直接相等判断失败。

通用的解决思路 #

处理浮点数比较的正确做法是不要直接判断相等/不等，而是判断差值是否在一个极小的容忍范围内，比如：

CASE 
  WHEN ABS(MOD1 - MOD0 - .015) > 1e-8 THEN 2 
  WHEN ABS(MOD2 - MOD1 - .015) > 1e-8 THEN 3 
  WHEN ABS(MOD3 - MOD2 - .015) > 1e-8 THEN 4 
  WHEN ABS(MOD4 - MOD3 - .015) > 1e-8 THEN 5 
  ELSE 0 
END

这里的1e-8是可根据数据精度调整的容忍阈值，只要差值小于这个值，就认为二者“相等”。

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