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

参考链接：

https://stat.ethz.ch/pipermail/bioconductor/2014-January/056992.html

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