pheatmap cytokine_annotation

Question

我试图在pheatmap中开发cytokine_annotation并获取错误消息。

Seq.int中的错误(rx1L，rx2L，length.out = nb)：'from‘必须是有限的

R版本3.3.3 pheatmap_1.0.8

可复制的例子：

#Using cytokine annotations
M<-matrix(rnorm(8*20),ncol=8)
row_annotation<-data.frame(A=gl(4,nrow(M)/4),B=gl(4,nrow(M)/4))
eg<-expand.grid(factor(c(0,1)),factor(c(0,1)),factor(c(0,1)))
colnames(eg)<-c("IFNg","TNFa","IL2")
rownames(eg)<-apply(eg,1,function(x)paste0(x,collapse=""))
rownames(M)<-1:nrow(M)
colnames(M)<-rownames(eg)
cytokine_annotation=eg
pheatmap(M,annotation=annotation,row_annotation=row_annotation,annotation_legend=TRUE,row_annotation_legend=TRUE,cluster_rows=FALSE,cytokine_annotation=cytokine_annotation,cluster_cols=FALSE)

Answer 1

在我的R3.3.3中使用pheatmap_1.0.8，下面的代码可以工作：

set.seed(1)
M <- matrix(rnorm(8*20),ncol=8)
row_annotation <- data.frame(A=gl(4,nrow(M)/4),B=gl(4,nrow(M)/4))
eg <- expand.grid(factor(c(0,1)),factor(c(0,1)),factor(c(0,1)))
colnames(eg) <- c("IFNg","TNFa","IL2")
rownames(eg) <- apply(eg,1,function(x)paste0(x,collapse=""))
rownames(M) <- 1:nrow(M)
colnames(M) <- rownames(eg)
cytokine_annotation <- eg
library(pheatmap)
pheatmap(M, annotation=cytokine_annotation, row_annotation=row_annotation,
            annotation_legend=TRUE, row_annotation_legend=TRUE, cluster_rows=FALSE,
            cytokine_annotation=cytokine_annotation, cluster_cols=FALSE)

​

​

Answer 2

我也有这个问题，发现简单地重新启动R解决了这个问题。