[問題] FDR P-value修正

[purplehui (生命的強度表現在不依賴)]

我想用FDR做P-value修正

將P-value排序之後，利用mafdr函數:

C=textread('CCP-value.txt');

[FDR,Q]=mafdr(C,'showplot',true);

跑出兩個圖

一是PI0(λ)對λ，另一個是q-value對p-value

但我需要的是p-value vs. FDR baseline

也就是 p-value對k/N

請問該做哪些修改

另外，想請問PI0(λ)是代表甚麼？　


謝謝！！

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.126.132.96
推 skylikewater:推一下　我也在困擾 FDR
→ purplehui:http://ppt.cc/-c7O 我有找到這個，裡面有FDR.m可以用，
→ purplehui:但....依舊畫不出圖來=口= 救命！