http://home.52brain.com/forum.php?mod=viewthread&tid=27066&page=1#pid170857

http://www.mathworks.com/matlabcentral/answers/144113-mafdr-interpreting-q-values-vs-bhfdr-adjusted-p-values

matlab自带函数mafdr,当ttest数较多时,可直接用[FDR, Q]=mafdr(P);但是Storey procedure在p值少于1000个时会崩溃,此时应改用BH FDR方法:mafdr(P,'BHFDR', true)。该方法对于少量ttest更稳健,但是更保守。

04-18 02:34