logsample = stats.norm.rvs(loc=10, scale=3, size=1000)
sample = np.exp(logsample)
shape, loc, scale = stats.lognorm.fit(sample, floc=0)
shape, loc, scale
np.log(scale), shape
scipy, lognormal distribution - parameters