如何在Python中绘制悬挂根图？

受此启发question https://mathematica.stackexchange.com/questions/120174/how-to-draw-a-hanging-rootogram-in-mathematica，如何在 python 中制作相同类型的图？该图旨在以直观的方式展示您的分布与预期分布的偏差。它将直方图的条形挂在预期分布线上，因此预期值的差异是在条形底部和 x 轴之间读取的，而不是在条形顶部和预期分布曲线之间读取的。

我找不到任何内置功能。

这个想法是仅移动直方图的每个条形，使条形顶部处于预期值：

import matplotlib.pyplot as plt
import numpy as np
import matplotlib.mlab as mlab

fig, ax = plt.subplots(1, 2)
mu = 10
sig = 0.3
my_data = np.random.normal(mu, sig, 200)
x = np.linspace(9, 11, 100)

# I plot the data twice, one for the histogram only for comparison,
# and one for the rootogram.
# The trick will be to modify the histogram to make it hang to
# the expected distribution curve:

for a in ax:
    a.hist(my_data, normed=True)
    a.plot(x, mlab.normpdf(x, mu, sig))
    a.set_ylim(-0.2)
    a.set_xlim(9, 11)
    a.hlines(0, 9, 11, linestyle="--")

for rectangle in ax[1].patches:

    # expected value in the middle of the bar
    exp = mlab.normpdf(rectangle.get_x() + rectangle.get_width()/2., mu, sig)

    # difference to the expected value
    diff = exp - rectangle.get_height()
    rectangle.set_y(diff)

    ax[1].plot(rectangle.get_x() + rectangle.get_width()/2., exp, "ro")

ax[0].set_title("histogram")
ax[1].set_title("hanging rootogram")
plt.tight_layout()

这使：

HTH