R09

library(pheatmap)

bk <- c(seq(-1,-0.01,by=0.01),seq(0,1,by=0.01))
mycolor <- c(colorRampPalette(colors=c("navy","white"))(length(bk)/2),colorRampPalette(color=c("white","firebrick3"))(length(bk)/2))

r <- read.table('../readcount.matrix',sep="\t",row.names=1,header=T)
rdata <- as.data.frame(r)
rdata <- rdata[apply(r, 1, function(x) sd(x)!=0),]
dpi=2500
p <- pheatmap(rdata,scale = "row", cellwidth = 20, cellheight = 0.5,cluster_cols = F,cluster_rows = T,treeheight_row=0, treeheight_col=0,show_rownames=T, color=mycolor,fontsize = 0.5,fontsize_col=12, breaks=bk,legend=T)

png('heatmap.png',w=12*dpi,h=10*dpi,units = "px",res = dpi,type='cairo')
p
dev.off()

order_row <- p$tree_row$order
#genename <- p$tree_row$labels
#table <- data.frame(order_row,genename)
#neworder <- table[order(table[,1]),]
neworder <- rownames(r[p$tree_row[["lable"]]])[order_row]
write.table(neworder,'neworder.txt',sep="\t")