numpy - 高い値を中心とした対数正規乱数

Question

numpy/scipy を使用して対数正規分布から乱数を作成しようとしています。

平均は 2000、シグマは 800 として与えられます。

numpy.random.lognormal(mean=2000, sigma=800, size=10000) を使用してランダムな値を作成すると、非常に高い数値または inf の数値しか得られません。

これを回避する方法はありますか?

Answer 1

注意:meanおよびsigma引数は、ディストリビューションのログのディストリビューションに対応しますlognormal。分布の実際の算術平均は であり、 と の場合、標準の倍精度浮動小数点exp(mean + sigma**2/2)で評価されます。infmean=2000sigma=800

詳細については 、 http://docs.scipy.org/doc/numpy/reference/generated/numpy.random.lognormal.html#numpy.random.lognormal および http://en.wikipedia.org/wiki/Log-normal_distribution を参照してください。詳細。