haskell / haskell/statistics

Problems with beta distrbution

Open
#161 0 comments 0 reactions 0 assignees View on GitHub
Bug
Dominant language
Haskell
Stars
309
Forks
76
PR merge metrics
No merged PRs in 30d

Description

Test for quantile inverse fails spectacularly:
```
Quantile is CDF inverse: FAIL
*** Failed! Falsified (after 99 tests):
improperBetaDistr 1.9354985700057857e-2 9.298452506189731
Double01 0.7873411995889938
Quantile = 9.435096862311606e-5
Probability = 0.7873411995889938
Probability' = 0.8812454856896215
Rel. error = 0.10655860100905087
Abs. error = 9.390428610062762e-2
Expected err. = 16.346348384470552
Distance = 845814615782558
Err/est = 5.1743337159392945e13 <---
Use --quickcheck-replay=315907 to reproduce.
```
This is really bug in incompleteBeta and/or invIncompleteBeta. Actual roundtrip error is so bad: 1e13 larger that expected that there's no point in clearing

Probably related: bos/math-functions#32

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.