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 !