Multivariate faceting
g = sb.FacetGrid(data = df, col = 'cat_var1', size = 4) g.map(sb.boxplot, 'cat_var2', 'num_var2')