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.