haskell / haskell/statistics

Problems with Student-T distribution

Abierto
#160 1 comentario 0 reacciones 0 asignados Ver en GitHub
Bug
Lenguaje dominante
Haskell
Estrellas
309
Forks
76
Métricas de merge de PR
Sin PR fusionados en 30 d

Descripción

So far following problems with Student-T distribution were identified:

1. `quantile` lose precision for p≈0.5 since x becomes close to 1, while p is near zero
```haskell
quantile :: StudentT -> Double -> Double
quantile (StudentT ndf) p
| p >= 0 && p <= 1 =
let x = invIncompleteBeta (0.5 * ndf) 0.5 (2 * min p (1 - p))
in case sqrt $ ndf * (1 - x) / x of
r | p < 0.5 -> -r
| otherwise -> r
```

2. Both precision and performance suffer greatly for large degrees of freedom. Reason is likely incomplete beta and its inverse perform poorly for large parameters. However distribution becomes close to normal so other approximations could work!

Guía de contribución

No hay ninguna guía de contribución indexada para este repositorio

Evaluación

Este issue todavía no se ha evaluado.

Recibe los nuevos issues en tu correo

Un resumen breve de issues de GitHub para principiantes.