Java中具有较小双精度值的反向累积beta分布函数

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的所有负输出视为零就足够了。