Skip to content
Snippets Groups Projects
Commit fc999941 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Add secondary analysis to DPr bash script

parent fad1aca5
Branches
Tags
No related merge requests found
Hide whitespace changes
Inline Side-by-side
bash.scripts/sc_TissueMapper-DPrF.sh
+ 3
− 1
View file @ fc999941
...@@ -10,4 +10,6 @@ ...@@ -10,4 +10,6 @@
module load R/3.4.1-gccmkl module load R/3.4.1-gccmkl
Rscript ../r.scripts/sc-TissueMapper_RUN.D.R Rscript ../r.scripts/sc-TissueMapper_RUN.D.R
\ No newline at end of file Rscript ../r.scripts/sc-TissueMapper_RUN.D.diy.R
Rscript ../r.scripts/sc-TissueMapper_RUN.D.pseudotime.R
\ No newline at end of file
r.scripts/sc-TissueMapper.R
+ 2
− 3
View file @ fc999941
...@@ -1133,6 +1133,5 @@ scPseudotime <- function(sc10x,i,ds){ ...@@ -1133,6 +1133,5 @@ scPseudotime <- function(sc10x,i,ds){
BEAM_res <- BEAM(scMonocle,branch_point=1,cores=50) BEAM_res <- BEAM(scMonocle,branch_point=1,cores=50)
BEAM_res <- BEAM_res[order(BEAM_res$qval),] BEAM_res <- BEAM_res[order(BEAM_res$qval),]
BEAM_res <- BEAM_res[,c("gene_short_name","pval","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) #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
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment