Heatmap in R / Excel과 Rstudio로 heatmap 그리기

1. Excel data sheet 만들기

Sheet 1 : Group A,B,C,D의 RNAseq data

Sheet2 : Grouping samples

 

2. Rstudio 실행

Rscript ▼

##Heatmap

install.packages("xlsx")

install.packages("pheatmap")

install.packages("realxl")

library("xlsx")

library("pheatmap")

library("readxl")

RNA <- read.xlsx("excelfile.xlsx",1)

rownames(RNA) <- RNA$Gene

RNA_2 <- RNA[,2:5]

list <-read.xlsx("excelfile.xlsx",2)

annotation <- data.frame(Subtype= list [,’Subtypes'],row.names= list$Samples)

heat_colors <- colorRampPalette(c("skyblue","white","Red"))(50)

pheatmap(RNA_2, color= heat_colors, cluster_col=T, cluster_row = F, show_rownames = T,

         show_colnames =F, border_color = NA, annotation= annotation, fontsize = 7, scale="row")

 

3. 결과 확인