RF+BK: change function signature for stats (make it dataset-only), support returning z-scores and p-values
RF+BK: use spherical neighborhood instead of voxelselection