Java中的反向累积beta分布函数，具有较小的Double值

Question

我正在使用apache-commons-math进行反向累积beta分布函数，但是对于小于2E-15的值，它会给出不正确的结果。最高为2E-15，结果是正确的，并且也从R。

验证

Value 2E-15
Result in R = -7.854929
Result in Java  = -7.8529

Value 2E-16
Result in R = -8.1385
Result in Java  = -35355339.059

任何人都可以建议我任何解决方案

Answer 1

Beta分布的逆CDF范围为[0,1]。由于R和Apache Commons Math都返回负值，因此在计算这种小概率的逆CDF时显然存在一些数值误差。

您的问题可能没有数字解决方案。相反，将反向CDF的所有负输出视为零应该足够了。