有没有办法在使用R pheatmap包时关闭群集?

时间:2013-08-12 09:12:10

标签: r heatmap hierarchical-clustering

我创建了一个矩阵,现在我想使用pheatmap绘制热图,同时保留矩阵行的顺序。我想关闭群集。目前,当我运行以下脚本时,pheatamp正在对行进行聚类:

    tissuedata<-read.table('OrderedByGeneCounts.csv', header=TRUE, sep=',')
    tissuedata2<-data.matrix(tissuedata)
    tissuedata3<-log10(tissuedata2)
    is.na(tissuedata3) <- do.call(cbind,lapply(tissuedata3, is.infinite))
    tissuedata3[is.na(tissuedata3)] <- 0
    tissue_exons<-scan(file="OrderedExonsList.txt", what="character", strip.white=TRUE)
    rownames(tissuedata3)<-tissue_exons

    pdf(file="pheatmapENSG00000026025.pdf", paper="a4", height=15)

    pheatmap(tissuedata3,
    color = colorRampPalette(rev(c("#D73027", "#FC8D59", "#FEE090", "#FFFFBF", "#E0F3F8", "#91BFDB", "#4575B4")))(100),
    cellwidth = 25, cellheight = 3.8,
    treeheight_row = 0,
    treeheight_column=0,
    kmeans_k = NA,
    show_rownames = T, show_colnames = T,
    fontsize=4,
    scale="none",
    legend=FALSE,
    )

    dev.off()

1 个答案:

答案 0 :(得分:1)

cluster_rows=FALSE, cluster_cols=FALSE一般来说,根据@blmoore