使用Seaborn FacetGrid从数据框中绘制错误条

时间:2014-07-22 02:48:56

标签: python matplotlib plot pandas seaborn

我想在Seaborn FacetGrid上的pandas数据框中的列中绘制误差条

import matplotlib.pyplot as plt
import pandas as pd
import seaborn as sns
df = pd.DataFrame({'A' : ['foo', 'bar', 'foo', 'bar']*2,
                   'B' : ['one', 'one', 'two', 'three',
                         'two', 'two', 'one', 'three'],
                  'C' : np.random.randn(8),
                  'D' : np.random.randn(8)})
df

示例数据框

    A       B        C           D
0   foo     one      0.445827   -0.311863
1   bar     one      0.862154   -0.229065
2   foo     two      0.290981   -0.835301
3   bar     three    0.995732    0.356807
4   foo     two      0.029311    0.631812
5   bar     two      0.023164   -0.468248
6   foo     one     -1.568248    2.508461
7   bar     three   -0.407807    0.319404

此代码适用于固定大小的错误栏:

g = sns.FacetGrid(df, col="A", hue="B", size =5)
g.map(plt.errorbar, "C", "D",yerr=0.5, fmt='o');

enter image description here

但我无法使用数据框中的值

来使用它
df['E'] = abs(df['D']*0.5)
g = sns.FacetGrid(df, col="A", hue="B", size =5)
g.map(plt.errorbar, "C", "D", yerr=df['E']);

g = sns.FacetGrid(df, col="A", hue="B", size =5)
g.map(plt.errorbar, "C", "D", yerr='E');

都会产生错误的痕迹

编辑:

经过大量的matplotlib doc阅读,以及各种stackoverflow的答案, 这是一个纯matplotlib解决方案

#define a color palette index based on column 'B'
df['cind'] = pd.Categorical(df['B']).labels

#how many categories in column 'A'
cats = df['A'].unique()
cats.sort()

#get the seaborn colour palette and convert to array
cp = sns.color_palette()
cpa = np.array(cp)

#draw a subplot for each category in column "A"
fig, axs = plt.subplots(nrows=1, ncols=len(cats), sharey=True)
for i,ax in enumerate(axs):
    df_sub = df[df['A'] == cats[i]]
    col = cpa[df_sub['cind']]
    ax.scatter(df_sub['C'], df_sub['D'], c=col)
    eb = ax.errorbar(df_sub['C'], df_sub['D'], yerr=df_sub['E'], fmt=None)
    a, (b, c), (d,) = eb.lines
    d.set_color(col)

除标签外,轴限制为OK。它为“A”列中的每个类别绘制了一个单独的子图,由“B”列中的类别着色。 (注意随机数据与上面的不同)

如果有人有任何想法,我仍然喜欢大熊猫/海豹的解决方案吗?

enter image description here

2 个答案:

答案 0 :(得分:9)

使用FacetGrid.map时,任何引用data DataFrame的内容都必须作为位置参数传递。这将适用于您的情况,因为yerrplt.errorbar的第三个位置参数,但为了证明我将使用提示数据集:

from scipy import stats
tips_all = sns.load_dataset("tips")
tips_grouped = tips_all.groupby(["smoker", "size"])
tips = tips_grouped.mean()
tips["CI"] = tips_grouped.total_bill.apply(stats.sem) * 1.96
tips.reset_index(inplace=True)

然后我可以使用FacetGriderrorbar

进行策划
g = sns.FacetGrid(tips, col="smoker", size=5)
g.map(plt.errorbar, "size", "total_bill", "CI", marker="o")

enter image description here

但是,请记住,有一个seaborn绘图功能,用于从完整数据集转到带有错误栏的图(使用自举),因此对于许多应用程序而言,这可能不是必需的。例如,您可以使用factorplot

sns.factorplot("size", "total_bill", col="smoker",
               data=tips_all, kind="point")

enter image description here

lmplot

sns.lmplot("size", "total_bill", col="smoker",
           data=tips_all, fit_reg=False, x_estimator=np.mean)

enter image description here

答案 1 :(得分:0)

您没有显示df['E']实际上是什么,以及它是否是与df['C']df['D']相同长度的列表。

yerr关键字参数(kwarg)采用单个值,该值将应用于数据帧中键C和D的列表中的每个元素,或者它需要一个与这些值相同长度的值列表列表。

因此,C,D和E必须都与相同长度的列表相关联,或者C和D必须是相同长度的列表,并且E必须与单个float或{{1 }}。如果单个intfloat位于列表中,则必须将其解压缩,例如int

包含df['E'][0]的示例matplotlib代码: http://matplotlib.org/1.2.1/examples/pylab_examples/errorbar_demo.html

描述yerr的条形图API文档: http://matplotlib.org/api/pyplot_api.html#matplotlib.pyplot.bar