matplotlib中的多步直方图

时间:2014-11-01 18:46:46

标签: python matplotlib plot histogram

亲爱的python / matplotlib社区,

我在matplotlib中遇到问题:我似乎无法使用以下内容在同一绘图空间中绘制多个重叠的直方图:

binsize = 0.05

min_x_data_sey, max_x_data_sey    = np.min(logOII_OIII_sey), np.max(logOII_OIII_sey)
num_x_bins_sey                    = np.floor((max_x_data_sey - min_x_data_sey) / binsize)

min_x_data_comp, max_x_data_comp  = np.min(logOII_OIII_comp), np.max(logOII_OIII_comp)
num_x_bins_comp                   = np.floor((max_x_data_comp - min_x_data_comp) / binsize)

min_x_data_sf, max_x_data_sf      = np.min(logOII_OIII_sf), np.max(logOII_OIII_sf)
num_x_bins_sf                     = np.floor((max_x_data_sf - min_x_data_sf) / binsize)

axScatter_farright = fig.add_subplot(gs_right[0,0])

axScatter_farright.tick_params(axis='both', which='major', labelsize=10)
axScatter_farright.tick_params(axis='both', which='minor', labelsize=10)
axScatter_farright.set_ylabel(r'$\mathrm{N}$', fontsize='medium')
axScatter_farright.set_xlim(-1.5, 1.0)
axScatter_farright.set_xlabel(r'$\mathrm{log([OII]/[OIII])}$', fontsize='medium')

axScatter_farright.hist(logOII_OIII_sey, num_x_bins_sey, ec='0.3', fc='none', histtype='step')
axScatter_farright.hist(logOII_OIII_comp, num_x_bins_comp, ec='0.3', fc='none', histtype='step')
axScatter_farright.hist(logOII_OIII_sf, num_x_bins_sf, ec='0.3', fc='none', histtype='step')

似乎轴类无法处理多个直方图?如果和/或我哪里出错了,请纠正我。

我的总体情节是1行,3列绘图空间。我想使用网格规格给图块一个很好的布局。

这是我的情节到目前为止的样子:

enter image description here

这就是我希望图形的直方图部分在步骤类型直方图叠加(带图例)方面看起来像:

enter image description here

我将数据集作为从csv文件生成的三个不同的元组类型数组。即使用x, y = np.genfromtext(datafile.csv)

如果有人能够解释如何做到这一点,我将非常感激。

1 个答案:

答案 0 :(得分:3)

你正在做的事应该完美无缺。是否有可能只有一个分布在-1.5到1的x范围内,你之前设置了几行? (即尝试删除手动set_xlim语句,看看其他发行版是否出现。)

作为一个快速,独立的例子来证明事情应该有效:

import numpy as np
import matplotlib.pyplot as plt

num = 1000
d1 = np.random.normal(-1, 1, num)
d2 = np.random.normal(1, 1, num)
d3 = np.random.normal(0, 3, num)

fig, ax = plt.subplots()
ax.hist(d1, 50, ec='red', fc='none', lw=1.5, histtype='step', label='Dist A')
ax.hist(d2, 50, ec='green', fc='none', lw=1.5, histtype='step', label='Dist B')
ax.hist(d3, 100, ec='blue', fc='none', lw=1.5, histtype='step', label='Dist C')
ax.legend(loc='upper left')
plt.show()

enter image description here

(如果您希望图例显示行而不是框,则您需要使用代理艺术家。如果您愿意,我可以添加一个示例。这不在此范围之内但问题是。)