SAS以增量形式循环运行

Question

为什么以下内容不能产生期望的结果？

data data1; 

do a = 0.0 to 1.0 by 0.1; 
    do b = 0.0 to 1.0 by 0.1; 
     do c = 0.0 to 1.0 by 0.1; 
      do d = 0.0 to 1.0 by 0.1; 
       if (a+b+c+d)=1 then output; 
      end; 
     end; 
    end; 
end; 
format a b c d 4.1; 
run; 

Answer 1

我不是很熟悉SAS，但总体上，当你有号码，如0.1，它在二进制文件的表示。由于.1不能完全用二进制表示，所以数学方程并不总是完全相加。例如，在浮点运算中，0.1乘以10不等于1.0。通常，在使用浮点时不要使用相等性。

参见http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html

Answer 2

号码在SAS被存储在二进制（因为它们是在大多数计算应用），并且常常不能精确地表示相应的十进制数。就像0.33333333！=（1/3）一样，许多1/10也不代表它们的精确十进制值。

data data1; 

do a = 0.0 to 1.0 by 0.1; 
    do b = 0.0 to 1.0 by 0.1; 
     do c = 0.0 to 1.0 by 0.1; 
      do d = 0.0 to 1.0 by 0.1; 
       if round(a+b+c+d,.1)=1 then output; 
      end; 
     end; 
    end; 
end; 
format a b c d 4.1; 
run; 

舍入解决了这个问题。

有关更多信息，您可以阅读this SAS technical paper。