diff --git a/bash.scripts/sc_TissueMapper-DPrF.sh b/bash.scripts/sc_TissueMapper-DPrF.sh
index c38f8a99a337eaec700828c452fcf444ca0e0ea0..ae34cb34d687d4355a68f014a245ae74757c1ba9 100644
--- a/bash.scripts/sc_TissueMapper-DPrF.sh
+++ b/bash.scripts/sc_TissueMapper-DPrF.sh
@@ -10,4 +10,6 @@
 
 module load R/3.4.1-gccmkl
 
-Rscript ../r.scripts/sc-TissueMapper_RUN.D.R
\ No newline at end of file
+Rscript ../r.scripts/sc-TissueMapper_RUN.D.R
+Rscript ../r.scripts/sc-TissueMapper_RUN.D.diy.R
+Rscript ../r.scripts/sc-TissueMapper_RUN.D.pseudotime.R
\ No newline at end of file
diff --git a/r.scripts/sc-TissueMapper.R b/r.scripts/sc-TissueMapper.R
index 1ee7a08ac449bc552afb8dff74843d6d2345f7a8..cf280c0b663c50b1fbbed2c67f86f4480b7fa95a 100644
--- a/r.scripts/sc-TissueMapper.R
+++ b/r.scripts/sc-TissueMapper.R
@@ -1133,6 +1133,5 @@ scPseudotime <- function(sc10x,i,ds){
   BEAM_res <- BEAM(scMonocle,branch_point=1,cores=50)
   BEAM_res <- BEAM_res[order(BEAM_res$qval),]
   BEAM_res <- BEAM_res[,c("gene_short_name","pval","qval")]
-  plot_genes_branched_heatmap(scMonocle[rownames(subset(BEAM_res,qval<1e-4)),],branch_point=1,num_clusters=5,cores=50,use_gene_short_name=TRUE,show_rownames=TRUE)
-  
-}
\ No newline at end of file
+  #plot_genes_branched_heatmap(scMonocle[rownames(subset(BEAM_res,qval<1e-4)),],branch_point=1,num_clusters=5,cores=50,use_gene_short_name=TRUE,show_rownames=TRUE)
+}