Юлия: BigFloat Нормальное распределение

В Julia кто-нибудь реализовал нормальные дистрибутивы pdf или cdf для поддержки BigFloats произвольной точности.

Например, этот код возвращает 0.0, когда на самом деле значения должны немного отличаться.

x = parse(BigFloat, "2.1")
x_small = float64(x)
pdf(Normal(), x) - pdf(Normal(), x_small)
Источник

1 ответ

Решение

Не напрямую. Мой конечный план состоит в том, чтобы сделать типы распределения параметрическими, что также позволило бы Float32 аргументы, но это еще далеко.

А пока есть неэкспортируемый φ который дает результат, который вы хотели:

Distributions.φ(x) - pdf(Normal(), x_small)
Источник
Другие вопросы по тегам