需求 和這一篇的需求一致耿焊，不同公司的芯片原始數(shù)據(jù)處理方式也不同殉挽，最終都是為了得到表達(dá)矩陣峻凫，今天介紹的是illumina beadchip 芯片秋泄，用神奇的limma搞定它香伴。找不...
IP屬地：廣東
-
illumina beadchip 芯片原始數(shù)據(jù)處理 -
快捷查找KEGG里的通路和基因
需求 1.快捷查找ID對(duì)應(yīng)的description揽趾，知道通路對(duì)應(yīng)的編號(hào)是多少小作。2.找出某一個(gè)/幾個(gè)通路里的全部基因亭姥，用來(lái)做單獨(dú)的下游分析。 如果是要做KEGG的富集分析顾稀，cl...
-
讓TCGA表達(dá)矩陣和臨床信息方便的匹配
0.背景 關(guān)于樣本id和分組信息：TCGA的樣本id里藏著分組信息[https://mp.weixin.qq.com/s/UohDh5fHkTPlwNp43ivWHg] TC...
-
簡(jiǎn)單高效滴處理R語(yǔ)言里的小數(shù)
舉個(gè)栗子 這是一組隨機(jī)數(shù)达罗，假如需要把這樣的一組數(shù)字放在圖里，是需要整齊統(tǒng)一的顯示出來(lái)静秆×溉啵基本上有兩種思路，一是保留相同的位數(shù)抚笔，二是用相同格式的科學(xué)計(jì)數(shù)法表示扶认。 方法一：統(tǒng)一保留...
-
哦，走開(kāi)殊橙，你這個(gè)可惡的離群值1.離群值是啥子[http://222.70.95.243:6893/files/tmp/%E7%A6%BB%E7%BE%A4%E5%80%BC%EF%BC%8C%E6%88...
-
按照熱圖聚類結(jié)果來(lái)給樣本分組1.需求 如圖辐宾，熱圖自帶聚類樹(shù)狱从，將樣本(列)分成了三個(gè)簇。想要按照聚類的結(jié)果叠纹，獲取到 的信息矫夯，并將這個(gè)信息添加到熱圖上。 2.作圖的數(shù)據(jù) 輸入數(shù)據(jù)是個(gè)表達(dá)矩陣吊洼，行為基因训貌，列為...
-
emapplot-原來(lái)能用的函數(shù)，現(xiàn)在報(bào)錯(cuò)了腫么辦emapplot是富集分析可視化的一種冒窍，展示有共同基因的通路递沪，詳見(jiàn)：http://www.reibang.com/p/c45cc2e3890f[https://www.ji...
-
讓我們蕩起雙槳，重新認(rèn)識(shí)一下p值在正文之前先分享一點(diǎn)生活故事： 上周天傍晚综液，我閃著腰了款慨，躺了一晚上不見(jiàn)好，發(fā)朋友圈尋醫(yī)問(wèn)藥了： 因?yàn)樽龅氖巧欧窒砻ǎ遗笥讶镝t(yī)學(xué)生和醫(yī)生非常多檩奠，我想看他們?cè)趺凑f(shuō)。我收到了幾...
-
相關(guān)性系數(shù)和它的朋友們本文內(nèi)容是基于statquest-協(xié)方差與相關(guān)系數(shù)的學(xué)習(xí)筆記附帽，介紹協(xié)方差埠戳、相關(guān)系數(shù)和R^1三個(gè)統(tǒng)計(jì)量，并添加了R語(yǔ)言實(shí)現(xiàn)蕉扮，我知道一寫(xiě)統(tǒng)計(jì)學(xué)就會(huì)閱讀量暴跌整胃，沒(méi)人看就當(dāng)寫(xiě)給自己看...
-
上下調(diào)基因分別富集，也能一起畫(huà)圖雙向富集分析圖 1.需求 上下調(diào)基因分別富集分析喳钟，把結(jié)果放在一張圖上: 2.示例數(shù)據(jù)和r包準(zhǔn)備 模擬的是kegg的富集結(jié)果屁使，編了一個(gè)數(shù)據(jù)。p是p值奔则。畫(huà)圖的時(shí)候用-log10(...
-
即使你不normal蛮寂，平均值也normal--神奇的中心極限定理0.寫(xiě)在前面 統(tǒng)計(jì)學(xué)知識(shí)刷起來(lái)。其實(shí)這又是一個(gè)系列易茬，晦澀難懂沒(méi)人看的那種酬蹋。但這樣的知識(shí)確實(shí)非常有價(jià)值，能夠帶給我以及認(rèn)真學(xué)習(xí)的讀者們實(shí)打?qū)嵉倪M(jìn)步疾呻，所以不管閱讀量高低除嘹，我都寫(xiě)写半。...
-
披著ggplot皮的pheatmap,深夜激動(dòng)更新我的包0.寫(xiě)在前面 從我學(xué)生信開(kāi)始岸蜗，畫(huà)熱圖用的就是pheatmap，有時(shí)復(fù)雜化的需求會(huì)用到ComplexHeatmap叠蝇。作為R語(yǔ)言畫(huà)圖的絕對(duì)霸主璃岳，ggplot2的熱圖卻無(wú)法在熱圖的...
-
KEGG年缎、GO富集分析與GSEA并不是二選一最廣為人知的富集分析做法是把上調(diào)、下調(diào)基因分別或者合并铃慷，拿來(lái)做GO和KEGG富集分析单芜。經(jīng)常有一些數(shù)據(jù)集，拿差異基因做得不到結(jié)果犁柜，那是因?yàn)榇_實(shí)富集不到任何通路洲鸠，是正常的。不妨試...
-
你知道韋恩圖管誰(shuí)叫大哥不馋缅？臨近國(guó)慶節(jié)啦扒腕，花花9月19號(hào)上完了9月的線上課，飛回自己家里當(dāng)廢物啦（不要問(wèn)我為什么不帶豆豆萤悴，回家當(dāng)廢物怎么能帶豆豆呢瘾腰，帶了豆豆還怎么當(dāng)廢物呢）何況豆豆是一個(gè)要按時(shí)按點(diǎn)上班的...
-
模型不如意十之八九，可拿出手無(wú)二三今天又來(lái)開(kāi)放性討論啦覆履。最近做了一個(gè)TCGA的數(shù)據(jù)分析蹋盆，走到最后多因素cox模型構(gòu)建這一步了，選擇了三個(gè)基因硝全，對(duì)數(shù)據(jù)集劃分了高低風(fēng)險(xiǎn)栖雾，C-index值和timeROC計(jì)算的AU...
-
棒棒糖圖，好吃也好用0. 需求 群里看到學(xué)員求助一個(gè)圖伟众，畫(huà)它岩灭。 其實(shí)就是個(gè)點(diǎn)圖加棍，棒棒糖圖赂鲤。找到了示例數(shù)據(jù)和代碼： https://www.r-graph-gallery.com/301-cu...
-
-
單細(xì)胞交響樂(lè)29-實(shí)戰(zhàn)十二 10X 小鼠嵌合體胚胎劉小澤寫(xiě)于2020.7.21為何取名叫“交響樂(lè)”？因?yàn)閱渭?xì)胞分析就像一個(gè)大樂(lè)團(tuán)泡孩，需要各個(gè)流程的協(xié)同配合單細(xì)胞交響樂(lè)1-常用的數(shù)據(jù)結(jié)構(gòu)SingleCellExperiment單...
-
弦圖车摄，管他中不中用，反正高大上0.準(zhǔn)備輸入數(shù)據(jù) 可以是矩陣 行是一邊仑鸥，列是一邊吮播，矩陣正文的數(shù)字表示連接的強(qiáng)度。 可以是tidy數(shù)據(jù)框 一行是一個(gè)連接關(guān)系眼俊，第一列和第二列是要連接的兩邊意狠，第三列是連接的強(qiáng)度。...
-
