Javaの逆累積ベータ分布関数が小さいDouble値を持つ

Question

逆累積ベータ分布関数にapache-commons-mathを使用していますが、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

逆CDFは、[0の範囲を持っています1 ]。 RとApache Commons Mathの両方が負の値を返すので、そのような小さな確率の逆CDFを計算する際には明らかにいくつかの数値エラーがあります。

問題の数値的解決策はありません。代わりに、逆CDFからのすべての負の出力をゼロとして扱うだけで十分です。