Hi Kersten,

I used this code to calculate the threshold:

P = [ F_lhstats1.pval(lhcortex) F_rhstats1.pval(rhcortex) ];

G = [ F_lhstats1.sgn(lhcortex) F_rhstats1.sgn(rhcortex) ];

[detvtx,sided_pval,pth] = lme_mass_FDR2(P,G,[],0.05,0);

pcor = -log10(pth)

I substituted 0.1, 0.2, 0.5, 0.8, and 1 where the 0.05 is, and I kept getting the same pcor value of 4.8943.

Jessica Hua