heat

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 3.2.5

library(grid)

setwd("D:/RData/heat/cortext_hip")
data = read.table("data.txt",head=T,sep = "\t")
tmp = data[order(data[,3]),]
tmp <- transform(tmp, gene = factor(gene, levels = unique(gene)))
classes = unique(tmp$class)

mylabs = c("4pcw<Age<8pcw","8pcw<Age<10pcw","10pcw<Age<13pcw","13pcw<Age<16pcw","16pcw<Age<19pcw","19pcw<Age<24pcw","24pcw<Age<38pcw","0M(birth)<Age<6M","6M<Age<12M","1Y<Age<6Y","6Y<Age<12Y","12Y<Age<20Y","20Y<Age<40Y","40Y<Age<60Y","60Y<Age")

ggplot(tmp, aes(period, gene))+
    geom_tile(aes(fill = expression), size = 2)+
    scale_fill_gradientn(colours=c("#052749","#ffffff","#CC1E11"))+
    scale_x_discrete(limits=factor(1:15), labels = mylabs)+
    facet_grid(class ~ ., scales = "free_y")+
    theme(legend.position="top",
          legend.key.width = unit(1,"cm"),
          legend.title = element_blank(),
          panel.background = element_blank(),
          axis.text = element_text(size = 12,color = "black"),
          axis.text.x = element_text(angle=90, hjust=1, vjust=1),
          axis.title = element_blank(),
          axis.ticks = element_blank())