Hi!

I am using 3dClustSim to conduct cluster correction for an analysis run using FSL’s FEAT. I need to estimate the smoothness of my data using 3dFWHMx.

So in Afni, when you complete first level analysis with afni_proc.py, you get a residual file typically called errts+tlrc. This is the residual of the raw time series after you account for all your model parameters, and what is recommended to use as input to estimate the smoothness of the data using 3dFWHMx

In our lab we use fsl FEAT, which also outputs a residual file called res4d.nii

BUT there is an important difference in how first level analyses with multiple runs are typically analyzed in afni and fsl. In afni, you concatenate the timeseries, then run the analysis. In fsl, you run the analysis separately for each run, then average those analyses in a gfeat. You therefore have a res4d.nii file for each run (in this task 4 runs) and one for the gfeat model that averages the models for each run.

The res4d.nii files from each of the 4 run are residual timeseries of that run.

The res4d.nii files from the gfeat are the residuals of the model averaging the 4 runs, so they just contain 4 data points per voxel that sum to 0.

My question is: what should I use as input to 3dFWHMx?

Thanks!

David