cosmo anova feature selector hdr

function selected_indices = cosmo_anova_feature_selector(dataset, how_many)