轉(zhuǎn)載自http://www.sohu.com/a/200219694_278730

接觸這么多組學(xué)數(shù)據(jù)，大家肯定做過各種相關(guān)性分析。大到幾個轉(zhuǎn)錄組樣本的整體相關(guān)性分析，小到挑選了一些候選基因看它們在不同樣本中的表達(dá)模式相關(guān)性。當(dāng)我們做完了相關(guān)性分析的時候，常常會用下面這種圖形來展示：

465.jpeg

那么問題來了，有沒有別的方法，可以讓我的相關(guān)性分析展現(xiàn)方式更高大上呢？

請出今天的主角，corrplot程序包！

老規(guī)矩，第一步是安裝并調(diào)用corrplot：

install.packages("corrplot")
library("corrplot")

第二步，整理數(shù)據(jù)。今天就拿一些女明星的身高、體重、年齡、微博粉絲與女神指數(shù)的數(shù)據(jù)來分析一下（純屬杜撰）

87954.png

第三步，我們讀取數(shù)據(jù)，并做相關(guān)性分析。

a<-read.table(file=file.choose(),row.names=1,header=T,sep="t")
b<-cor(a)

第四步，非常簡單，可以直接作圖了。

corrplot(b) #這里默認(rèn)是圓形顯示

564123.png

還可以用餅圖（pie），顏色（color）等等來顯示。

corrplot(b,method="pie")

897465.png

corrplot(b,method="color",addCoef.col="grey") 
#用顏色顯示，同時顯示相關(guān)系數(shù)，是不是跟開頭綠綠的圖一樣啦。

56123.png

第五步，是時候展示高端操作了！我們把顏色對調(diào)一下，正相關(guān)為紅色，負(fù)相關(guān)為海軍藍(lán)色，同時把右上角用圓形展示。

col=colorRampPalette(c("navy", "white", "firebrick3")) #設(shè)置顏色
corrplot(b,type="upper",col=col(10),tl.pos="d") #tl.pos="d"即不顯示周邊各列名字

465123.png