function ds_avg = cosmo_average_samples(ds, varargin)
cosmo anova feature selector hdr
cosmo balance dataset hdr