Well, now I know that '3dinfo' provides more information about the 'stats.${subject}_REML+orig' images/bricks than its AFNI header. Thus, I see that some bricks are coefficient maps rather than statistic maps. Yesterday I was thinking that all 9 bricks in my 'stats.${subject}_REML+orig' are F- and t- statistic maps, some of which referred to the motion covariates... Now I know that there are no images/bricks referring to the motion covariates.

I am only interested in extracting the statistic map corresponding to the HRF amplitude. Yesterday I called this statistic map the "ultimate" one. I am still not sure whether its the statistic map in brick [2]:

-- At sub-brick #0 'activation_stimulus#0_Tstat' datum type is float: -6.67307 to 22.7179
statcode = fitt; statpar = 88

or the one in brick [7]:

-- At sub-brick #0 'activity#0_Tstat' datum type is float: -8.30011 to 19.8636
statcode = fitt; statpar = 88

They are very similar to each other. Now I also know where the degrees of freedom are saved and how to convert a t-statistic map to a z-statistic map in AFNI:

[afni.nimh.nih.gov]

I am not sure what is the difference between bricks [2] and [7] in my 'stats.${subject}_REML+orig'. Could you please give a hint?

Best wishes,

Wiktor Olszowy