The bar plots in steps 3 and 5 could have been created with the more advanced seaborn PairGrid constructor, which can plot a bivariate relationship. Using PairGrid is a two-step process. The first call to PairGrid prepares the grid by alerting it to which variables will be x and which will be y. The second step applies a plot to all of the combinations of x and y columns:

>>> g = sns.PairGrid(diamonds,size=5,
                 x_vars=["color", "cut", "clarity"],
                 y_vars=["price"])
>>> g.map(sns.barplot)
>>> g.fig.suptitle('Replication of Step 3 with PairGrid', y=1.02)