Python 数据科学中的 Seaborn 绘图可视化

2023-02-27

数据 kde sns

 Python中文社区（ID:python-china）本篇文章主要研究Seaborn库——Seaborn是一个统计绘图库，建立在Matplotlib基础之上。它具有非常漂亮的默认绘图样式，并且也与PandasDataFrames配合得很好。Seaborn可以像安装任何其他Python包一

Python中文社区（ID:python-china）

本篇文章主要研究 Seaborn 库——Seaborn 是一个统计绘图库，建立在 Matplotlib 基础之上。它具有非常漂亮的默认绘图样式，并且也与 Pandas DataFrames 配合得很好。

Seaborn 可以像安装任何其他 Python 包一样使用“pip”进行安装。

pip install seaborn 
1.

Seaborn 的官方文档地址如下：

https://seaborn.pydata.org/

另一个重要的地址是官方 API，它引用了各种可用的绘图类型。

https://seaborn.pydata.org/api.html

我将尝试将 Seaborn 的功能分解为不同的类别——并从使我们能够可视化数据集分布的图开始。

让我们从导入开始并加载数据 - 我将使用“Financial Sample.xlsx”数据。excel文件可以长按扫码文末二维码后进入宽客量化俱乐部下载：

Financial Sample.xlsx

import pandas as pd  
import seaborn as sns  
#如果使用 Jupyter Notebooks，下面的行允许我们在浏览器中显示图表  
%matplotlib inline  
#在 Pandas DataFrame 中加载我们的数据  
df = pd.read_excel('Financial Sample.xlsx')   
#打印前 5 行数据以确保正确加载  
df.head() 
1.
2.
3.
4.
5.
6.
7.
8.

让我们首先看看“distplot”——这让我们可以看到一组单变量观测值的分布——单变量只是意味着一个变量。

#绘制 DataFrame "Profit" 列的分布  
sns.displot(df['Profit']) 
1.
2.

我们现在有一个感兴趣的分布图——但作为一个快速入门，风格看起来有点平淡。让我们给它一个更常见的“Seaborn”样式，试图让它看起来更好一点……

#设置我们希望用于绘图的样式  
sns.set_style("darkgrid")  
#绘制 DataFrame "Profit" 列的分布  
sns.displot(df['Profit']) 
1.
2.
3.
4.

我们已经设法用一行代码绘制 DataFrame 数据的直方图以及“KDE”线——即核密度估计图。如果我们在 plot 调用中添加“kde=False”，我们可以删除 KDE。我们还可以按如下方式更改直方图中“bins”的数量——在本例中，它们被设置为 50：

sns.displot(df['Profit'],kde=False,bins=50) 
1.

现在让我们看一个“联合图”——这允许我们组合两个 distplots 并处理双变量数据。让我们创建一个快速的联合图。为此，我们需要通过传入列名来指定我们想要绘制的 DataFrame 列，以及我们从中提取列的实际 DataFrame。这可以按如下方式完成：假设我想绘制“Profit”列与“Units Sold”列。

sns.jointplot(x='Profit',y='Units Sold',data=df) 
1.

我们现在有一个图，显示了两个变量列之间的散点图，以及它们在任一侧的相应分布图（它甚至在右上角为我们提供了皮尔逊相关系数和 p 分数。）

Jointplot 还允许我们设置一个名为“kind”的附加参数。这允许您影响主图表的表示方式。目前它是一个“散点”，因为这是默认值，但是如果我们将其更改为“十六进制”，例如，我们将得到以下图，它将图表上的点表示为密度六边形 - 即包含更多数据点的六边形显示为比包含较少点的那些更暗。

sns.jointplot(x='Profit',y='Units Sold',data=df,kind='hex') 
1.

我们可以为“kind”添加的另一个参数是“reg”，它代表回归。这看起来很像散点图，但这次将添加线性回归线。

sns.jointplot(x='Profit',y='Units Sold',data=df,kind='reg') 
1.

我们可以规定的另一种类型是“kde”，它将绘制一个二维 KDE 图，它基本上只显示数据点最常出现的位置的密度。

sns.jointplot(x='Profit',y='Units Sold',data=df,kind='kde') 
1.

让我们从jointplots继续看“pairplots”。这些使我们能够查看整个数据帧（对于数值数据）的成对关系，并且还支持分类数据点的“色调”参数。所以 pairplot 本质上是为 DataFrame 中数字列的每个可能组合创建一个联合图。我将快速创建一个新的 DataFrame，它删除“Month Number”和“Year”列，因为这些并不是我们连续数字数据的一部分，例如“利润”和“COGS”（销售成本）。我还将删除其他几列以缩小我们的 DataFrame，这样我们的输出图就不会过于拥挤。

#删除不需要的列  
new_df = df.drop(['Month Number','Year','Manufacturing Price','Sale Price'],axis=1)  
sns.pairplot(new_df) 
1.
2.
3.

请注意，我们基本上对每对列都有一个配对图，并且在对角线上我们有一个分布的直方图，因为将数据与自身进行联合图是没有意义的。这是快速可视化数据的好方法。我们还可以添加一个“色调”——这是我们指定一个用于分割数据的分类变量的地方。让我们添加“Segment”列作为我们的“色调”。

sns.pairplot(new_df,hue='Segment') 
1.

现在数据点根据分类数据着色——颜色图例显示在图的右侧边缘。我们还可以通过设置“调色板”参数来更改绘图使用的调色板。以下是使用“岩浆”配色方案的示例。所有可用的方案都可以在 Matplotlib 站点上找到。

sns.pairplot(new_df,hue='Segment',palette='magma') 
1.

我们将看到的下一个图是一个“rugplot”——这将帮助我们构建和解释我们之前创建的“kde”图是什么——无论是在我们的 distplot 中还是当我们传递“kind=kde”作为我们的参数时。

sns.rugplot(df['Profit']) 
1.

如上所示，对于 rugplot，我们将要绘制的列作为参数传递 - rugplot 的作用是为分布中的每个点绘制一个破折号。所以 rugplot 和 distplot 之间的区别在于 distplot 涉及“bins”的概念，并将把每个 bin 中的所有数据点相加，并绘制这个数字，而 rugplot 只是在每个数据点绘制一个标记。

所以现在让我们将 rugplot 转换为 KDE 图。KDE 代表“核密度估计”。下图是解释如何将 rugplots 构建到 KDE 图中的。

如果我们愿意的话，我们可以从一组数据和 rugplot 中构建我们自己的 KDE 图，看看它是否与使用内置的“kdeplot”直接创建的 KDE 图相匹配.

#设置一组 30 个取自正态分布的数据点  
x = np.random.normal(0, 1, size=30)  
#设置 KDE 点的带宽  
bandwidth = 1.06* x.std() * x.size ** (-1/ 5.)  
#设置 y 轴的限制  
support = np.linspace(-4, 4, 200)  
#遍历数据点并为每个点创建内核，然后绘制内核  
kernels = []  
for x_i in x:  
    kernel = stats.norm(x_i, bandwidth).pdf(support)  
    kernels.append(kernel)  
    plt.plot(support, kernel, color="r")  
sns.rugplot(x, color=".2", linewidth=3) 
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.

#使用复合梯形规则沿给定轴积分并创建 KDE 图  
from scipy.integrate import trapz  
density = np.sum(kernels, axis=0)  
density /= trapz(density, support)  
plt.plot(support, density) 
1.
2.
3.
4.
5.