現在、ある化合物を処理/未処理した癌細胞のRNA-seqデータに対して、MSigDBの遺伝子セットを用いてGSEAを行っています
そこで例えば"HALLMARK_KRAS_SIGNALING_DN"のように「KRASが活性化していると発現が『低下する』遺伝子セット」のエンリッチメントスコアが『正の値』を取る場合、結局その遺伝子群は発現低下(= KRASが活性化)しているのかどうかがいまいちよくわかりません
その遺伝子セットの中で発現が低下している遺伝子が多ければエンリッチメントスコアは負になるような気がするのですが、どう判別すればいいでしょうか |
|