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更穩健，但是更保守。