15

我正在创建直方图

pylab.hist(data,weights,histtype='step',normed=False,bins=150,cumulative=True)

得到(还有其他情节,现在无关紧要)紫线

直方图

为什么直方图最后又归零了?累积函数通常应该是非递减的。有没有办法解决这个问题,无论是错误还是功能?

编辑:解决方案(黑客):

# histtype=step returns a single patch, open polygon
n,bins,patches=pylab.hist(data,weights,histtype='step',cumulative=True)
# just delete the last point
patches[0].set_xy(patches[0].get_xy()[:-1])
4

2 回答 2

1

如果您不喜欢 OP 不错的简单解决方案,这里有一个过于复杂的解决方案,我们手动构建绘图。如果您只能访问直方图计数并且不能使用 matplotlib 的 hist 函数,那么它可能会很有用。

import numpy as np
import matplotlib.pyplot as plt

data = np.random.randn(5000)
counts, bins = np.histogram(data, bins=20)
cdf = np.cumsum(counts)/np.sum(counts)

plt.plot(
    np.vstack((bins, np.roll(bins, -1))).T.flatten()[:-2],
    np.vstack((cdf, cdf)).T.flatten()
)
plt.show()

输出

于 2017-04-28T12:51:34.817 回答
0

这是默认行为。把它想象成柱状图的轮廓作为条形图。至于快速解决方法,我不知道。一种解决方案是自己计算直方图:python histogram one-liner

于 2012-05-21T18:10:11.443 回答