我试图拟合对数正态曲线,但我不知道如何将y(和x轴)更改为对数刻度?
我尝试使用logspace而不是linspace,但是遇到了一个错误OverflowError: cannot convert float infinity to integer
我试过了:
^{pr2}$
但这似乎也不管用。我附上了我的图表来告诉你我的意思。我真的看不到400左右发生了什么,所以我想用对数标度y(也可以用对数标度x来比较)from collections import Counter
import matplotlib.pyplot as plt
from matplotlib.backends.backend_pdf import PdfPages
from scipy.stats import lognorm
import numpy as np
data = list(pre_data)
params = (0.40951093774076597, 5.1802137214177684, 60.158303995566413)
shape, loc, scale = params[0], params[1], params[2]
print params
prob = 1-lognorm.cdf(388,shape,loc=params[1], scale=params[2])
print prob * 2994
count, bins, ignored = plt.hist(data,100,normed=True)
mu = np.mean(np.log(data))
sigma = np.std(np.log(data))
x = np.linspace(min(bins),max(bins),10000)
pdf = (np.exp(-(np.log(x)-mu)**2 / (2 * sigma**2)) / (x * sigma * np.sqrt(2*np.pi)))
plt.plot(x,pdf,color='r',linewidth= 2)