        for c in range(0, nr_cols):
            i = r * nr_cols + c
            if i < len(li_cat_feats):
                sns.boxplot(x=li_cat_feats[i], y=target, data=df, ax = axs[r][c])

    plt.tight_layout()