使用图像（）绘图并注释行

Question

我有一个很大的数据框（下面有一个类似的小例子），我想使用image（）显示矩阵中的值。 首先，它绘制了从右边第一列开始并且从右向左的列，但是，我希望它在左边的第一列进行绘图，然后从左向右进行绘制。

其次，我想添加一个'颜色条'的情节表明行1和行2具有相同的$ taxid，行3和行4具有相同的$ taxid。类似于热图中的“ColSideColors”。

有什么建议吗？

谢谢

df<-data.frame(taxid=c("A","A","B","B"),c2=c(1,2,3,4),c3=c(3,4,5,6)) 
nr<-dim(df)[1] 
nc<-dim(df)[2] 
z<-as.matrix(log(df[,2:3])) 
image(x=1:nc, y=1:nr, z=t(z), axes=FALSE, frame.plot=TRUE) 

好了，我现在已经决定尝试，并为colorstip一个“瘦” barplot。但是，它不是'叠加'barplot（虽然高度是一个向量，并且旁边= FALSE。因此，当我使用布局将它们组合在一起时，它抱怨没有足够大的边距。

任何关于堆叠barplot，或其他更复杂的解决方案 感谢

layout(matrix(c(1,2),1,2,byrow=FALSE),heights=c(1,1),widths=c(1,9)) 
barplot(as.vector(table(df$taxid)),width=0.1,beside=FALSE,col=rainbow(length(table(df$taxid)))) 
image(x=1:nc, y=1:nr, z=t(z), axes=FALSE, frame.plot=TRUE) 

Answer 1

可以使用rasterImage功能色带添加到现有的情节：

tmp <- as.raster( rev(c('#00ff00','#0000ff')[df$taxid])) 
par(xpd=NA) 
rasterImage(tmp, 3, 0.5, 3.1, 4.5, interpolate=FALSE) 

如果没有rasterImage选项，那么你可以覆盖2张图片。在下面，我将2列NA添加到原始矩阵并扩展x=参数，以便在右侧添加一个小空白区域，然后创建一个新的矩阵，并在第4列（所有其他NA）中使用taxid信息并使用它使用相同的x=，但使用add=TRUE，以便它在旧的顶部创建新的图像图（但只在先前的空白区域添加颜色，并且不会重叠绘制感兴趣的图像，因为这些值为NA）。

z2 <- cbind(z, NA, NA) 
image(x=c(1:3,3.05, 3.1), y=1:nr, z=t(z2), axes=FALSE) 

z3 <- cbind(NA,NA,NA, as.numeric(df$taxid)) 
image(x=c(1:3,3.05, 3.1), y=1:nr, z=t(z3), add=TRUE, col=c('#00ff00','#0000ff'))