3

Julia では、正規分布の pdf または cdf を実装して、任意精度の BigFloat をサポートする人がいます。

たとえば、このコードは 0.0 を返しますが、実際には値がわずかに異なるはずです。

x = parse(BigFloat, "2.1")
x_small = float64(x)
pdf(Normal(), x) - pdf(Normal(), x_small)
4

1 に答える 1

4

直接ではありません。私の最終的な計画は、ディストリビューション タイプをパラメトリックにすることです。これにより、Float32引数も可能になりますが、それはまだしばらく先のことです。

φそれまでの間、あなたが望む結果を与えるエクスポートされていないものがあります:

Distributions.φ(x) - pdf(Normal(), x_small)
于 2015-05-26T10:49:31.083 に答える