R画corplot

相关性系数是判别两个样品之间相似度的一种度量方式，如果有多个样品，你打算看看他们之间的相似性，比较好的一种办法就是画corplot

这里给出来一个R语言的ggplot版本的方法，

install.packages("corrplot")

library(corrplot)

rna<-read.table("merge.xls",header=TRUE)

data<-subset(rna,sum1>0)

rnacor<-cor(data[2:7])

#这里用了col1的概念，其实是corplot官方说明文档上的一段代码，挺好用，大家也可以直接拷过去

col1 <- colorRampPalette(c("#7F0000","red","#FF7F00","yellow","white", 
			"cyan", "#007FFF", "blue","#00007F"))
col2 <- colorRampPalette(c("#67001F", "#B2182B", "#D6604D", "#F4A582", "#FDDBC7",
			"#FFFFFF", "#D1E5F0", "#92C5DE", "#4393C3", "#2166AC", "#053061"))	
col3 <- colorRampPalette(c("red", "white", "blue"))	
col4 <- colorRampPalette(c("#7F0000","red","#FF7F00","yellow","#7FFF7F", 
			"cyan", "#007FFF", "blue","#00007F"))	
wb <- c("white","black")


par(ask = TRUE)

corrplot(rnacor,method="pie",shade.col="red",tl.col="red",tl.offset=0.8,tl.srt=45,order="hclust",addCoef.col="red")