function selected_indices = cosmo_anova_feature_selector(dataset, how_many)
cosmo align hdr
cosmo average samples hdr