diff --git a/r.scripts/sc-TissueMapper.R b/r.scripts/sc-TissueMapper.R index 533a67e0f120b78aa0c2a5afe13d5efeaf2341da..196bb868929a39fa62e848ef433c1fd63dcb4cd5 100644 --- a/r.scripts/sc-TissueMapper.R +++ b/r.scripts/sc-TissueMapper.R @@ -103,14 +103,17 @@ scCellCycle <- function(sc10x,sub=FALSE){ sc10x <- ScaleData(object=sc10x,display.progress=FALSE,do.par=TRUE,num.cores=45) sc10x <- CellCycleScoring(object=sc10x,s.genes=genes.s,g2m.genes=genes.g2m,set.ident=TRUE) - postscript(paste0(folder,"Ridge_cc.Raw.eps")) - plot <- RidgePlot(object=sc10x,features.plot=c("PCNA","TOP2A","MCM6","MKI67"),y.log=TRUE,nCol=2,do.return=TRUE) - plot(plot) - dev.off() - postscript(paste0(folder,"Violin_cc.Raw.eps")) - plot <- VlnPlot(object=sc10x,features.plot=c("PCNA","TOP2A","MCM6","MKI67"),nCol=2,point.size.use=1,size.title.use=20,x.lab.rot=TRUE) - plot(plot) - dev.off() + tryCatch({ + postscript(paste0(folder,"Ridge_cc.Raw.eps")) + plot <- RidgePlot(object=sc10x,features.plot=c("PCNA","TOP2A","MCM6","MKI67"),y.log=TRUE,nCol=2,do.return=TRUE) + plot(plot) + dev.off() + postscript(paste0(folder,"Violin_cc.Raw.eps")) + plot <- VlnPlot(object=sc10x,features.plot=c("PCNA","TOP2A","MCM6","MKI67"),nCol=2,point.size.use=1,size.title.use=20,x.lab.rot=TRUE) + plot(plot) + dev.off() + },error=function(e){cat("ERROR : ",conditionMessage(e),"/\n")}) + sc10x <- RunPCA(object=sc10x,pc.genes=c(genes.s,genes.g2m),do.print=FALSE,pcs.store=2) postscript(paste0(folder,"PCA_cc.Raw.eps")) plot <- PCAPlot(object=sc10x,do.return=TRUE)