需求 和這一篇的需求一致，不同公司的芯片原始數(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.需求 如圖，熱圖自帶聚類樹褐捻，將樣本(列)分成了三個(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)，我知道一寫統(tǒng)計(jì)學(xué)就會(huì)閱讀量暴跌爱致，沒(méi)人看就當(dāng)寫給自己看...
-
上下調(diào)基因分別富集烤送，也能一起畫圖雙向富集分析圖 1.需求 上下調(diào)基因分別富集分析，把結(jié)果放在一張圖上: 2.示例數(shù)據(jù)和r包準(zhǔn)備 模擬的是kegg的富集結(jié)果糠悯，編了一個(gè)數(shù)據(jù)帮坚。p是p值。畫圖的時(shí)候用-log10(...
-
即使你不normal互艾，平均值也normal--神奇的中心極限定理0.寫在前面 統(tǒng)計(jì)學(xué)知識(shí)刷起來(lái)试和。其實(shí)這又是一個(gè)系列，晦澀難懂沒(méi)人看的那種纫普。但這樣的知識(shí)確實(shí)非常有價(jià)值阅悍，能夠帶給我以及認(rèn)真學(xué)習(xí)的讀者們實(shí)打?qū)嵉倪M(jìn)步，所以不管閱讀量高低昨稼，我都寫节视。...
-
披著ggplot皮的pheatmap,深夜激動(dòng)更新我的包0.寫在前面 從我學(xué)生信開(kāi)始，畫熱圖用的就是pheatmap悦昵，有時(shí)復(fù)雜化的需求會(huì)用到ComplexHeatmap肴茄。作為R語(yǔ)言畫圖的絕對(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è)圖逻谦，畫它。 其實(shí)就是個(gè)點(diǎn)圖加棍陪蜻，棒棒糖圖邦马。找到了示例數(shù)據(jù)和代碼： https://www.r-graph-gallery.com/301-cu...
-
-
單細(xì)胞交響樂(lè)29-實(shí)戰(zhàn)十二 10X 小鼠嵌合體胚胎劉小澤寫于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)度医寿。...
-
