作業要求: 我們統一選擇p<0.05而且abs(logFC)大於一個與眾的基因為顯著差異表達基因集,對這個基因集用R包做KEGG/GO超幾何分布檢驗分析。 然后把表達矩陣和分組信息分別作出cls和gct文件,導入到GSEA軟件分析。 基本任務是完成這個分析。 【1】環境 ...
作業要求: 使用R語言,載入表達矩陣,然后設置好分組信息,統一用DEseq 進行差異分析,當然也可以走走edgeR或者limma的voom流程。 基本任務是得到差異分析結果,進階任務是比較多個差異分析結果的異同點。 安裝DESeq DESeq 對於輸入數據的要求: .DEseq 要求輸入數據是由整數組成的矩陣。 .DESeq 要求矩陣是沒有標准化的。 DESeq 進行差異表達分析 DESeq 分析 ...
2018-07-03 21:47 0 5049 推薦指數:
作業要求: 我們統一選擇p<0.05而且abs(logFC)大於一個與眾的基因為顯著差異表達基因集,對這個基因集用R包做KEGG/GO超幾何分布檢驗分析。 然后把表達矩陣和分組信息分別作出cls和gct文件,導入到GSEA軟件分析。 基本任務是完成這個分析。 【1】環境 ...
引入clusterProfiler與注釋數據 GO(gene ontology)分析 GO,Gene Ontology,是基因功能國際標准分類體系。它旨在建立一個適用於各種物種的,對基因和蛋白質功能進行限定和描述的,並能隨着研究不斷深入而更新的語言詞匯標准。GO分為分子功能 ...
edgeR的介紹 背景 RNA-seq表達譜與生物復制的差異表達分析。 實現一系列基於負二項分布的統計方法,包括經驗貝葉斯估計,精確檢驗,廣義線性模型和准似然檢驗。 與RNA-seq一樣,它可用於產生計數的其他類型基因組數據的差異信號分析,包括ChIP-seq,Bisulfite-seq ...
任務列表 1.在UCSC下載hg19參考基因組; 2.從gencode數據庫下載基因注釋文件,並且用IGV去查看感興趣的基因的結構,比如TP53,KRAS,EGFR等等。 3.截圖幾個基因的IGV可視化結構 4.下載ENSEMBL,NCBI的gtf,也導入IGV看看,截圖 ...
轉錄組差異表達分析小實戰(二) Posted: 八月 14, 2017 Under: Transcriptomics By Kai no Comments 差異基因表達分析 我按照前面的流程轉錄組差異表達分析小實戰(一),將小鼠的4個樣本又重新跑了一遍,從而獲得一個新 ...
Ballgown是分析轉錄組差異表達的R包。 軟件安裝: 運行R, source(“http://bioconductor.org/biocLite.R”) biocLite(“ballgown”) R會自動安裝Ballgown,及相應的依賴包。 Ballgown的輸入文件 ...
轉錄組差異表達分析小實戰(一) Posted: 七月 28, 2017 Under: Transcriptomics By Kai no Comments 讀文獻獲取數據 文獻名稱:AKAP95 regulates splicing through ...
作業要求: 在UCSC下載hg19參考基因組,我博客有詳細說明,從gencode數據庫下載基因注釋文件,並且用IGV去查看你感興趣的基因的結構,比如TP53,KRAS,EGFR等等。 作業,截圖幾個基因的IGV可視化結構!還可以下載ENSEMBL,NCBI的gtf,也導入IGV看看,截圖基因結構 ...