Dear all,
I would like to run 3dGroupInCorr for a comparison between 2 groups (2-sample t-test with no covariates). For that, I did as shown in slide 23 of
[
afni.nimh.nih.gov]
# Create GroupInCorr structures for each Group
ggg=
for subject in $subjects_FOGp
do
ggg=$(echo $ggg pb05.${subject}.anaticor.blur+tlrc.HEAD )
done
3dSetupGroupInCorr -mask full_mask.frac0.85_TT_caez_ml_18+tlrc. -prefix FOGp -byte -prep DEMEAN $ggg
ggg=
for subject in $subjects_FOGm
do
ggg=$(echo $ggg pb05.${subject}.anaticor.blur+tlrc.HEAD )
done
3dSetupGroupInCorr -mask full_mask.frac0.85_TT_caez_ml_18+tlrc. -prefix FOGm -byte -prep DEMEAN $ggg
where, in my opinion, the .grpincorr.niml files were created correctly:
[ccaballero@cajal02 GroupInCorr]$ 3dSetupGroupInCorr -mask full_mask.frac0.85_TT_caez_ml_18+tlrc. -prefix CONt -byte -prep DEMEAN $ggg
++ 3dSetupGroupInCorr: AFNI version=AFNI_16.0.01 (Jan 27 2016) [64-bit]
++ Authored by: RW Cox
++ Number of voxels in mask = 39948
++ Starting scan through 24 datasets
++ Dataset 1/24: ./pb05.001_Park_Pda_Cm.anaticor.blur+tlrc.BRIK input = 3679 ms; DEMEAN/normalize = 143 ms; output = 22 ms
++ Dataset 2/24: ./pb05.002_Park_Pjo_Cb.anaticor.blur+tlrc.BRIK input = 4067 ms; DEMEAN/normalize = 168 ms; output = 8 ms
++ Dataset 3/24: ./pb05.002_Park_Mab_Cm.anaticor.blur+tlrc.BRIK input = 2925 ms; DEMEAN/normalize = 163 ms; output = 9 ms
++ Dataset 4/24: ./pb05.003_Park_Cue_Cb.anaticor.blur+tlrc.BRIK input = 2554 ms; DEMEAN/normalize = 167 ms; output = 9 ms
++ Dataset 5/24: ./pb05.004_Park_Jibm_Cb.anaticor.blur+tlrc.BRIK input = 3221 ms; DEMEAN/normalize = 170 ms; output = 8 ms
++ Dataset 6/24: ./pb05.004_Park_Czs_Cm.anaticor.blur+tlrc.BRIK input = 2655 ms; DEMEAN/normalize = 171 ms; output = 8 ms
++ Dataset 7/24: ./pb05.005_Park_Cgg_Cm.anaticor.blur+tlrc.BRIK input = 2575 ms; DEMEAN/normalize = 165 ms; output = 7 ms
++ Dataset 8/24: ./pb05.005_Park_Trc_Cb.anaticor.blur+tlrc.BRIK input = 2668 ms; DEMEAN/normalize = 166 ms; output = 7 ms
++ Dataset 9/24: ./pb05.007_Park_Jalv_Cm.anaticor.blur+tlrc.BRIK input = 2486 ms; DEMEAN/normalize = 170 ms; output = 9 ms
++ Dataset 10/24: ./pb05.008_Park_Sci_Cb.anaticor.blur+tlrc.BRIK input = 2677 ms; DEMEAN/normalize = 170 ms; output = 7 ms
++ Dataset 11/24: ./pb05.009_Park_Aee_Cb.anaticor.blur+tlrc.BRIK input = 2410 ms; DEMEAN/normalize = 167 ms; output = 8 ms
++ Dataset 12/24: ./pb05.010_Park_Mlzf_Cm.anaticor.blur+tlrc.BRIK input = 3054 ms; DEMEAN/normalize = 135 ms; output = 7 ms
++ Dataset 13/24: ./pb05.011_Park_Evl_Cb.anaticor.blur+tlrc.BRIK input = 2764 ms; DEMEAN/normalize = 161 ms; output = 8 ms
++ Dataset 14/24: ./pb05.011_Park_Mcg_Cm.anaticor.blur+tlrc.BRIK input = 3802 ms; DEMEAN/normalize = 132 ms; output = 7 ms
++ Dataset 15/24: ./pb05.012_Park_Ccp_Cm.anaticor.blur+tlrc.BRIK input = 3786 ms; DEMEAN/normalize = 162 ms; output = 5 ms
++ Dataset 16/24: ./pb05.012_Park_Mcai_Cb.anaticor.blur+tlrc.BRIK input = 2977 ms; DEMEAN/normalize = 160 ms; output = 7 ms
++ Dataset 17/24: ./pb05.013_Park_Msg_Cm.anaticor.blur+tlrc.BRIK input = 2655 ms; DEMEAN/normalize = 167 ms; output = 8 ms
++ Dataset 18/24: ./pb05.013_Park_Jlul_Cb.anaticor.blur+tlrc.BRIK input = 2669 ms; DEMEAN/normalize = 168 ms; output = 7 ms
++ Dataset 19/24: ./pb05.014_Park_Bga_Cb.anaticor.blur+tlrc.BRIK input = 2813 ms; DEMEAN/normalize = 166 ms; output = 8 ms
++ Dataset 20/24: ./pb05.015_Park_Fao_Cb.anaticor.blur+tlrc.BRIK input = 3098 ms; DEMEAN/normalize = 136 ms; output = 6 ms
++ Dataset 21/24: ./pb05.016_Park_Amea_Cb.anaticor.blur+tlrc.BRIK input = 2795 ms; DEMEAN/normalize = 167 ms; output = 7 ms
++ Dataset 22/24: ./pb05.017_Park_Amg_Cb.anaticor.blur+tlrc.BRIK input = 2591 ms; DEMEAN/normalize = 161 ms; output = 7 ms
++ Dataset 23/24: ./pb05.019_Park_Cll_Cb.anaticor.blur+tlrc.BRIK input = 3157 ms; DEMEAN/normalize = 132 ms; output = 6 ms
++ Dataset 24/24: ./pb05.020_Park_Jmer_Cb.anaticor.blur+tlrc.BRIK input = 3726 ms; DEMEAN/normalize = 135 ms; output = 8 ms
++ Wrote data file CONt.grpincorr.data = 284,070,228 bytes (about 284 million [mega])
+ Total elapsed time = 75887 ms = 1.26 min
++ Wrote head file CONt.grpincorr.niml
Finally, I run 3dGroupInCorr and then open AFNI -niml in another terminal,
3dGroupInCorr -setA FOGp.grpincorr.niml -setB FOGm.grpincorr.niml -verb
However, although I can overlay the GroupInCorr, the values are all zero for all created datasets, i.e. both max and min in the overlay are 0. What could it be the issue here?
Thanks very much in advance.
Best,
Cesar