PHPでの二項分布の計算について

Question

前提・実現したいこと

PHPの二項分布を計算できるライブラリを使い二項分布の計算を行うシステムを作成していますが、大きい数字の計算になるとNaNとなり正常に計算ができません。

発生している問題・エラーメッセージ

前提・実現したいことと同様

該当のソースコード

stats_dens_pmf_binomial(x, n, pi)

試したこと

stats_dens_pmf_binomial(x, n, pi)関数に引数を渡して計算しており基本的には動作しているのですが、大きい数字で計算すると計算結果がNaNとなります。
（例 stats_dens_pmf_binomial(2000, 8000, 0.157480314960630）等）

補足情報（FW/ツールのバージョンなど）

■使用ライブラリ
https://www.php.net/manual/ja/book.stats.php