@果蠅飼養(yǎng)員的生信筆記 哦哦 好嘞 學(xué)習(xí)老師解答
使用clusterProfiler進(jìn)行富集分析clusterProfiler是業(yè)界大神Y叔寫的R包偷拔,看過了它的繪圖就會(huì)使人沉迷于其顏值無法自拔。本文很多內(nèi)容直接翻譯自官方文檔亏钩,加上了我的測試數(shù)據(jù)(還是以果蠅為例)莲绰。 一、...
@wo_monic 好的 我試試看 謝謝
構(gòu)建自己的R包--KEGG.db參考地址在使用clusterProfiler時(shí)姑丑,KEGG注釋包蛤签,可以自定義。畢竟許多物種的數(shù)據(jù)庫是不完善的≌ぐВ現(xiàn)在就自定義一個(gè)KEGG.db的包震肮,用于KEGG注釋。 簡單快捷方...
老師 寫得好好留拾!不過我疑惑就是:為什么一開始獲得geneList之后就沒有看到這個(gè)使用這個(gè)變量了 geneList是后面用到的mydata嘛戳晌? 另外GSEA不是用所有基因而非差異基因嘛 為什么DEgenelist要設(shè)置logFC>2呀 謝謝老師
使用clusterProfiler進(jìn)行富集分析clusterProfiler是業(yè)界大神Y叔寫的R包,看過了它的繪圖就會(huì)使人沉迷于其顏值無法自拔痴柔。本文很多內(nèi)容直接翻譯自官方文檔沦偎,加上了我的測試數(shù)據(jù)(還是以果蠅為例)。 一咳蔚、...
你好 為什么我在第二步就報(bào)錯(cuò)了 createKEGG::create_kegg_db("mmu")后報(bào)錯(cuò):不能給維度小于二的對象設(shè)‘colnames’
構(gòu)建自己的R包--KEGG.db參考地址在使用clusterProfiler時(shí)豪嚎,KEGG注釋包,可以自定義谈火。畢竟許多物種的數(shù)據(jù)庫是不完善的〕扪現(xiàn)在就自定義一個(gè)KEGG.db的包,用于KEGG注釋堆巧。 簡單快捷方...
批次效應(yīng)除去-GEO案例分析 最近分析來自不同平臺(tái)的GEO數(shù)據(jù)妄荔,GEO分析推薦學(xué)習(xí)Jimmy大神的B站視頻泼菌,也可以學(xué)習(xí)我之前的教程。limma和sva兩種方法進(jìn)行批次效應(yīng)去除...
首先是學(xué)習(xí)優(yōu)秀的帖子 最近處理數(shù)據(jù)發(fā)現(xiàn)同一個(gè)平臺(tái)的RNA-seq或者芯片數(shù)據(jù)啦租,不同的樣本之間有時(shí)候竟然存在很大的差異哗伯,也就是所謂的批次效應(yīng)。因此學(xué)習(xí)了一下篷角,并推薦一下優(yōu)秀的推...
本文講的是aglient芯片原始數(shù)據(jù)的處理,參考資料是limma 的userguide文檔嘉蕾。GEO數(shù)據(jù)庫下載的表達(dá)矩陣不符合預(yù)期贺奠,比如是空的,或者是有負(fù)值的错忱,那我們就處理一下...
通常情況下儡率,通過以下幾種指標(biāo)來對模型進(jìn)行評價(jià)。1)區(qū)分度:采用指標(biāo)C-index和ROC曲線來評價(jià)區(qū)分度以清,一般文章都是二選一儿普。 C-index即一致性指數(shù)(index of ...
說來慚愧,跟師姐聊天的時(shí)候掷倔,我按照傳統(tǒng)的生信思路把臨床特征做風(fēng)險(xiǎn)打分眉孩,然后定義高低風(fēng)險(xiǎn)打分,接著就是標(biāo)準(zhǔn)的生存圖勒葱、森林圖了浪汪。在做到森林圖的時(shí)候,忽然被問到凛虽,單因素與多因素的區(qū)...
實(shí)例:如何拿到KEGG數(shù)據(jù)庫中多巴胺通路相關(guān)的基因集 一吟宦、確定目標(biāo)通路 打開KEGG選擇pathway,在搜索框前輸入物種涩维,框內(nèi)填入關(guān)鍵詞。 篩選結(jié)果顯示僅有hsa04728...
一袁波、TCGA轉(zhuǎn)錄組數(shù)據(jù)下載及整理 ①進(jìn)去先看cart是否為空瓦阐,沒有要清空。 ②下載轉(zhuǎn)錄組數(shù)據(jù) 基因表達(dá)量包括LncRNA和mRNA篷牌,用矯正后的FPKM下載3個(gè)文件 ③下載數(shù)據(jù)...