240 發(fā)簡信
IP屬地:上海
  • 120
    GO富集分析之弦圖

    本文使用 GOplot 包繪制弦圖 首先要準(zhǔn)備兩個(gè)文件俊戳,一個(gè)是GO富集分析的結(jié)果,一個(gè)是基因差異表達(dá)分析的結(jié)果瞬逊。一般我們都是用差異表達(dá)的基因進(jìn)行GO富集分析碱呼,所以這兩個(gè)文件還...

  • @果蠅飼養(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注釋堆巧。 簡單快捷方...

  • 120
    來自不同平臺(tái)GEO數(shù)據(jù)批次效應(yīng)除去

    批次效應(yīng)除去-GEO案例分析 最近分析來自不同平臺(tái)的GEO數(shù)據(jù)妄荔,GEO分析推薦學(xué)習(xí)Jimmy大神的B站視頻泼菌,也可以學(xué)習(xí)我之前的教程。limma和sva兩種方法進(jìn)行批次效應(yīng)去除...

  • 批次效應(yīng)學(xué)習(xí)

    首先是學(xué)習(xí)優(yōu)秀的帖子 最近處理數(shù)據(jù)發(fā)現(xiàn)同一個(gè)平臺(tái)的RNA-seq或者芯片數(shù)據(jù)啦租,不同的樣本之間有時(shí)候竟然存在很大的差異哗伯,也就是所謂的批次效應(yīng)。因此學(xué)習(xí)了一下篷角,并推薦一下優(yōu)秀的推...

  • 120
    aglient芯片原始數(shù)據(jù)處理

    本文講的是aglient芯片原始數(shù)據(jù)的處理,參考資料是limma 的userguide文檔嘉蕾。GEO數(shù)據(jù)庫下載的表達(dá)矩陣不符合預(yù)期贺奠,比如是空的,或者是有負(fù)值的错忱,那我們就處理一下...

  • 120
    TCGA+biomarker——ROC曲線

    通常情況下儡率,通過以下幾種指標(biāo)來對模型進(jìn)行評價(jià)。1)區(qū)分度:采用指標(biāo)C-index和ROC曲線來評價(jià)區(qū)分度以清,一般文章都是二選一儿普。 C-index即一致性指數(shù)(index of ...

  • 單因素與多因素cox分析區(qū)別

    說來慚愧,跟師姐聊天的時(shí)候掷倔,我按照傳統(tǒng)的生信思路把臨床特征做風(fēng)險(xiǎn)打分眉孩,然后定義高低風(fēng)險(xiǎn)打分,接著就是標(biāo)準(zhǔn)的生存圖勒葱、森林圖了浪汪。在做到森林圖的時(shí)候,忽然被問到凛虽,單因素與多因素的區(qū)...

  • 120
    R-下載某一條通路的所有基因名字(KEGG)

    實(shí)例:如何拿到KEGG數(shù)據(jù)庫中多巴胺通路相關(guān)的基因集 一吟宦、確定目標(biāo)通路 打開KEGG選擇pathway,在搜索框前輸入物種涩维,框內(nèi)填入關(guān)鍵詞。 篩選結(jié)果顯示僅有hsa04728...

  • 120
    TCGA轉(zhuǎn)錄組數(shù)據(jù)及臨床數(shù)據(jù)下載及整理

    一袁波、TCGA轉(zhuǎn)錄組數(shù)據(jù)下載及整理 ①進(jìn)去先看cart是否為空瓦阐,沒有要清空。 ②下載轉(zhuǎn)錄組數(shù)據(jù) 基因表達(dá)量包括LncRNA和mRNA篷牌,用矯正后的FPKM下載3個(gè)文件 ③下載數(shù)據(jù)...

亚洲A日韩AV无卡,小受高潮白浆痉挛av免费观看,成人AV无码久久久久不卡网站,国产AV日韩精品