2
我想要適合對數正態曲線,我不知道如何將y(和x軸)更改爲對數刻度? 我嘗試使用LOGSPACE代替linspace,但我得到一個錯誤Python中的日誌規模mathplotlib?
OverflowError: cannot convert float infinity to integer
我想:
x = np.linspace(np.log(min(bins)),np.log(max(bins)),10000)
但似乎沒有任何工作。 。 。我附上我的圖表來向你展示我的意思。我真的不能看到發生了什麼約有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)
也請選擇此答案公認。 – sebix 2014-12-09 20:31:58