具有较小Double值的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的所有负输出视为零即可。