作業要求: 我們統一選擇p<0.05而且abs(logFC)大於一個與眾的基因為顯著差異表達基因集,對這個基因集用R包做KEGG/GO超幾何分布檢驗分析。 然后把表達矩陣和分組信息分別作出cls和gct文件,導入到GSEA軟件分析。 基本任務是完成這個分析。 【1】環境 ...
引入clusterProfiler與注釋數據 GO gene ontology 分析 GO,Gene Ontology,是基因功能國際標准分類體系。它旨在建立一個適用於各種物種的,對基因和蛋白質功能進行限定和描述的,並能隨着研究不斷深入而更新的語言詞匯標准。GO分為分子功能 Molecular Function 生物過程 Biological Process 和細胞組成 Cellular Comp ...
2020-08-10 12:28 0 580 推薦指數:
作業要求: 我們統一選擇p<0.05而且abs(logFC)大於一個與眾的基因為顯著差異表達基因集,對這個基因集用R包做KEGG/GO超幾何分布檢驗分析。 然后把表達矩陣和分組信息分別作出cls和gct文件,導入到GSEA軟件分析。 基本任務是完成這個分析。 【1】環境 ...
作業要求: 使用R語言,載入表達矩陣,然后設置好分組信息,統一用DEseq2進行差異分析,當然也可以走走edgeR或者limma的voom流程。 基本任務是得到差異分析結果,進階任務是比較多個差異分析結果的異同點。 【1】安裝DESeq2 DESeq2對於輸入數據 ...
參考基因組下載 基因組各版本的對應關系http://www.bio-info-trainee.com/1469.html GRCh36 (hg18): ENSEMBL release_52. GRCh37 (hg19): ENSEMBL release_59/61/64/68/69/75. ...
任務列表 1.在UCSC下載hg19參考基因組; 2.從gencode數據庫下載基因注釋文件,並且用IGV去查看感興趣的基因的結構,比如TP53,KRAS,EGFR等等。 3.截圖幾個基因的IGV可視化結構 4.下載ENSEMBL,NCBI的gtf,也導入IGV看看,截圖 ...
作業要求: 在UCSC下載hg19參考基因組,我博客有詳細說明,從gencode數據庫下載基因注釋文件,並且用IGV去查看你感興趣的基因的結構,比如TP53,KRAS,EGFR等等。 作業,截圖幾個基因的IGV可視化結構!還可以下載ENSEMBL,NCBI的gtf,也導入IGV看看,截圖基因結構 ...
參考:http://www.biotrainee.com/thread-558-1-1.html http://bioconductor.org/packages/3.7/bioc/ http:/ ...
edgeR的介紹 背景 RNA-seq表達譜與生物復制的差異表達分析。 實現一系列基於負二項分布的統計方法,包括經驗貝葉斯估計,精確檢驗,廣義線性模型和准似然檢驗。 與RNA-seq一樣,它可用於產生計數的其他類型基因組數據的差異信號分析,包括ChIP-seq,Bisulfite-seq ...
基因組注釋主要包括四個研究方向:重復序列的識別;非編碼RNA的預測;基因結構預測和基因功能注釋。我們將分別對這四個領域進行闡述。 1 重復序列的識別。 1.1 重復序列的研究背景和意義:重復序列可分為串聯重復序列(Tendam repeat)和散在重復序列 ...