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

最新推荐文章于 2023-08-29 02:51:26 发布

原创最新推荐文章于 2023-08-29 02:51:26 发布 · 5k 阅读

1 ·

CC 4.0 BY-SA版权

文章标签：

#R语言

R语言专栏收录该内容

5 篇文章

订阅专栏

探讨了R语言中处理浮点数时遇到的精度问题，包括比较和运算中的异常结果，并提供了几个典型示例。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

举个例子：

> 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里面的数据尤其是判断时，需要特别注意精度问题。

解决方案有待补充。