时间戳

首页 > 解决方案 > 使用 R：如何将两个数字与与选项（数字 = n）相关的浮点问题进行比较，以及数字是如何引入的？

r - 使用 R：如何将两个数字与与选项（数字 = n）相关的浮点问题进行比较，以及数字是如何引入的？

问题描述

因此，让我们从常规 R 设置中的正常值开始：

n1 = 0.15;
n2 = 0.15;

n1 == n2;  # we expect and get TRUE

现在，让我们使用选项将我们的数字更新为 22

options(digits=22);

n1 == n2; # we expect and still get TRUE but begin to worry

# our data has changed

> n1
[1] 0.14999999999999999
> n2
[1] 0.14999999999999999

我在运行 64 位的 Windoze 10 上。

现在让我们创建一个向量或导入一个 data.frame 或做一些理论上元素之一是 0.15 的事情......

vec =  seq(0,1, by=0.05);

vec;

 [1] 0.000000000000000000 0.050000000000000003 0.100000000000000006
 [4] 0.150000000000000022 0.200000000000000011 0.250000000000000000
 [7] 0.300000000000000044 0.350000000000000033 0.400000000000000022
[10] 0.450000000000000011 0.500000000000000000 0.550000000000000044
[13] 0.600000000000000089 0.650000000000000022 0.700000000000000067
[16] 0.750000000000000000 0.800000000000000044 0.850000000000000089
[19] 0.900000000000000022 0.950000000000000067 1.000000000000000000

最后，将向量与 n 之一进行比较

vec == n1;

 [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE
[13] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE

哦，快！经过几个小时的调试，我发现了这一点（在子集上下文中）。

有时，根据向量的构建方式，比较会正常工作。有时，就像在这个例子中一样，它不会。

> n1;
[1] 0.14999999999999999
> n2;
[1] 0.14999999999999999
> vec[4];
[1] 0.15000000000000002


> n1 == vec[4];
[1] FALSE
> n1 == n2;
[1] TRUE

如果我们可以应用 GLOBAL ENV 中的 Math.tolerance，它们是相等的。我的解决方法是手动建立一个公差并使用我自己的公差==用关节替换操作。>lowerbound | <upperbound

由于我认为这是由于浮点问题造成的，因此我希望从其他人那里获得一些见解。

x==y即使浮点问题表明它们可能不是，我如何比较两个我得到真实的数字？

标签： rfloating-pointprecisionfloating-accuracy

解决方案

all.equal在这里可能有用：

n1 <- 0.14999999999999999
n2 <- .15
n3 <- 0.15000000000000002

all.equal(n1,n2)
# [1] TRUE

all.equal(n1,n3)
# [1] TRUE

如果愿意，您可以手动指定容差，例如，

all.equal(n1, n3, tolerance = 1.5e-16)
# [1] "Mean relative difference: 1.850372e-16"

最后，正如帮助页面all.equal所说，如果您需要返回 bool，请将其包装在isTRUE(all.equal(...))or中identical。

推荐阅读