Dear Yuki,

The work by Manuel was implemented at the voxel level. 

For the region growing, I would think that if you use 

opt.thre_size = 0;

in the options of the BASC pipeline, the region growing will simply give one voxel per roi, so you will run a voxel level analysis. 

Make sure you specify a `files_in.mask` that is restricted to your area of interest. If you try to run BASC on the full brain at the voxel level you will saturate the memory (unless you have very big voxels).

If you don't want to use the AAL to reduce the memory load, specify your binary mask of interest in `files_in.areas`. See http://niak.simexp-lab.org/pipe_basc.html

Re the block length, check the help of niak_bootstrap_tseries
% BLOCK_LENGTH (if OPT.DGP == 'CBB')
% (integer, default [2*ceil(sqrt(T)) 3*ceil(sqrt(T))]) window width
% used in the circular block bootstrap. If multiple values are
% specified, a random parameter is selected in the list.

In other words, by default the bootstrap uses a block length of 2*the square root of the number of time samples, or 3* the square root of the number of time samples, and this is selected randomly for each bootstrap sample. You can change that parameter using

opt.stability_tseries.sampling.opt.block_length

but you will need to set a single value that applies to all you dataset (i.e. you cannot adjust the block length based on the number of time samples using that method). 

I hope this helps,

Pierre