ロジスティック分布をPythonで

numpyを使ってロジスティック分布を書いてみました。

import numpy as np
from  matplotlib import pyplot as plt

sigma = 1.0
mu = 0

x = np.arange(-5., 5., 0.001)
y = np.exp(-(x-mu)/sigma) / (sigma*(1+np.exp(-(x-mu)/sigma))**2)
    plt.plot(x, y)

plt.show()

こんな感じ。

累積分布も

import numpy as np
from matplotlib import pyplot as plt

sigma = 1.0
mu = 0

x = np.arange(-5., 5., 0.001)
y = 1. /(1 + (np.exp(-(x-mu)/sigma)))
plt.plot(x, y)

plt.show()

参考文献

numpy+matplotlibで正規分布 - Soleil cou coupé
ロジスティック分布 - Wikipedia