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  

|
June 05, 2021 08:50AM
Hello experts:
I performed 3dMEMA on five beta effects of 20 people ,the code was shown as below:
3dMEMA \
-prefix $outputdir/all_${taste[$i]} \
-verb 1 \
-max_zeros 0.25 \
-mask $outputdir/all_rigid_mask.Pilot1+tlrc \
-HKtest \
-model_outliers \
-residual_Z \
-jobs 12 \
-set ${taste[$i]}Resp \
1 ${test[0]}[$((3*$i+1))] ${test[0]}[$((3*$i+2))] \
2 ${test[1]}[$((3*$i+1))] ${test[1]}[$((3*$i+2))] \
3 ${test[2]}[$((3*$i+1))] ${test[2]}[$((3*$i+2))] \
4 ${test[3]}[$((3*$i+1))] ${test[3]}[$((3*$i+2))] \
5 ${test[4]}[$((3*$i+1))] ${test[4]}[$((3*$i+2))] \
6 ${test[5]}[$((3*$i+1))] ${test[5]}[$((3*$i+2))] \
7 ${test[6]}[$((3*$i+1))] ${test[6]}[$((3*$i+2))] \
8 ${test[7]}[$((3*$i+1))] ${test[7]}[$((3*$i+2))] \
9 ${test[8]}[$((3*$i+1))] ${test[8]}[$((3*$i+2))] \
10 ${test[9]}[$((3*$i+1))] ${test[9]}[$((3*$i+2))] \
11 ${test[10]}[$((3*$i+1))] ${test[10]}[$((3*$i+2))] \
12 ${test[11]}[$((3*$i+1))] ${test[11]}[$((3*$i+2))] \
13 ${test[12]}[$((3*$i+1))] ${test[12]}[$((3*$i+2))] \
14 ${test[13]}[$((3*$i+1))] ${test[13]}[$((3*$i+2))] \
15 ${test[14]}[$((3*$i+1))] ${test[14]}[$((3*$i+2))] \
16 ${test[15]}[$((3*$i+1))] ${test[15]}[$((3*$i+2))] \
17 ${test[16]}[$((3*$i+1))] ${test[16]}[$((3*$i+2))] \
18 ${test[17]}[$((3*$i+1))] ${test[17]}[$((3*$i+2))] \
19 ${test[18]}[$((3*$i+1))] ${test[18]}[$((3*$i+2))] \
20 ${test[19]}[$((3*$i+1))] ${test[19]}[$((3*$i+2))]
as the "-residual_Z: Output residuals and their Z values used in identifying
outliers at voxel level" told us,but I found the output of 'all_bitter_resZ+tlrc ' files only have the sub-brick of 'res:Z'.As shown below:
-- At sub-brick #0 '01-Res:Z' datum type is short: -26935 to 32767 [internal]
[* 0.00019549] -5.26553 to 6.40563 [scaled]
statcode = fizt
-- At sub-brick #1 '02-Res:Z' datum type is short: -20885 to 32767 [internal]
[* 0.00021256] -4.43932 to 6.96497 [scaled]
statcode = fizt
-- At sub-brick #2 '03-Res:Z' datum type is short: -22407 to 32767 [internal]
[* 0.00019495] -4.36824 to 6.38792 [scaled]
statcode = fizt
-- At sub-brick #3 '04-Res:Z' datum type is short: -25425 to 32767 [internal]
[* 0.000197025] -5.00935 to 6.45591 [scaled]
statcode = fizt
-- At sub-brick #4 '05-Res:Z' datum type is short: -26215 to 32767 [internal]
[* 0.000175355] -4.59692 to 5.74585 [scaled]
statcode = fizt
-- At sub-brick #5 '06-Res:Z' datum type is short: -23350 to 32767 [internal]
[* 0.000192913] -4.50452 to 6.32118 [scaled]
statcode = fizt
-- At sub-brick #6 '07-Res:Z' datum type is short: -26924 to 32767 [internal]
[* 0.00016495] -4.44112 to 5.40492 [scaled]
statcode = fizt
-- At sub-brick #7 '08-Res:Z' datum type is short: -30550 to 32767 [internal]
[* 0.000151262] -4.62105 to 4.9564 [scaled]
statcode = fizt
-- At sub-brick #8 '09-Res:Z' datum type is short: -26838 to 32767 [internal]
[* 0.000167578] -4.49745 to 5.49102 [scaled]
statcode = fizt
-- At sub-brick #9 '10-Res:Z' datum type is short: -21694 to 32767 [internal]
[* 0.000197313] -4.2805 to 6.46535 [scaled]
statcode = fizt
-- At sub-brick #10 '11-Res:Z' datum type is short: -19951 to 32767 [internal]
[* 0.000218496] -4.35921 to 7.15945 [scaled]
statcode = fizt
-- At sub-brick #11 '12-Res:Z' datum type is short: -20442 to 32767 [internal]
[* 0.000235805] -4.82032 to 7.72661 [scaled]
statcode = fizt
-- At sub-brick #12 '13-Res:Z' datum type is short: -20568 to 32767 [internal]
[* 0.000216673] -4.45653 to 7.09973 [scaled]
statcode = fizt
-- At sub-brick #13 '14-Res:Z' datum type is short: -21314 to 32767 [internal]
[* 0.000229566] -4.89298 to 7.5222 [scaled]
statcode = fizt
-- At sub-brick #14 '15-Res:Z' datum type is short: -18763 to 32767 [internal]
[* 0.000227736] -4.27301 to 7.46223 [scaled]
statcode = fizt
-- At sub-brick #15 '16-Res:Z' datum type is short: -26444 to 32767 [internal]
[* 0.000216411] -5.72277 to 7.09114 [scaled]
statcode = fizt
-- At sub-brick #16 '17-Res:Z' datum type is short: -27842 to 32767 [internal]
[* 0.00016686] -4.64572 to 5.46751 [scaled]
statcode = fizt
-- At sub-brick #17 '18-Res:Z' datum type is short: -29683 to 32767 [internal]
[* 0.000169653] -5.03582 to 5.55903 [scaled]
statcode = fizt
-- At sub-brick #18 '19-Res:Z' datum type is short: -32767 to 29361 [internal]
[* 0.000193574] -6.34286 to 5.68354 [scaled]
statcode = fizt
-- At sub-brick #19 '20-Res:Z' datum type is short: -25542 to 32767 [internal]
[* 0.000198609] -5.07287 to 6.50781 [scaled]
statcode = fizt

Is there some problems with my code ,why is there no residuals results of 3dmema output? Because I want to use the Output residuals to run 3dClustSim.Thanks for your help !
Subject Author Posted

why is there no residuals results of 3dmema output?

zdlu June 05, 2021 08:50AM

Re: why is there no residuals results of 3dmema output?

gang June 05, 2021 01:07PM