具有較小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的所有負輸出視爲零即可。