R在浮点运算中的精度问题

    举个例子：

> x = seq(0,1,by = 0.2)
> y = seq(0,1,by = 0.2)
> x[3]
<span style="color:#999900;">[1] 0.4</span>
> y[4]
<span style="color:#999900;">[1] 0.6</span>
> 1 - x[3]
<span style="color:#999900;">[1] 0.6</span>
> y[4] > 1-x[3]
<span style="color:#999900;">[1] TRUE</span>

    
    这里竟然出现了0.6>0.6的情况。这样的结果应算是R的缺陷，一般的数学运算软件都会通过一些方法来克服这个问题。简单的说，浮点数在计算机中的表达有限制，不能以任意精度存储，尤其是微小的或巨大的数字在运算时经常会有些意外情况。

    再举几个例子：

> seq(0,1,0.1)
<span style="color:#999900;">[1] 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0</span>

> seq(0,1,0.1) == c(0.0, 0.1, 0.2, 0.3, 0.4 ,0.5, 0.6, 0.7, 0.8, 0.9, 1.0)
<span style="color:#999900;">[1]  TRUE  TRUE  TRUE FALSE  TRUE  TRUE FALSE FALSE  TRUE  TRUE  TRUE</span>
> 0.3-0.7+0.4 == 0
<span style="color:#999900;">[1] FALSE</span>
> sqrt(2)^2 == 2
<span style="color:#999900;">[1] FALSE</span>

    因此，在处理R里面的数据尤其是判断时，需要特别注意精度问题。

    解决方案有待补充。