Hello-
We are trying to use easythresh to do family-wise error correction based for level two analyses. The stat maps were generated in another program for an experiment that had 27 subjects. We would appreciate some feedback to see if the following approach is valid.
1) First we obtain smoothness estimates for each subject’s residual time-series following first level analyses (time-series = 149 images)
smoothest –d 148 –r sub1.resid.nii.gz –m ave_mask -V
…..
…..
…..
smoothest –d 148 –r sub27.resid.nii.gz –m ave_mask –V
2) We then average the temporally corrected DLH and RESELS values from all subjects to obtain a single DLH ($ave_DLH) and RESEL ($ave_RESEL) for the experiment. Is this similar to what is done in FEAT for level two analyses?
3) Finally, we use easythresh to correct for false positives, after first hard-coding our derived DLH and RESEL values from step 2 into the easythresh script.
|