R中的概率之和是一个浮点问题

问题描述

我正面临一个与概率之和有关的浮点问题。由于有些微小差异,很难证明概率之和为1。 MWE如下。

p1 = 0.99999999

p2 = 0.00000000003

p1 + p2 == 1

[1]错误

概率之和约为1。与1的差是1-p1-p2 = 9.97e-09,这很小。我需要在许多函数中应用概率条件的总和。但是由于浮点,执行被暂停。

有人可以指导我吗?

谢谢

致谢

解决方法

比较浮点数时的一般原则不是测试相等性,而是测试差异是否小于某个ε,即您将接受的最大差异为“相等”。就您而言,您需要测试abs(p1+p2-1) < epsilon。搜索浮点比较,网络上有很多文档和讨论,例如How should I do floating point comparison?