如何解决使用R:如何将两个数字与与选项digits = 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
在哪里我都正确?
解决方法
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
的帮助页面上所述,如果您需要返回布尔值,请将其包装在isTRUE(all.equal(...))
或identical
中。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。