for audio_feature in ['ACOUSTIC', 'DANCE', 'DURATION', 'ENERGY', 'INSTRUMENTAL', 'LIVE', 'LOUD', 'SPEECH', 'TEMPO', 'VALENCE']:
    count, bin_edges = np.histogram(features_df[audio_feature])
    g = sns.FacetGrid(features_df, col='PREFERENCE')
    g.map(plt.hist, audio_feature, color="steelblue", bins=bin_edges) #color="coral"
    plt.show()