Hi all,
I've got a spearman's correlation file generated by 3dfim+, which I then use in a 3dcalc command as follows
3dcalc -a grey_thr2.1_s2.5_tstat1.0.95thresholded.ordered.nii.gz -b sp+tlrc.HEAD -expr "step(a)*b" -prefix killme1
However, the output of this command (killme1) looses its statcode = fico, so I tried to add it back in (I'd like the p-value to be available in the afni window) with 3drefit as follows
3drefit -fico killme1+tlrc.HEAD
but it reports
++ 3drefit: AFNI version=AFNI_2008_07_18_1710 (Jul 23 2008) [64-bit]
++ Authored by: RW Cox
++ Processing AFNI dataset killme1+tlrc.HEAD
** ERROR: Can't change dataset to new type:
* mismatch in number of sub-bricks!
*+ WARNING: Over-writing dataset ./killme1+tlrc.BRIK
++ 3drefit processed 1 datasets
I tried adding the -statpar argument to 3drefit but it made no difference, 3drefit still complains.
3dinfo on the two brik files provided to 3dcalc reports the following:
++ 3dinfo: AFNI version=AFNI_2008_07_18_1710 (Jul 23 2008) [64-bit]
Dataset File: sp+tlrc
Identifier Code: XYZ_rFwYSjLLwu5zKtSs49ScTg Creation Date: Wed Dec 3 14:53:52 2008
Dataset Type: Func-Bucket (-fbuc)
Byte Order: LSB_FIRST [this CPU native = LSB_FIRST]
Storage Mode: BRIK file
Storage Space: 1805258 (1.8 million) bytes
Geometry String: "MATRIX(2,0,0,-90,0,-2,0,126,0,0,2,-72):91,109,91"
Data Axes Tilt: Plumb
Data Axes Orientation:
first (x) = Right-to-Left
second (y) = Posterior-to-Anterior
third (z) = Inferior-to-Superior [-orient RPI]
R-to-L extent: -90.000 [R] -to- 90.000 [L] -step- 2.000 mm [ 91 voxels]
A-to-P extent: -90.000 [A] -to- 126.000 [P] -step- 2.000 mm [109 voxels]
I-to-S extent: -72.000
-to- 108.000 [S] -step- 2.000 mm [ 91 voxels]
Number of values stored at each pixel = 1
-- At sub-brick #0 'Spearman CC' datum type is short: -32767 to 29741 [internal]
[* 2.27563e-05] -0.745655 to 0.676795 [scaled]
statcode = fico; statpar = 33 1 2
----- HISTORY -----
[colmconn@cortex.tchpc.tcd.ie: Wed Dec 3 14:53:52 2008] 3dfim+ -input users_GM_mod_merg_s2.5.nii.gz -mask /usr/support/mrisoftware/fsl-4.0/etc/standard/MNI152_T1_2mm_brain_mask.nii.gz -bucket sp -out 'Spearman CC' -ideal_file abdur.1D
[colmconn@cortex.tchpc.tcd.ie: Wed Dec 3 14:53:52 2008] Output prefix: sp
Dataset File: grey_thr2.1_s2.5_tstat1.0.95thresholded.ordered.nii.gz
Identifier Code: XYZ__X5zYyxalrFcqZ792ejZGQ Creation Date: Wed Jan 7 14:24:05 2009
Dataset Type: Anat Bucket (-abuc)
Byte Order: LSB_FIRST {assumed} [this CPU native = LSB_FIRST]
Storage Mode: NIFTI file
Storage Space: 3610516 (3.6 million) bytes
Geometry String: "MATRIX(2,0,0,-90,0,-2,0,126,0,0,2,-72):91,109,91"
Data Axes Tilt: Plumb
Data Axes Orientation:
first (x) = Right-to-Left
second (y) = Posterior-to-Anterior
third (z) = Inferior-to-Superior [-orient RPI]
R-to-L extent: -90.000 [R] -to- 90.000 [L] -step- 2.000 mm [ 91 voxels]
A-to-P extent: -90.000 [A] -to- 126.000 [P] -step- 2.000 mm [109 voxels]
I-to-S extent: -72.000 -to- 108.000 [S] -step- 2.000 mm [ 91 voxels]
Number of values stored at each pixel = 1
-- At sub-brick #0 '?' datum type is float: 0 to 3
----- HISTORY -----
[colmconn@cortex.tchpc.tcd.ie: Wed Jan 7 14:24:05 2009] 3dmerge -1clust_order 3.6 0 -prefix thresholded/2.1/grey_thr2.1_s2.5_tstat1.0.95thresholded.ordered.nii thresholded/2.1/grey_thr2.1_s2.5_tstat1.0.95thresholded.nii.gz
Does anyone know why 3drefit will not do what I want?
Bye,
Colm.