提问人：mshaffer 提问时间：9/8/2020 更新时间：9/8/2020 访问量：716
使用 R：如何比较两个数字与与选项（digits=n）相关的浮点问题以及如何引入这些数字？[复制]
Using R: How to compare two numbers with floating-point issues related to the options(digits=n) and how the numbers were introduced? [duplicate]
问：
因此，让我们从常规 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
我在 Windoze 10 上运行 64 位。
现在让我们创建一个向量或导入一个 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
答：
1赞
heds1
9/8/2020
#1
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.equalisTRUE(all.equal(...))identical
评论
0赞
mshaffer
9/9/2020
谢谢。根据手册：所以我想有时可能会返回 NULL。isTRUE(x) is the same as { is.logical(x) && length(x) == 1 && !is.na(x) && x }all.equal
评论