學習筆記49—matlab FDR校正

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