11.5-9 RNA-seq 數(shù)據(jù)差異基因的得到和可視化學(xué)習(xí)(一)

很倒霉昼接,電腦壞了衣屏。之前未完成的rna-seq數(shù)據(jù)下載和mcscan不能繼續(xù)進(jìn)行了。

最近要在網(wǎng)上直接下載rna-seq的txt數(shù)據(jù)辩棒,所以主要學(xué)習(xí)一下如何挑選出差異基因,之后進(jìn)行可視化膨疏。

http://www.reibang.com/p/3a0e1e3e41d0

1一睁、挑選出差異基因。

This tutorial shows an example of RNA-seq data analysis with DESeq2, followed by KEGG pathway analysis using GAGE.:

http://www.gettinggeneticsdone.com/2015/12/tutorial-rna-seq-differential.html

2、可視化

由于我們的下載數(shù)據(jù)是rpkm。所以我們直接log2 transform髓棋。

http://www.reibang.com/p/807cf4a969fb

1、MA plot

在DESeq2中,函數(shù)plotMA顯示了一個給定變量對DESeqDataSet中所有樣本的歸一化計(jì)數(shù)均值的log2倍變化豺旬。如果調(diào)整后的p值小于0.1扛拨,點(diǎn)會被涂成紅色。從窗口掉下來的點(diǎn)被畫成向上或向下的開放三角形渴频。

wikipedia 公式

https://en.wikipedia.org/wiki/MA_plot

Several Bioconductor packages, for the R software, provide the facility for creating MA plots. These include affy (ma.plot, mva.pairs), limma (plotMA), marray (maPlot), and edgeR(maPlot)

2场钉、Plot counts

DESeq2提供了一個plotCounts()函數(shù)來查看某一個感興趣的gene在組間的差別批钠。

?plotCounts()來看如何使用

3坡氯、熱圖

我們主要根據(jù)http://www.reibang.com/p/398115d2d2e8來學(xué)習(xí)

heatmap():用于繪制簡單熱圖的函數(shù)

heatmap.2():繪制增強(qiáng)熱圖的函數(shù)

d3heatmap:用于繪制交互式熱圖的R包

ComplexHeatmap:用于繪制、注釋和排列復(fù)雜熱圖的R&bioconductor包(非常適用于基因組數(shù)據(jù)分析)

1、首先,最簡單的熱圖

將id作為列名

使用as.matrix不可以將list轉(zhuǎn)化為矩陣,所以:

轉(zhuǎn)化為數(shù)值矩陣

usetest<-matrix(unlist(test),6,14)

?heatmap(usetest,scale="column")得到

這里的數(shù)字是可以用rowname和colname語句換成想要的。

簡單換顏色唯鸭,可以加上一句 ?

col <- colorRampPalette(c("red", "white", "blue"))(256)? 在heatmap前

2缭付、增強(qiáng)版

gplots中的heatmap.2

heatmap.2(usetest,scale="column",col=bluered(20),trace = "none",density.info = "none")

heatmap


對于文本大小的調(diào)節(jié):

https://www.plob.org/article/10045.html/comment-page-1/

heatmap.2(x, srtCol=0, adjCol = c(0.5,1) )

heatmap.2(x, srtCol=45, adjCol = c(1,1) )舅巷。。媚值。嫡良。。等

這篇文章還講了如何運(yùn)用多種聚類的方法京闰。

今天又畫了一次熱圖

baifen<-read.csv("baifen.csv",header=T)

rownames(baifen)<-baifen[,1]

baifen<-baifen[,-1]

mabaifen<-as.matrix(scale(baifen))

heatmap.2(mabaifen, keysize=1.5,key.title = NA,col=colorRampPalette(c("green","white","red")),scale="none",trace = "none",dendrogram="none",Colv = FALSE,cexCol=1,cexRow=0.8)

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市萌丈,隨后出現(xiàn)的幾起案子辆雾,更是在濱河造成了極大的恐慌坛梁,老刑警劉巖政鼠,帶你破解...
    沈念sama閱讀 217,185評論 6 503
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件遏佣,死亡現(xiàn)場離奇詭異,居然都是意外死亡敞曹,警方通過查閱死者的電腦和手機(jī)抓歼,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,652評論 3 393
  • 文/潘曉璐 我一進(jìn)店門湃窍,熙熙樓的掌柜王于貴愁眉苦臉地迎上來役衡,“玉大人,你說我怎么就攤上這事棵介《质觯” “怎么了?”我有些...
    開封第一講書人閱讀 163,524評論 0 353
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經(jīng)常有香客問我痴腌,道長腕侄,這世上最難降的妖魔是什么毁枯? 我笑而不...
    開封第一講書人閱讀 58,339評論 1 293
  • 正文 為了忘掉前任蹂窖,我火速辦了婚禮,結(jié)果婚禮上恩敌,老公的妹妹穿的比我還像新娘瞬测。我一直安慰自己,他們只是感情好纠炮,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,387評論 6 391
  • 文/花漫 我一把揭開白布月趟。 她就那樣靜靜地躺著,像睡著了一般恢口。 火紅的嫁衣襯著肌膚如雪孝宗。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 51,287評論 1 301
  • 那天耕肩,我揣著相機(jī)與錄音因妇,去河邊找鬼。 笑死猿诸,一個胖子當(dāng)著我的面吹牛婚被,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播梳虽,決...
    沈念sama閱讀 40,130評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼址芯,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起谷炸,我...
    開封第一講書人閱讀 38,985評論 0 275
  • 序言:老撾萬榮一對情侶失蹤删顶,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后淑廊,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 45,420評論 1 313
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡特咆,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,617評論 3 334
  • 正文 我和宋清朗相戀三年季惩,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片腻格。...
    茶點(diǎn)故事閱讀 39,779評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡画拾,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出菜职,到底是詐尸還是另有隱情青抛,我是刑警寧澤,帶...
    沈念sama閱讀 35,477評論 5 345
  • 正文 年R本政府宣布酬核,位于F島的核電站蜜另,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏嫡意。R本人自食惡果不足惜举瑰,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,088評論 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望蔬螟。 院中可真熱鬧此迅,春花似錦、人聲如沸旧巾。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,716評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽鲁猩。三九已至坎怪,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間廓握,已是汗流浹背芋忿。 一陣腳步聲響...
    開封第一講書人閱讀 32,857評論 1 269
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留疾棵,地道東北人戈钢。 一個月前我還...
    沈念sama閱讀 47,876評論 2 370
  • 正文 我出身青樓,卻偏偏與公主長得像是尔,于是被迫代替她去往敵國和親殉了。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,700評論 2 354

推薦閱讀更多精彩內(nèi)容