4

我正在尝试使用 Matplotlib 的 2d 散点图函数绘制一些数据，同时在 x 和 y 轴上生成投影直方图。我找到的示例来自 matplotlib 图像库（pylab_examples 示例代码： scatter_hist.py）。

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import NullFormatter

# the random data
x = np.random.randn(1000)
y = np.random.randn(1000)

nullfmt   = NullFormatter()         # no labels

# definitions for the axes
left, width = 0.1, 0.65
bottom, height = 0.1, 0.65
bottom_h = left_h = left+width+0.02

rect_scatter = [left, bottom, width, height]
rect_histx = [left, bottom_h, width, 0.2]
rect_histy = [left_h, bottom, 0.2, height]

# start with a rectangular Figure
plt.figure(1, figsize=(8,8))

axScatter = plt.axes(rect_scatter)
axHistx = plt.axes(rect_histx)
axHisty = plt.axes(rect_histy)

# no labels
axHistx.xaxis.set_major_formatter(nullfmt)
axHisty.yaxis.set_major_formatter(nullfmt)

# the scatter plot:
axScatter.scatter(x, y)

# now determine nice limits by hand:
binwidth = 0.25
xymax = np.max( [np.max(np.fabs(x)), np.max(np.fabs(y))] )
lim = ( int(xymax/binwidth) + 1) * binwidth

axScatter.set_xlim( (-lim, lim) )
axScatter.set_ylim( (-lim, lim) )

bins = np.arange(-lim, lim + binwidth, binwidth)

axHistx.hist(x, bins=bins)
axHisty.hist(y, bins=bins, orientation='horizontal')

axHistx.set_xlim( axScatter.get_xlim() )
axHisty.set_ylim( axScatter.get_ylim() )

plt.show()

唯一的问题是该示例不起作用。我收到以下错误：

~$ python ~/Desktop/scatter_and_hist.py 
Traceback (most recent call last):
  File "/Users/username/Desktop/scatter_and_hist.py", line 45, in <module>
    axHisty.hist(y, bins=bins, orientation='horizontal')
  File "//anaconda/lib/python2.7/site-packages/matplotlib/axes.py", line 8180, in hist
    color=c, bottom=bottom)
TypeError: barh() got multiple values for keyword argument 'bottom'

我已经浏览了代码并隔离了问题。导致问题的是第 45 行 (axHisty.hist(y, bins=bins,orientation='horizo​​ntal'))。在他们的图片库中看到你想要的情节真是令人沮丧，但是这个例子不起作用。第二组眼睛将不胜感激！

4

1 回答 1

5

您在 v1.2.1 ( https://github.com/matplotlib/matplotlib/pull/1985 )中遇到了一个错误。您可以升级您的 matplotlib，使用错误修复对您的版本进行猴子补丁，或者使用您自己的正确参数顺序进行np.histogram调用barh。

作为旁注，此问题所需的唯一代码是：

x = np.random.rand(100)
plt.hist(x, orientation='horizontal')
plt.show()

您发布的其他所有内容都是噪音。

于 2013-08-26T21:14:55.633 回答