Hi all

I came across a puzzling behavior of 3dclust. When I run the following:

3dclust -1thresh 3.30 wo.irf.t.thr_0.001_2 \
4.3 54 'stats.irf+orig[1]' > wo_irf_clust.txt

it says:
++ Writing dataset ./wo.irf.t.thr_0.001_2+orig.HEAD

but the above output dataset is nowhere to be found! (although the text file with the cluster table is correctly written)


Identical results if I specify 1dindex and 1tindex:

3dclust -1thresh 3.30 -1dindex 0 -1tindex 0 -prefix wo.irf.t.thr_0.001_2 \
4.3 54 'stats.irf+orig[1]' > wo_irf_clust.txt


On the other hand, if I don't include the subbrick specification in the input dataset, but indicate it in 1dindex and 1tindex, as in:

3dclust -1thresh 3.30 -1dindex 1 -1tindex 1 -prefix wo.irf.t.thr_0.001_2 \
4.3 54 stats.irf+orig[1] > wo_irf_clust.txt

I do get a dataset, but instead of a single-brik dataset I get a bucket dataset with the same number of subbricks as the input dataset, even though in the help of 3dmerge, with regard to the -1dindex and -1tindex options, is stated:

"With these, you can operate on any given sub-brick of the inputs dataset(s) to produce *** as output a 1 brick dataset. ***" (asterisks added).


any comments or suggestions?

thanks!

giuseppe