#データの読み込み #TN_mRNA_expression.tsv library(data.table) expr <- as.matrix(fread("data/TN_mRNA_expression.tsv", header = T, sep = "\t", drop = 1)) rownames(expr) <- as.vector(as.matrix(fread("data/TN_mRNA_expression.tsv", select = 1))) #expr <- read.table("data/TN_mRNA_expression.tsv",header = T, sep = "\t") #遺伝子のフィルタリング #遺伝子ごとに分散を計算 v <- apply(expr,1,var) #分散の大きい上位100遺伝子の発現データを取り出す idx <- order(v, decreasing = TRUE)[1:100] topvar <- expr[idx, ] #ヒートマップの描画 pdf("data/test_heatmap01.pdf") heatmap.2(topvar,trace = "none", scale = "row", col= bluered(256)) dev.off() pdf("data/test_heatmap02.pdf", width = 10, height = 10) cols <- colorpanel(n = 256, low = "white", high = "red") heatmap.2(topvar,trace = "none", scale = "none", col= cols) dev.off()