On 01/06/2016 11:00 AM, stdp82@virgilio.it wrote:

Hi list,

I would like to ask you a number of questions regarding the FDR which is performed following seed-based fcMRI using FS-FAST. 1-In the previous mail you have suggested me to use mir_fdr command line (version 6.0). Because I have some problem with new version license (as reported in previous mail) and, overall I have run my data using version 5.3, I would ask you whether FDR integrated in tksurfer (version 5.3) is however reliable as the FDR included in 6.0?

Yes, they call the same code.

2-Becouse the FDR analysis was performed separately for the right and left hemispheres and the subcortical structures, should I set the FDR threshold to 0.05/3 in order to correct for multiple comparison? (I have one group)

I'm not sure. You can try to run both hemispheres in matlab, somethiing like lh = MRIread('sig.lh.mgh'); lhp = 10^-abs(lh.vol); rh = MRIread('sig.rh.mgh'); rlhp = 10^-abs(rh.vol);

p = [lhp rhp]; fdrthresh = fast_fdrthresh(p,fdr);

3- Is easy for me to follow the guidelines and visualize the cortical areas which are functionally connected with the omolateral seed. How can I do to visualize also the the cortical areas which are functionally connected with the controlateral seed? If change only the rh/lh in command line the visualization of the results is wrong.

I don't understand.

4- I'd like to perform mri_surfcluster on a group to eliminate cluster < 5 voxels from sig.nii.gz after FDR correction. I'm running mri_surfcluster --in sig.nii.gz --hemi rh --sign abs --fdr 0.05 --minarea 5 --sum summaryfile --ocn ocnid but --subject ... must be included. How should do I do? I would like to run my analysis on group results, not on a single subject.

If this is from a group analysis, the subject is probablly fsaverage

Thanks for your work. Regards,

Stefano

---

Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer