function compute_tfce2( DS, nbrh, Opt, tempath)


for nmeth = 1 : length( Opt)
    
    % run tfce (across ROI sizes)
    test_DS     =   cosmo_montecarlo_cluster_stat( DS, nbrh, Opt(nmeth));
    pvals       =   1 - normcdf( test_DS.samples);
    
    % save data
    save( tempath{nmeth}, 'pvals')

end