genes<-read.csv("E:/Rdata/data/zhangjia.csv")
library(RColorBrewer)
#rownames(genes)<-genes[1,]
data_matrix<-data.matrix(genes)
pal=brewer.pal(6,"Greys")
breaks<-seq(0,1,0.2)
layout(matrix(data=c(1,2),nrow=1,ncol=2), widths=c(8,1),
heights=c(1,1))
par(mar=c(3,7,12,2),oma=c(0.2,0.2,0.2,0.2), mex=0.5)
image(x=1:nrow(data_matrix),y=1:ncol(data_matrix),
z=data_matrix,xlab="",ylab="",breaks=breaks,
col=pal[1:(length(breaks)-1)],axes=F)
#text(x=1:nrow(data_matrix)+0.75, y=par("usr")[4]+1.75,
#srt = 45,adj=1, labels=colnames(data_matrix),
#xpd=T)
#axis(2, at=1:ncol(data_matrix),labels=colnames(data_matrix),col="white",las=1)
abline(h=c(1:ncol(data_matrix))+0.5, v=c(1:nrow(data_matrix))+0.5,
col="white",lwd=2,xpd=F)
title("Correlation between genes",line=8,adj=0)
breaks2<-breaks[-length(breaks)]
par(mar = c(22,1,22,13))
image(x=0.5,y=0:length(breaks2),z=t(matrix(breaks2))*1.001,
col=pal[1:(length(breaks)-1)],axes=F,
breaks=breaks,xlab="",ylab="",
xaxt="n")
axis(4,at=0:(length(breaks2)),labels=breaks,col="white",las=1)
abline(h=c(1:length(breaks2)),col="white",lwd=2,xpd=F)
这是我做热图的代码，你参考下
