cole-trapnell-lab / monocle-release

277 stars 116 forks source link

plot_genes_branched_heatmap is stuck #496

Open myshu2017-03-14 opened 1 year ago

myshu2017-03-14 commented 1 year ago

Hi, I recently got stuck (for 3 days) while using plot_genes_branched_heatmap in Monocle2. But other data works fine. I ran the following code:

library(monocle)
cds_data <- readRDS("combined.rds")
load("BEAM_res.Rdata")
BEAM_res2 <- BEAM_res[order(BEAM_res$qval),]

pdf(file="heatmap.pdf")
print(plot_genes_branched_heatmap(cds_data[row.names(head(BEAM_res2, 10)),],branch_point = 1,num_clusters = 4,cores =2,             use_gene_short_name = TRUE,show_rownames = TRUE,return_heatmap = F))
dev.off()

I would like to ask what is the cause of this problem and how can I send the test data to you?