subset_by function

Subset spectra by factor