対数正規分布の等価性を正規分布の指数で見たいと思っています。私は、次のパラメータ化次に、対応するドキュメントからscipyのダウンロードパッケージを使用します。lognormalとexp(normal)の確率密度の等価性
s = sigma
scale = exp(mu)
ここthe linkです。私はこのパラメータ化を使用するとき、対数正規分布をプロットすることから正規分布と対数正規分布の対数をプロットすることから偏差を見つけることはないと思います。私はフォローしようとしました:
from scipy.stats import lognorm
from scipy.stats import norm
import matplotlib.pyplot as plt
import numpy as np
mu = 0.5
sigma = 0.2
x1 = np.linspace(norm.ppf(0.01, loc = mu, scale = sigma),
norm.ppf(0.99, loc = mu, scale = sigma), 100)
y1 = norm.pdf(x1, loc = mu, scale = sigma)
x2 = np.linspace(lognorm.ppf(0.01, s = sigma, scale = np.exp(mu)),
lognorm.ppf(0.99, s = sigma, scale = np.exp(mu)), 100)
y2 = lognorm.pdf(x2, s = sigma, scale = np.exp(mu))
fig, ax = plt.subplots(2, 1)
ax[0].plot(x1,y1, label = "normal")
ax[0].plot(np.log(x2), y2, label = "lognormal")
ax[0].grid()
ax[0].legend()
ax[1].plot(np.exp(x1),y1, label = "normal")
ax[1].plot(x2, y2, label = "lognormal")
ax[1].grid()
ax[1].legend()
fig.show()
ご覧のとおり、曲線は同じではありません。 SciPyの文書は間違っていますか?別のパラメータを使用する必要がありますか、あるいは間違いがないと私は考えていますか?