Hi All,
I've been doing whole-brain functional analyses using glmfit and mri_volcluster. As I understand it, correction for multiple comparisons is done on a cluster level with this command, by specifying a minimum cluster size and then applying a threshold within each cluster. This is working ok, but I'd like to compare it with another method of multiple comparison correction. Is there a command that will allow me to do multiple comparison correction on a voxel level, such as family-wise error, by adjusting the threshold for each voxel?
Thanks, Joe Andreano
_________________________________
Joseph M. Andreano, PhD Research Fellow Interdisciplinary Affective Science Lab Massachusetts General Hospital _________________________________
Nothing that is easy to do. You can compute your own voxel-wise Bonferroni threshold based on the search space and then use mri_binarize to do the thresholding. I don't have a command-line that computes FDR, but you can get the threshold for a given FDR in matlab with fast_fdrthresh, then apply the threshold using mri_binarize.
doug
Joe Andreano wrote:
Hi All,
I've been doing whole-brain functional analyses using glmfit and mri_volcluster. As I understand it, correction for multiple comparisons is done on a cluster level with this command, by specifying a minimum cluster size and then applying a threshold within each cluster. This is working ok, but I'd like to compare it with another method of multiple comparison correction. Is there a command that will allow me to do multiple comparison correction on a voxel level, such as family-wise error, by adjusting the threshold for each voxel?
Thanks, Joe Andreano
Joseph M. Andreano, PhD Research Fellow Interdisciplinary Affective Science Lab Massachusetts General Hospital _________________________________
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
freesurfer@nmr.mgh.harvard.edu