FDR值的计算方法如下:1)对每个基因进行p-value的计算 假设观测到基因A对应的reads数为x,已知在一个大文库中,每个基因的表达量只占所有基因表达量的一小部分,在这种情况下,p(x)的分布服从泊松分布。已知样本一中唯一比对到基因组的总reads数为N1,样本二中唯一比对到基因组的总reads数为N2,样本...
安装fdrci或fdrtoolpackages运行p.adjust(p,method=‘fdr’,n=length(p))
Q-value,log2 fold change的意思是取log2,fold change的意思是样本质检表达量的差异倍数。