>test = matrix(rnorm(200), 20, 10)
>test[1:10, seq(1, 10, 2)] = test[1:10, seq(1, 10, 2)] + 3
> test[11:20, seq(2, 10, 2)] = test[11:20, seq(2, 10, 2)] + 2
> test[15:20, seq(2, 10, 2)] = test[15:20, seq(2, 10, 2)] + 4
> colnames(test) = paste("Test", 1:10, sep = "")
> rownames(test) = paste("Gene", 1:20, sep = "")
此时test是一个20*10列的矩阵,行名为Gene开头,列名为Test开头
1、对于热力图来说:
> library(pheatmap)
> pheatmap(test,cluster_row = FALSE,cluster_col=FALSE)
图形为:(行列和原矩阵的位置一一对应)
2、对于等高线图来说
> filled.contour(test, color=rainbow,plot.axes=contour(test, add=T))
图形为:(图形刚好为上述热力图逆时针旋转90度,即横轴从左到右依次为Gene1-Gene20,纵轴从上到下依次为Test10-Test1,对应矩阵为原矩阵逆时针旋转90度)