function ds_stat = cosmo_montecarlo_phase_stat(ds, varargin)
cosmo montecarlo cluster stat hdr
cosmo naive bayes classifier searchlight hdr