home <- "~/projects/proj_hgc2018/" setwd(home) library(data.table) data1 <- fread("data/data1.txt",drop=1,header=TRUE, data.table = F) rownames(data1) <- fread("data/data1.txt", select=1,header=TRUE, data.table = F)[,1] data1 <- as.matrix(data1) data1<-t(data1) d <- dist(data1) hc <- hclust(d) plot(hc) cutree(hc,2) library(gplots) pdf("output/sample_heatmap.pdf",width = 8,height = 8) heatmap.2(data1,col = bluered(256),trace="none") dev.off()