我用对数正态pdf拟合一些二进制数据,但我的曲线与数据不太匹配,见下图。我的代码是:
import numpy as np
import matplotlib.pyplot as plt
from scipy.stats import lognorm
data = genfromtxt('data.txt')
data = np.sort(data)
# plot histogram in log space
ax.hist(data, bins=np.logspace(0,5,200),normed=1)
ax.set_xscale("log")
shape,loc,scale = lognorm.fit(data)
print shape, loc, scale
pdf = sp.stats.lognorm.pdf(data, shape, loc, scale)
ax.plot(data,pdf)
plt.show()
它看起来是这样的:
我是否需要对形状、位置和比例进行合理的猜测?在
谢谢!在
您试图拟合的数据看起来不像对数正态分布。对数正态分布,当以对数x标度绘制时,应该看起来像正态分布。你展示的情节不是这样的。当分布不适合数据时,会得到奇怪的参数。在
您需要先了解数据是如何分布的(严格地说,这与SO无关),然后才能尝试适合某些内容。在
这是我们使用从对数正态分布随机抽取的数据时得到的结果:
相关问题 更多 >
编程相关推荐