基本
#-*- coding:utf-8 -*-
from __future__ import division
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy import stats, integrate
import seaborn as sns
1、快速查看分布
distplot
x = np.random.randint(0,10,100)
sns.distplot(x, color="b", bins=10, rug=True)
sns.distplot(x, hist=False)
03_01.png
2、畫單變量圖
x = np.random.normal(0, 1, size=30)
bandwidth = 1.06 * x.std() * x.size ** (-1 / 5.) # 帶寬計(jì)算公式
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="#C55458")
sns.rugplot(x, color=".2", linewidth=3)
density = np.sum(kernels, axis=0)
density /= integrate.trapz(density, support)
plt.plot(support, density)
03_02.png
sns.kdeplot(x, shade=True)
03_03.png
sns.kdeplot(x)
sns.kdeplot(x, bw=.2, label="bw: 0.2")
sns.kdeplot(x, bw=2, label="bw: 2")
plt.legend()
03_04.png
x = np.random.gamma(6, size=200)
sns.distplot(x)
sns.distplot(x, label='kde', hist=False)
sns.distplot(x, kde=False, fit=stats.gamma, label='fit-gamma', hist=False)
plt.legend()
03_05.png
3咧七、畫二元變量分布圖
mean, cov = [0, 1], [[1, .5], [.5, 1]]
data = np.random.multivariate_normal(mean, cov, 200)
df = pd.DataFrame(data, columns=["x", "y"])
（1）散點(diǎn)圖
sns.jointplot(x="x", y="y", data=df, stat_func=stats.pearsonr)
03_06.png
（2）六邊圖
data2 = np.random.multivariate_normal(mean, cov, 1000)
df2 = pd.DataFrame(data2, columns=["x", "y"])
with sns.axes_style("white"):
sns.jointplot(x=df2["x"], y=df2["y"], kind="hex", color="k")
03_061.png
（3）相交的 KDE 圖
sns.jointplot(x="x", y="y", data=df2, kind="kde")
03_07.png
（4）地毯圖的 KDE 圖
f, ax = plt.subplots(figsize=(6, 6))
sns.kdeplot(df2.x, df2.y, ax=ax)
sns.rugplot(df2.x, color="#55A868", ax=ax)
sns.rugplot(df2.y, vertical=True, ax=ax, color="#4C72B0")
03_08.png
（5）星云圖
f2, ax = plt.subplots(figsize=(6, 6))
cmap = sns.cubehelix_palette(as_cmap=True, dark=0, light=1, reverse=True)
sns.kdeplot(df2.x, df2.y, cmap=cmap, n_levels=60, shade=True)
03_09.png
（6）加入+圖
g = sns.jointplot(x="x", y="y", data=df, kind="kde", color="#988CBE")
g.plot_joint(plt.scatter, c="w", s=30, linewidth=1, marker="+")
g.ax_joint.collections[0].set_alpha(0)
g.set_axis_labels("$X$", "$Y$")
03_10.png
（7）畫出兩兩變量散點(diǎn)圖
iris = sns.load_dataset("iris")
sns.pairplot(iris)
03_11.png
（8）畫聯(lián)合核函數(shù)
g = sns.PairGrid(iris)
g.map_diag(sns.kdeplot)
g.map_offdiag(sns.kdeplot, cmap="Blues_d", n_levels=6)
03_12.png
