Python：seaborn的散点图矩阵（Pairs Plots）可视化数据

矩阵式决策: 利用矩阵图评估项目风险与收益 #生活技巧# #团队建设技巧# #团队决策方法#

散点图矩阵让我们看到了两个变量之间的关系，分析趋势。

在本文中，我们将通过使用seaborn可视化库在Python中进行对图的绘制和运行。我们将看到如何创建默认配对图以快速检查我们的数据，以及如何自定义可视化以获取更深入的洞察力。该项目的代码在GitHub上以Jupyter Notebook的形式提供。在这个项目中，我们将探索一个真实世界的数据集，由GapMinder收集的国家级社会经济数据组成。

Seaborn的散点图矩阵（Pairs Plots）

在开始之前，我们需要知道我们有什么数据。我们可以将社会经济数据用熊猫（Pandas）数据框加载并查看列：

每行数据代表一个国家在一年内的结果，列中包含变量（这种格式的数据称为整洁数据）。有2个分类专栏（国家和大陆）和4个数字专栏。这些专栏包括：life_exp是几年出生时的预期寿命，pop是人口，gdp_per_cap是以国际美元为单位的人均国内生产总值。

虽然后面我们将使用分类变量进行着色，但seaborn中的默认对图仅绘制了数字列。创建默认的散点图矩阵很简单：我们加载到seaborn库并调用pairplot函数，将它传递给我们的数据框：

import seaborn as sns

sns.pairplot(df)

'

我仍然惊讶于一行简单的代码就可以完成我们整个需求！散点图矩阵建立在两个基本图形上，直方图和散点图。对角线上的直方图允许我们看到单个变量的分布，而上下三角形上的散点图显示了两个变量之间的关系。例如，第二行中最左边的图表显示life_exp与年份的散点图。

默认的散点图矩阵图经常给我们提供有价值的见解。我们看到人均预期寿命和人均GDP是正相关的，这表明高收入国家的人们倾向于更长寿（尽管这当然不能证明导致其他人也是如此）。它也似乎看起来是全世界的预期寿命随着时间的推移而上升。为了在未来的图中更好地显示这些变量，我们可以通过取这些值的对数来转换这些列：

df['log_pop'] = np.log10(df['pop'])

df['log_gdp_per_cap'] = np.log10(df['gdp_per_cap'])

df = df.drop(columns = ['pop', 'gdp_per_cap'])

虽然这种制图本身可以用于分析，但我们可以发现，通过对诸如大陆这样的分类变量进行数字着色，使其更有价值。这在seaborn中非常简单！我们所需要做的就是在hue中使用sns.pairplot函数调用使用关键字：

sns.pairplot(df, hue = 'continent')

'

现在我们看到大洋洲和欧洲的人均预期寿命最高，亚洲人口最多。请注意，我们对人口和gdp的日志转换使这些变量正态分布，从而更全面地表示值。

上图更具信息性，但仍然存在一些问题：找不到叠加的直方图，就像在对角线上那样，它非常易于理解。显示来自多个类别的单变量分布的更好方法是密度图。我们可以在函数调用中交换柱状图的密度图。当我们处理它时，我们会将一些关键字传递给散点图，以更改点的透明度，大小和边缘颜色。

sns.pairplot(df, hue = 'continent', diag_kind = 'kde',

plot_kws = {'alpha': 0.6, 's': 80, 'edgecolor': 'k'},

size = 4)

'

对角线上的密度图比堆积条更容易比较各大洲之间的分布。改变散点图的透明度可以提高可读性，因为这些数字有相当多的重叠（称为重叠绘图）。

作为pairplot默认的最后一个例子，让我们通过绘制2000年后的年份来减少数据混乱。我们仍然会按照大陆分布着色，但现在我们不会绘制年份列。为了限制绘制的列，我们将一个列表传递vars给函数。为了说明情节，我们还可以添加标题。

sns.pairplot(df[df['year'] >= 2000],

vars = ['life_exp', 'log_pop', 'log_gdp_per_cap'],

hue = 'continent', diag_kind = 'kde',

plot_kws = {'alpha': 0.6, 's': 80, 'edgecolor': 'k'},

size = 4);

plt.suptitle('Pair Plot of Socioeconomic Data for 2000-2007',

size = 28);

'

这开始看起来很不错！如果我们要进行建模，我们可以使用这些图中的信息来帮助我们进行选择。例如，我们知道log_gdp_per_cap与life_exp正相关，所以我们可以创建一个线性模型来量化这种关系。对于这篇文章，我们将坚持绘图，如果我们想要更多地探索我们的数据，我们可以使用PairGrid类自定义散点图矩阵。

使用PairGrid进行自定义

与sns.pairplot函数相反，sns.PairGrid是一个类，它意味着它不会自动填充我们的网格plot。相反，我们创建一个类实例，然后将特定函数映射到网格的不同部分。要用我们的数据创建一个PairGrid实例，我们使用下面的代码，这也限制了我们将显示的变量：

grid = sns.PairGrid(data= df_log[df_log['year'] == 2007],

vars = ['life_exp', 'log_pop',

'log_gdp_per_cap'], size = 4)

如果我们要显示它，我们会得到一个空白图，因为我们没有将任何函数映射到网格部分。有三个网格部分填写PairGrid：上三角形、下三角形和对角线。要将网格映射到这些部分，我们使用grid.map 部分中的方法。例如，要将散点图映射到我们使用的上三角形：

grid = grid.map_upper(plt.scatter, color = 'darkred')

该map_upper方法接受任何两个变量数组（如plt.scatter）和关联的关键字（如color）的函数。该map_lower方法完全相同，但填充网格的下三角形。因为它需要在接受单个阵列（记住对角线仅示出了一个变量）的函数略有不同。一个例子是plt.hist我们用来填写下面的对角线部分：

grid = grid.map_diag(plt.hist, bins = 10, color = 'darkred',

edgecolor = 'k')

grid = grid.map_lower(sns.kdeplot, cmap = 'Reds')

在这种情况下，我们在下三角形中使用2-D（密度图）的核密度估计值。放在一起，这段代码给了我们下面的图：

使用PairGrid类的真正好处在于我们想要创建自定义函数来将不同的信息映射到图上。例如，我可能想要将两个变量之间的Pearson相关系数添加到散点图中。为此，我会编写一个函数，它接受两个数组、计算统计量，然后在图上绘制它。下面的代码显示了这是如何完成的（归功于这个Stack Overflow答案）：

def corr(x, y, **kwargs):

coef = np.corrcoef(x, y)[0][1]

label = r'$\rho$ = ' + str(round(coef, 2))

ax = plt.gca()

ax.annotate(label, xy = (0.2, 0.95), size = 20, xycoords = ax.transAxes)

grid = sns.PairGrid(data= df[df['year'] == 2007],

vars = ['life_exp', 'log_pop', 'log_gdp_per_cap'], size = 4)

grid = grid.map_upper(plt.scatter, color = 'darkred')

grid = grid.map_upper(corr)

grid = grid.map_lower(sns.kdeplot, cmap = 'Reds')

grid = grid.map_diag(plt.hist, bins = 10, edgecolor = 'k', color = 'darkred');

我们的新函数映射到上三角形，因为我们需要两个数组来计算相关系数（还要注意我们可以将多个函数映射到网格部分）。这会产生以下图形：

相关系数现在出现在散点图上方。这是一个相对直接的例子，但我们可以使用PairGrid将我们想要的任何函数映射到图上。我们可以根据需要添加尽可能多的信息，只要我们能够弄清楚如何编写函数！作为最后一个例子，这里是一个显示对角线而不是网格的汇总统计图。

它显示了我们只做图标的总体思路，除了使用库中的任何现有功能（例如matplotlib将数据映射到图上）之外，我们还可以编写自己的函数来显示自定义信息。

结论

散点图矩阵是快速探索数据集中的分布和关系的强大工具。Seaborn提供了一个简单的默认方法，可以通过Pair Grid类来定制和扩展散点图矩阵。在一个数据分析项目中，价值的主要部分往往不在于浮华的机器学习，而在于直观的数据可视化。散点图举证为我们提供了全面的数据分析，是数据分析项目的一个很好的起点。

文章原标题《visualizing-data-with-pair-plots-in-python》，

转自：https://yq.aliyun.com/articles/581263

网址：Python：seaborn的散点图矩阵（Pairs Plots）可视化数据 https://www.yuejiaxmz.com/news/view/259522

相关内容

正交矩阵； 实对称矩阵； 为什么实对称矩阵一定可以对角化； AB=0 r(A)+r(B)＜=n 证明； 初等矩阵； 初等矩阵的逆矩阵； 矩阵的左除右除；
短视频矩阵一键批量发布，引爆流量新高度！
线性代数中，从矩阵AB=E可以推出AB=BA吗
数据分享|R语言、SPSS基于主成分PCA的中国城镇居民消费结构研究可视化分析
谁给点提示==囧？求证凸函数的Hessian矩阵H是半正定的. 爱问知识人
相对熵与KL散度：推荐系统与个性化优化1.背景介绍 随着互联网的普及和数据的崛起，人工智能技术已经成为了我们生活中不可或
蓝桥杯 矩阵键盘的使用
Python数据分析：对饮食与健康数据的分析与可视化
《python数据分析与挖掘》
Python数据可视化

随便看看