带误差线的 Matplotlib 直方图

Question

我使用pyplot.hist() 函数创建了带有matplotlib 的直方图。我想在条形图中添加一个 bin 高度 (sqrt(binheight)) 的中毒误差平方根。我该怎么做？

.hist() 的返回元组包括return[2] -> 1 个 Patch 对象的列表。我只能发现可以将错误添加到通过pyplot.bar() 创建的条形图上。

Answer 1

确实需要使用 bar。您可以使用hist 的输出并将其绘制为条形：

import numpy as np
import pylab as plt

data       = np.array(np.random.rand(1000))
y,binEdges = np.histogram(data,bins=10)
bincenters = 0.5*(binEdges[1:]+binEdges[:-1])
menStd     = np.sqrt(y)
width      = 0.05
plt.bar(bincenters, y, width=width, color='r', yerr=menStd)
plt.show()

Answer 2

替代方案

您还可以使用pyplot.errorbar() 和drawstyle 关键字参数的组合。下面的代码使用阶梯线图创建直方图。每个 bin 的中心都有一个标记，每个 bin 都有必要的泊松误差条。

import numpy
import pyplot

x = numpy.random.rand(1000)
y, bin_edges = numpy.histogram(x, bins=10)
bin_centers = 0.5*(bin_edges[1:] + bin_edges[:-1])

pyplot.errorbar(
    bin_centers,
    y,
    yerr = y**0.5,
    marker = '.',
    drawstyle = 'steps-mid-'
)
pyplot.show()

我的个人意见

在同一张图上绘制多个直方图的结果时，线图更容易区分。此外，使用yscale='log' 绘图时它们看起来更好。