import numpy as np

g = sns.PairGrid(iris_df, hue="variety")
g.map_upper(sns.scatterplot)
g.map_lower(sns.kdeplot,fill=True,)
g.map_diag(sns.kdeplot,fill=True)