Add DIY bwCol tSNEs

r.scripts/sc-TissueMapper.R

@@ -867,6 +867,16 @@ sctSNECustCol <- function(sc10x,i,bl,rd){
   dev.off()
 }
 
+sctSNEbwCol <- function(sc10x,i){
+  sc10x <- SetAllIdent(object=sc10x,id=i)
+  postscript(paste0("./analysis/diy/tSNE_",i,".bwCol.eps"))
+  plot <- TSNEPlot(object=sc10x,pt.size=5,do.label=FALSE,label.size=10,do.return=TRUE,vector.friendly=FALSE,colors.use=colorRampPalette(c("white","black"))(length(unique(sc10x@ident))+2)[3:(length(unique(sc10x@ident))+2)])
+  plot <- plot+theme(axis.text.x=element_text(size=20),axis.text.y=element_text(size=20),axis.title.x=element_text(size=20),axis.title.y=element_text(size=20),legend.text=element_text(size=20))
+  plot <- plot+guides(colour=guide_legend(override.aes=list(size=10)))
+  plot(plot)
+  dev.off()
+}
+
 
 sctSNEHighlight <- function(sc10x,i,g){
   sc10x <- SetAllIdent(object=sc10x,id=i)

r.scripts/sc-TissueMapper_RUN.D17.R

@@ -360,6 +360,8 @@ sctSNECustCol(sc10x,i="Merge_Epi.dws_St.go",bl=c("BE","LE","OE1","OE2"),rd=c("Fi
 sctSNECustCol(sc10x.Epi,i="Epi.dws.sub",bl=c("BE","LE","OE1","OE2"),rd="")
 sctSNECustCol(sc10x.St,i="St.go",bl="",rd=c("Fib","SM","Endo","Leu"))
 
+sctSNEbwCol(sc10x,i="res0.2")
+
 for (g in c("Epi","St")){
   sctSNEHighlight(sc10x,i="Lin",g=g)
 }