Dear all,
I used 3dMEMA but I would like to correct the output.
Can I use ETAC? I am not sure if its only recommended for 3dttest++
If not, which other method should I use?
Many thanks!
Myriam
Dear all,
I used 3dMEMA but I would like to correct the output.
Can I use ETAC? I am not sure if its only recommended for 3dttest++
If not, which other method should I use?
Many thanks!
Myriam
Myriam, ETAC is only available in 3dttest++. You can use 3dClustSim with ACF estimated through individual level or 3dMEMA -residual_Z
Thanks your answer, Gang!
However, when I run it with -residual_Z. The following error pop up
:** option -subs_tstats has illegal dashed parameter: -residual_Z
→ maybe parameter is a mis-typed option?
This is my script:
gen_group_command.py -command 3dMEMA
-write_script cmd.mema.6
-prefix mema++.6_E-S
-dsets path/stats.bil_REML+tlrc.HEAD
-dset_index1_list ‘1,3,4,5,6,7,8,9,10,11,12,13,16,19,20,21,23,27’
-dsets path/stats.bil_REML+tlrc.HEAD
-dset_index1_list ‘1,2,3,4,5,10,11,13,15,16,20,23,30,35,36,37,44,52’
-set_labels SP EN
-subs_betas ‘TASK#0_Coef’
-subs_tstats ‘TASK#0_Tstat’
-residual_Z
Any thoughts what should I do?
This is my afni version:Precompiled binary linux_openmp_64: Feb 10 2017 (Version AFNI_17.0.09)
Try first creating a 3dMEMA script using gen_group_command.py without -residual_Z, and then add -residual_Z to the script.
Thanks Gang! It worked!
However, I don’t understand the output. Now, I have three output files:
-mema+tlrc
-mema_ICC+tlrc
-mema_resZ+tlrc
The first one (mema+tlrc) is still the one uncorrected. But ICC & resZ shows sub-bricks of each participant individually.
How can get the comparison between groups (b & t values) but corrected by multiple comparisons?
Thanks!
Half of the content in resZ is the residuals (the other half is the Z-stats for the residuals). Use 3dFWHM to estimate ACF on the residuals, and then run 3dClustSim with the ACF values to obtain the cluster threshold.
Thank you Gang!
One more question,
-Just to double-check, I run 3dFWHM. are the acf values the ones written by: vary, vary and varzz? so, in the example below will be:
varxx=0.537552 varyy=0.504485 varzz=0.308338
3dFWHM -dset mema++.6_E-S_resZ+tlrc -out E-S_FWHM
user@dancer:./3dFWHM
++ 3dFWHM: AFNI version=AFNI_17.0.09 (Feb 10 2017) [64-bit]
++ Authored by: B. Douglas Ward
*+ WARNING: This program is obsolete! Use 3dFWHMx instead!!
count=1050624
var =0.774703
varxx=0.537552 varyy=0.504485 varzz=0.308338
sx=1.531972 sy=1.593274 sz=2.122862
3dClustSim -acf 0.537552 0.504485 0.308338 -LOTS -prefix 3dClustsim_E_S
-Also, as it says the ‘This program is obsolete. Use 3dFWHMx’
I did try it too by:
3dFWHMx -dset mema++.6_E-S_resZ+tlrc -out E-S_FWHMx
However, I don’t know how which ones are the -acf values, as I thought they have to be values between 0 and 1, and the 3dFWHMx output values are between 3 and 5.
thank you for your time!
You should use 3dFWHMx with option -acf. In addition, I believe that the residual file mema++.6_E-S_resZ+tlrc contrains both the residuals and the Z-stat in an alternating order. So, you need to use only the residuals without the Z-stat part. Try something like the following:
3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0:X:(2)]’ -out E-S_FWHMx
and replace the letter ‘X’ in the above command with the last sub-brick index in the file. Let me know if it works.
Thanks!!
So, I tried what you suggested plus adding the -acf command . But did not work (see below)
I have a total of 35 sub-bricks (see image attached)
3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0:34:(2)]’ -acf -out E-S_FWHMx
++ 3dFWHMx: AFNI version=AFNI_17.0.09 (Feb 10 2017) [64-bit]
++ Authored by: The Bob
** ERROR: selector syntax error 5 ‘:(2)]’
** ERROR: bad sub-brick selector [0:34:(2)]
** FATAL ERROR: Can’t open dataset ‘mema++.6_E-S_resZ+tlrc[0:34:(2)]’
** Program compile date = Feb 10 2017
However, I try this (see below) and it worked:
3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0…$(2)]’ -acf -out E-S_FWHMx
++ 3dFWHMx: AFNI version=AFNI_17.0.09 (Feb 10 2017) [64-bit]
++ Authored by: The Bob
*+ WARNING: removed 131033 voxels from mask because they are constant in time
++ start FWHM calculations
So, my questions are:
-Why is not working when I write: 3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0:34:(2)]’ -acf -out E-S_FWHMx
Am I choosing the wrong sub-bricks number? (see image below)
-is this command ok? 3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0…(2)]' -acf -out E-S_FWHMx I mean, writing [0..(2)] instead of '[0:34:(2)]
-If the command above is ok, the output is a file called: 3dFWHMx.1D
I assume the values in this file are the acf values. Which ones do I need to use for 3dClustSim?
The table looks like:
0 1 1 1
2 0.721748 0.721851 0.742904
2.82843 0.541156 0.545359 0.551906
3.4641 0.417775 0.423102 0.410013
4 0.34554 0.337426 0.3046
4.47214 0.279764 0.276865 0.226289
4.89898 0.233509 0.233667 0.168111
5.65685 0.178671 0.179776 0.0927813
6 0.162576 0.162904 0.0689275
6.32456 0.15554 0.150159 0.0512065
6.63325 0.14018 0.140327 0.0380415
6.9282 0.12457 0.132567 0.0282612
7.2111 0.129147 0.126291 0.0209953
7.48331 0.118861 0.121088 0.0155975
8 0.11475 0.112838 0.00860836
8.24621 0.109159 0.109446 0.00639519
8.48528 0.111247 0.106392 0.00475101
8.7178 0.0983605 0.103604 0.00352954
8.94427 0.101126 0.101029 0.00262211
9.16515 0.0982441 0.0986284 0.00194797
9.38083 0.0904522 0.0963737 0.00144716
9.79796 0.091374 0.0922211 0.000798694
10 0.0875434 0.0902942 0.000593353
10.198 0.0892459 0.0884527 0.000440804
10.3923 0.0866942 0.0866883 0.000327475
10.7703 0.082545 0.0833654 0.000180736
10.9545 0.0827532 0.0817966 0.000134269
11.3137 0.0766195 0.0788235 7.41039e-05
11.4891 0.0783004 0.0774124 5.50521e-05
11.6619 0.0762587 0.0760476 4.08984e-05
11.8322 0.0766266 0.0747266 3.03835e-05
12 0.0749448 0.073447 2.25721e-05
12.1655 0.0792722 0.0722067 1.67688e-05
many many thanks!
My bad! My memory cheated me… Sorry about this!
3dMEMA currently does not provide any output for the residuals, but the Z-stat for the residuals. So, I recommend that you estimate the ACF by using the residuals from the individual level.
Hi Chang,
What do you mean by ‘estimate the ACF by using the residuals from the individual level’?
How should I do that?
Thanks
Myriam
Hi Gang,
Do you mean I should not specify the residuals on the file. So instead of, 3dFWHMx -dset mema++.6_E-S_resZ+tlrc’[0:X:(2)]’ -out E-S_FWHMx
Should be:
3dFWHMx -dset mema++.6_E-S_resZ+tlrc -acf -out E-S_FWHMx
thanks,
Myriam,
Since the residuals at the group level currently cannot be obtained through from 3dMEMA (the option -residual_Z only provides the Z-stat for the residuals), I was suggesting that you estimate the ACF using the residuals from the individual subject analysis from 3dDeconvolve or 3dREMLfit. Specifically, when you perform individual subject analysis, there is an option -errts in 3dDeconvolve that provides the residuals for each subject. Then you run 3dFWHMx with the output from the option -errts in 3dDeconvolve, and the ACF estimates should not vary too much across subjects. You can use the average ACF values to get the cluster thresholds with 3dClustSim. Does this make sense?
Great, Gang! My very last question…
Actually, double checking and correct me if I am wrong, the afni_proc.py provides the ACF values on the blur.errts.1D file, right?
Something like:
3.88582 4.19455 5.32402 4.42727 # epits FWHM blur estimates
0.831582 2.56576 7.90266 6.34748 # epits ACF blur estimates
3.84754 4.1605 5.35565 4.40941 # errts FWHM blur estimates
0.82605 2.54857 7.59102 6.30045 # errts ACF blur estimates
3.83204 4.14254 5.32463 4.38863 # err_reml FWHM blur estimates
0.824958 2.54073 7.48632 6.27801 # err_reml ACF blur estimates
So, I will have to copy the “errts ACF blur estimates” row (0.82605 2.54857 7.59102) onto the 3dClustSim, is this correct?
Can I use the ACF values for 1 participant only, or should I do the average of all participants?
If needs the average, is there any ‘quickly’ way to get the average?
Best,
Thank you for your patience!!
I will have to copy the “errts ACF blur estimates” row (0.82605 2.54857 7.59102) onto the 3dClustSim, is this correct?
Correct.
Can I use the ACF values for 1 participant only, or should I do the average of all participants?
If needs the average, is there any ‘quickly’ way to get the average?
Use the average across the subjects since it is a little more reliable.
great!!
thank you soooo much!
Myriam
Thanks for the command.(tu) I will go through the document of the 3dMEMA and test it on my data first.
The National Institute of Mental Health (NIMH) is part of the National Institutes of Health (NIH), a component of the U.S. Department of Health and Human Services.
Freedom of Information Act |
HHS Vulnerability Disclosure
National Institute of Mental
Health | National Institutes of Health | U.S.
Department of Health and Human Services |