您可以将四分位距绘制为 seaborn 线图上的误差带吗？

Question

我正在使用 seaborn lineplot ( https://seaborn.pydata.org/generated/seaborn.lineplot.html )绘制时间序列数据，并绘制中位数而不是平均值。示例代码：

import seaborn as sns; sns.set()
import matplotlib.pyplot as plt

fmri = sns.load_dataset("fmri")
ax = sns.lineplot(x="timepoint", y="signal", estimator = np.median, data=fmri)

我希望误差带显示四分位距而不是置信区间。我知道我可以ci = "sd"用于标准偏差，但是有没有一种简单的方法来添加 IQR？我想不明白。

谢谢！

Answer 1

我不知道这是否可以单独使用 seaborn 来完成，但这是使用 matplotlib 做到这一点的一种方法，保持 seaborn 风格。该describe()方法方便地提供了 DataFrame 的汇总统计信息，其中包括四分位数，我们可以用它来绘制具有四分位间距的中位数。

import seaborn as sns; sns.set()
import matplotlib.pyplot as plt

fmri = sns.load_dataset("fmri")
fmri_stats = fmri.groupby(['timepoint']).describe()

x = fmri_stats.index
medians = fmri_stats[('signal', '50%')]
medians.name = 'signal'
quartiles1 = fmri_stats[('signal', '25%')]
quartiles3 = fmri_stats[('signal', '75%')]

ax = sns.lineplot(x, medians) 
ax.fill_between(x, quartiles1, quartiles3, alpha=0.3); 

Answer 2

您可以像之前一样计算 lineplot 中的中位数，将 ci 设置为 none 并使用 ax.fill_between()

import numpy as np
import seaborn as sns; sns.set()
import matplotlib.pyplot as plt

fmri = sns.load_dataset("fmri")
ax = sns.lineplot(x="timepoint", y="signal", estimator = np.median, 
                  data=fmri,ci=None)

bounds = fmri.groupby('timepoint')['signal'].quantile((0.25,0.75)).unstack()
ax.fill_between(x=bounds.index,y1=bounds.iloc[:,0],y2=bounds.iloc[:,1],alpha=0.1)