function ds_z = cosmo_montecarlo_cluster_stat(ds, nbrhood, varargin)
cosmo meta feature selection classifier hdr
cosmo montecarlo phase stat hdr