我用以下代码制作了50个最重要基因的热图:
resTC <- results(ddsTC)
topGenes <- head(order(resTC$padj),50)
betas <- coef(ddsTC)
mat <- betas[topGenes, -c(1,2)]
thr <- 5
mat[mat < -thr] <- -thr
mat[mat > thr] <- thr
pheatmap(mat, breaks=seq(from=-thr, to=thr, length=101),
cluster_col=TRUE, cellheight=5, cellwidth = 15, fontsize_row = 5, fontsize_col = 6,
na_col= "black", treeheight_row = 80, treeheight_col = 15, legend_breaks = -5:5)
我遇到的问题是,来自ddsTC分析(默认情况下)的比较名称可能有些混乱。是否可以通过热图重命名它们?