#######个人理解
没有人为设定cut off值,而是首先进行差异分析,将所有基因Fold Change的log值进行从正到负排序,建立一个gene list(假设有10000个基因)。
选择要研究的基因(eg:FBXO7),假设与FBXO7相关的基因有50个(基因F),将上述gene list中的基因一个个拿过来比对,如果list中的这个基因L在基因F1家族之中,那就给这个基因F1加分,如果不是,就减分。
因为list中的基因已经排好序了,处于同一家族的基因位置靠近,当这一段list在基因F1的家族中时,基因F1的打分会不断增加,直至处于峰值(>0),这段list gene比对完之后,由于之后的list中gene在基因F家族中的不断减少,F的打分也会不断降低。
这样就可以分析基因F1正向富集在gene list的哪一段,同理也可知道F1负向富集(抑制)的情况。
#########################
网友评论