afni.nimh.nih.gov

Hi Gang,

I have similar questions. I ran 3dMVM with -resid flag on n=74. Here is my script.

3dMVM -prefix GvsHC.nii.gz \
    -jobs 4 \
    -bsVars "Group*CIRSGcenter+Education+Age+Gender" \
    -wsVars "Valence*Magnitude" \
    -qVars  "Age,Education,CIRSGcenter" \
    -qcenter 63.5, 16.6, 0 \
    -resid   GvsHC_resid.nii.gz \
    -num_glt 18 \
    -gltLabel 1 3Way_Group_by_ValenceandMagnitude   -gltCode 1  "Group : 1*GRIEF -1*HC    Valence : 1*Reward -1*Punishment     Magnitude : 1*Lg -1*Small" \
    -gltLabel 2 2Way_Group_Valence                  -gltCode 2  "Group : 1*GRIEF -1*HC    Valence : 1*Reward -1*Punishment     Magnitude : 1*All" \
    -gltLabel 3 posthoc_Group_Reward                -gltCode 3  "Group : 1*GRIEF -1*HC    Valence : 1*Reward                   Magnitude : 1*All" \
    -gltLabel 4 posthoc_Group_Punishment            -gltCode 4  "Group : 1*GRIEF -1*HC    Valence : 1*Punishment               Magnitude : 1*All" \
    -gltLabel 5 posthoc_Group_RewardMagnitude       -gltCode 5  "Group : 1*GRIEF -1*HC    Valence : 1*Reward                   Magnitude : 1*Lg -1*Small" \
    -gltLabel 6 posthoc_Group_PunishmentMagnitude   -gltCode 6  "Group : 1*GRIEF -1*HC    Valence : 1*Punishment               Magnitude : 1*Lg -1*Small" \
    -gltLabel 7 posthoc_Grief_Reward		    -gltCode 7  "Group : 1*GRIEF          Valence : 1*Reward                   Magnitude : 1*All" \
    -gltLabel 8 posthoc_Grief_Punishment	    -gltCode 8  "Group : 1*GRIEF          Valence : 1*Punishment               Magnitude : 1*All" \
    -gltLabel 9 posthoc_HC_Reward		    -gltCode 9  "Group : 1*HC             Valence : 1*Reward                   Magnitude : 1*All" \
    -gltLabel 10 posthoc_HC_Punishment	            -gltCode 10  "Group : 1*HC            Valence : 1*Punishment               Magnitude : 1*All" \
    -gltLabel 11 posthoc_Grief_LgReward		    -gltCode 11  "Group : 1*GRIEF         Valence : 1*Reward                   Magnitude : 1*Lg" \
    -gltLabel 12 posthoc_Grief_LgPunishment	    -gltCode 12  "Group : 1*GRIEF         Valence : 1*Punishment               Magnitude : 1*Lg" \
    -gltLabel 13 posthoc_HC_LgReward		    -gltCode 13  "Group : 1*HC            Valence : 1*Reward                   Magnitude : 1*Lg" \
    -gltLabel 14 posthoc_HC_LgPunishment	    -gltCode 14  "Group : 1*HC            Valence : 1*Punishment               Magnitude : 1*Lg" \
    -gltLabel 15 posthoc_Grief_SmallReward          -gltCode 15  "Group : 1*GRIEF         Valence : 1*Reward                   Magnitude : 1*Small" \
    -gltLabel 16 posthoc_Grief_SmallPunishment	    -gltCode 16  "Group : 1*GRIEF         Valence : 1*Punishment               Magnitude : 1*Small" \
    -gltLabel 17 posthoc_HC_SmallReward		    -gltCode 17  "Group : 1*HC            Valence : 1*Reward                   Magnitude : 1*Small" \
    -gltLabel 18 posthoc_HC_SmallPunishment	    -gltCode 18  "Group : 1*HC            Valence : 1*Punishment               Magnitude : 1*Small" \
    -dataTable \

Then I ran 3dFWHMx with -acf, and the residual from 3dMVM as an input.

3dFWHMx -acf -mask $mask -input GvsHC_resid.nii.gz -out acf_3dMVM2

Here is the output (444 lines)
10.6753 9.84336 9.79558
10.1204 9.35333 9.63601
10.6 10.5259 10.0315
.
.
.
11.3669 9.93517 10.3598

When I plugged these number in 3dClustsim. I got an error message "[-acf: 'a' value should be between 0 and 1 " b]

3dClustSim -mask $mask -acf mean_col1 mean_col2 mean_col3 -pthr 0.05 0.01 00.005 0.001 -iter 10000 >>ClustSim_3dMVMresid_${mask}_n74.txt

Could you please advise?

Thank you,
Nuttaon