Change stress score binning to 100

r.scripts/sc-TissueMapper.R

@@ -350,8 +350,8 @@ scStress <- function(sc10x,stg="go",res.use=1,pc.use=10){
   
   #Find first turning point on a KDE of Stress1
   postscript("./analysis/pca/stress/Histo_Stress.eps")
-  histo <- hist(GetCellEmbeddings(object=sc10x,reduction.type="Stress",dims.use=1),breaks=50,prob=TRUE,plot=TRUE)
-  d1 <- density(GetCellEmbeddings(object=sc10x,reduction.type="Stress",dims.use=1),n=50)
+  histo <- hist(GetCellEmbeddings(object=sc10x,reduction.type="Stress",dims.use=1),breaks=100,prob=TRUE,plot=TRUE)
+  d1 <- density(GetCellEmbeddings(object=sc10x,reduction.type="Stress",dims.use=1),n=100)
   lines(d1,col="blue")
   ts_1 <- ts(d1$y)
   tp1 <- turnpoints(ts_1)