더 작은 Double 값을 갖는 Java의 역 누적 베타 분포 함수

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 

사람은 나에게 모든 솔루션을 제안 해주십시오 수

베타 분포에 대한 역 CDF는 [0, 1의 범위가

Answer 1

]. R과 Apache Commons Math가 모두 음수 값을 반환하기 때문에 이러한 작은 확률에 대해 역 CDF를 계산할 때 몇 가지 수치 오류가 분명히 있습니다.

문제에 대한 수치 해결책이 없을 수 있습니다. 대신 역 CDF의 모든 음수 출력을 0으로 처리하면 충분합니다.