求教大神：False Discovery Rate（fdr）怎么取有意义值

来源：蜘蛛抓取(WebSpider) 时间：2020-11-11 16:52 标签： 40n120fdr
两列样本数据的差异基因筛选方法：

1）对每个基因进行p-value的计算假设观测到基因A对应的reads数为x已知在一个大文库中，每个基因的表达量只占所有基因表达量的一小部分在這种情况下，p(x)的分布服从泊松分布已知样本一中唯一比对到基因组的总reads数为N1，样本二中唯一比对到基因组的总reads数为N2样本一中唯一比对箌基因A的总reads数为x，样本二中唯一比对到基因A的总reads数为y则基因A在两样本中表达量相等的概率可由以下公式计算：


假设你挑选了R个差异表达嘚基因，其中有S个是真正有差异表达的另外有V个其实是没有差异表达的，是假阳性的实践中希望错误比例Q＝V/R平均而言不能超过某个预先设定的值（比如0.05），在统计学上这也就等价于控制DR不能超过5％.
对所有候选基因的p值进行从小到大排序，则若想控制dr不能超过q则只需找到最大的正整数i，使得 p(i)<= (i*q)/m.然后挑选对应p(1),p(2),...,p(i)的基因做为差异表达基因，这样就能从统计学上保证dr不超过q因此，DR的计算公式如下：
  • 政治敏感、违法虚假信息

我要回帖

更多关于 40n120fdr 的文章

 

随机推荐