Hi AFNI folks,
I discovered (I think) a small bug when using the 3dANOVA programs with niml surface files. The error comes at the end of the std output in cases where I've included a -bucket option for putting all calculated sub-bricks into one bucket. Here's the command I ran and the associated error:
3dANOVA2 -type 3 -alevels 6 -blevels 9 \
-dset 1 1 BK_ColrL_LH.niml.dset \
-dset 1 2 CMS_ColrL_LH.niml.dset \
...
-dset 6 8 PM_HoldR_LH.niml.dset \
-dset 6 9 RS_HoldR_LH.niml.dset \
-acontr 1 0 0 0 -1 0 ColrL-HL_LH.niml.dset \
-acontr 0 1 0 0 0 -1 ColrR-HR_LH.niml.dset \
-acontr 0 0 1 0 -1 0 LocnL2R-HL_LH.niml.dset \
-acontr 0 0 0 1 0 -1 LocnR2L-HR_LH.niml.dset \
-bucket AttnSW_grpsrf_LH.niml.dset
++ 3dANOVA: AFNI version=AFNI_2008_02_01_1144 (Feb 19 2008) [32-bit]
++ Authored by: B. Douglas Ward
Data set dimensions: nx = 168118 ny = 1 nz = 1 nxyz = 168118
** Changes have been made for 3dANOVA2 computations of type 3.
For details, please see:
[
afni.nimh.nih.gov]
++ Writing combined dataset into ./ColrL-HL_LH.niml.dset+orig.HEAD
+ created 1 FDR curves in header
++ Writing combined dataset into ./ColrR-HR_LH.niml.dset+orig.HEAD
+ created 1 FDR curves in header
++ Writing combined dataset into ./LocnL2R-HL_LH.niml.dset+orig.HEAD
+ created 1 FDR curves in header
++ Writing combined dataset into ./LocnR2L-HR_LH.niml.dset+orig.HEAD
+ created 1 FDR curves in header
Writing `bucket' dataset into AttnSW_grpsrf_LH.niml.dset
RUNNING COMMAND: 3dbucket -prefix AttnSW_grpsrf_LH.niml.dset ./ColrL-HL_LH.niml.dset+orig.HEAD ./ColrR-HR_LH.niml.dset+orig.HEAD ./LocnL2R-HL_LH.niml.dset+orig.HEAD ./LocnR2L-HR_LH.niml.dset+orig.HEAD
++ 3dbucket: AFNI version=AFNI_2008_02_01_1144 (Feb 19 2008) [32-bit]
can't open dataset ./ColrL-HL_LH.niml.dset+orig.HEAD
RUNNING COMMAND: 3drefit -sublabel 0 ColrL-HL_LH.niml.dset:Contr -sublabel 1 ColrL-HL_LH.niml.dset:t-stat -sublabel 2 ColrR-HR_LH.niml.dset:Contr -sublabel 3 ColrR-HR_LH.niml.dset:t-stat -sublabel 4 LocnL2R-HL_LH.niml.dset:Contr -sublabel 5 LocnL2R-HL_LH.niml.dset:t-stat -sublabel 6 LocnR2L-HR_LH.niml.dset:Contr -sublabel 7 LocnR2L-HR_LH.niml.dset:t-stat ./AttnSW_grpsrf_LH.niml.dset+orig.HEAD
++ 3drefit: AFNI version=AFNI_2008_02_01_1144 (Feb 19 2008) [32-bit]
++ Authored by: RW Cox
** ERROR: Can't open dataset ./AttnSW_grpsrf_LH.niml.dset+orig.HEAD
++ 3drefit processed 0 datasets
Looks like the 3dbucket command has the standard AFNI '+orig.HEAD' appended to the 'niml.dset' files and, of course, can't find anything like that. As a result, no bucket is written and the 3drefit command fails, as well.
Thanks,
Adam