AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

|
December 20, 2022 03:19PM
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 sad smiley" 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
Subject Author Posted

3dMVM thresholding

Matthew Hoptman November 08, 2022 05:00PM

Re: 3dMVM thresholding

gang November 10, 2022 12:17PM

Re: 3dMVM thresholding

Matthew Hoptman November 10, 2022 01:55PM

Re: 3dMVM thresholding

Matthew Hoptman December 13, 2022 09:23AM

Re: 3dMVM thresholding

gang December 15, 2022 09:17AM

Re: 3dMVM thresholding

nuttaon December 20, 2022 03:19PM

Re: 3dMVM thresholding

gang December 21, 2022 08:43PM

Re: 3dMVM thresholding

nuttaon December 22, 2022 12:41PM

Re: 3dMVM thresholding

nuttaon December 22, 2022 03:51PM

Re: 3dMVM thresholding

gang December 23, 2022 01:37PM