什么是GSEA?如何解讀結果

什么是GSEA?

Gene Set Enrichment Analysis (基因集富集分析)用來評估一個預先定義的基因集的基因在與表型相關度排序的基因表中的分布趨勢,從而判斷其對表型的貢獻。其輸入數據包含兩部分,一是已知功能的基因集 (可以是GO注釋、MsigDB的注釋或其它符合格式的基因集定義),一是表達矩陣 (也可以是排序好的列表),軟件會對基因根據其與表型的關聯度(可以理解為表達值的變化)從大到小排序,然后判斷基因集內每條注釋下的基因是否富集于表型相關度排序后基因表的上部或下部,從而判斷此基因集內基因的協同變化對表型變化的影響。

分析結果如下:

GS:基因集的名字,GO條目的名字

SIZE:GO條目中包含表達數據集文中的基因數目(經過條件篩選后的值);

ES:富集評分;

NES:校正后的歸一化的ES值。由于不同用戶輸入的基因數據庫文件中的基因集數目可能不同,富集評分的標準化考慮了基因集個數和大小。其絕對值大于1為一條富集標準。

NOM p-val:即p-value,是對富集得分ES的統計學分析,用來表征富集結果的可信度;

FDR q-val:即q-value,是多重假設檢驗校正之后的p-value,即對NES可能存在的假陽性結果的概率估計,因此FDR越小說明富集越顯著;

RANK AT MAX:當ES值**時,對應基因所在排序好的基因列表中所處的位置;

(注:GSEA采用p-value<5%,q-value<25%進行數據過濾)

LEADING EDGE:該處有3個統計值,tags=59%表示核心基因占該基因集中基因總數的百分比;list=21%表示核心基因占所有基因的百分比;

對于分析結果中,我們一般認為|NES|>1,NOM p-val<0.05,FDR q-val<0.25的通路是顯著富集的。

最后還有一個該KEGG基因集下每個基因的詳細統計信息表,RANK IN GENE LIST表示在排序好的基因集中所處的位置;RANK METRIC SCORE是基因排序評分,我們這里是Signal2noise;RUNNING ES是分析過程中動態的ES值;CORE ENRICHMENT是對ES值有主要貢獻的基因,即Leading edge subset,在表中以綠色標記。

轉自R語言分析

文章分類: 科研頭條數據分析

?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 228,333評論 6 531
  • 序言:濱河連續發生了三起死亡事件,死亡現場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發現死者居然都...
    沈念sama閱讀 98,491評論 3 416
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 176,263評論 0 374
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 62,946評論 1 309
  • 正文 為了忘掉前任,我火速辦了婚禮,結果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 71,708評論 6 410
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上,一...
    開封第一講書人閱讀 55,186評論 1 324
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 43,255評論 3 441
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 42,409評論 0 288
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當地人在樹林里發現了一具尸體,經...
    沈念sama閱讀 48,939評論 1 335
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 40,774評論 3 354
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發現自己被綠了。 大學時的朋友給我發了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 42,976評論 1 369
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,518評論 5 359
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質發生泄漏。R本人自食惡果不足惜,卻給世界環境...
    茶點故事閱讀 44,209評論 3 347
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,641評論 0 26
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 35,872評論 1 286
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 51,650評論 3 391
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 47,958評論 2 373

推薦閱讀更多精彩內容

  • 上一篇GSEA可以做什么[http://www.lxweimin.com/p/1cc9f6c20ef6]之后,繼續...
    Y大寬閱讀 38,259評論 2 27
  • Gene set enrichment analysis, 基因集富集分析 1、算法思想 對某基因Knock do...
    DobbyBunny閱讀 7,386評論 0 21
  • GSEA是一種無閾值方法,可根據其差異表達等級或其他分數對所有基因進行分析,無需事先進行基因過濾。當基因組中的所有...
    Mingyan_C閱讀 11,090評論 1 17
  • GSEA縮寫 GSEA的全稱是Gene Set Enrichment Analysis,中文翻譯就是基因集富集分析...
    backup備份閱讀 16,950評論 0 66
  • 網頁上關于GSEA富集分析的教程很多,對自己這個生信小白來說比較簡單通俗易懂的當數生信寶典和生信技能樹的教程了。感...
    阿糖胞苷_SYSU閱讀 17,175評論 0 39