Dear Professor Reuter,

Thank you for your reply and recommendation for computing the single FDR. I could solve my problem as follow your explanation.

Best Regards,

Han.

On Fri, Jan 1, 2016 at 12:59 PM, Martin Reuter mreuter@nmr.mgh.harvard.edu wrote:

Hi Han,

I think detvtx is the detected vertices (a binary vector with vertices were you can reject null after FDR). The 3rd argument (pth) is the FDR threshold. You can compute pcor = -log10(pth) from it and use it as a threshold when plotting the p-value map to limit the plot to only significant values.

I also recommend to concat the left and right hemi results and compute a single FDR threshold. That way you won't need different color bars for each hemisphere:

P = [ F_lhstats.pval(lhcortex) F_rhstats.pval(rhcortex) ]; G = [ F_lhstats.sgn(lhcortex) F_rhstats.sgn(rhcortex) ]; [detvtx,sided_pval,pth] = lme_mass_FDR2(P,G,[],0.05,0); pcor = -log10(pth)

Best, Martin

On 12/31/2015 05:33 PM, Hanbyul Cho wrote:

Dear FreeSurfer Team,

My name is Han. I am wondering how I can save the FDR2 corrected results.

I processed in MATLAB for Linear mixed effects models analysis, as follow:

[lhTh01,lhRe01] = lme_mass_fit_EMinit(X,[1],Y,ni,lhcortex,3,4) [lhRgs01,lhRgMeans01] = lme_mass_RgGrow(lhsphere, lhRe01, lhTh01,lhcortex,2,95); lhstats01 = lme_mass_fit_Rgw(X,[1],Y,ni, lhTh01, lhRgs01, lhsphere,[],'euc','exp',4)

F_lhstats01=lme_mass_F(lhstats01,CM,4); dvtx = lme_mass_FDR2(F_lhstats01.pval, F_lhstats01.sgn, lhcortex, 0.05, 0); length(dvtx) ans = 8000

In this case, I search the corrected results that were 8000 vertices. I think the dvtx value was for saving the significant each vertex numbering label, is it right? However, I don't know the next step for save the 8000 corrected vertices.

fs_write_fstats(F_lhstats01,mri,'%%$.mgh','sig'); I can save the whole results 'F_lhstats01' as the mgh format. However these were uncorrected results, is it right?

Could you let me know the method for saving the FDR2 corrected results as a mgh format?

Best Regards,

Han.

---

Freesurfer mailing listFreesurfer@nmr.mgh.harvard.eduhttps://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

-- Martin Reuter, PhD Assistant Professor of Radiology, Harvard Medical School Assistant Professor of Neurology, Harvard Medical School A.A.Martinos Center for Biomedical Imaging Massachusetts General Hospital Research Affiliate, CSAIL, MIT Phone: +1-617-724-5652 Web : http://reuter.mit.edu

---

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

The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.