When use DWS NE geneset, use 100 bins (instead of 50) for NE-Score (PC1) turning point choice

r.scripts/sc_PC.Score.NE.R

@@ -89,7 +89,11 @@ dev.off()
 
 #find first turning point on a KDE of NE1
 #hist(GetCellEmbeddings(object=sc10x.Group.Epi,reduction.type="NE",dims.use=1),prob=TRUE)
-d1 <- density(GetCellEmbeddings(object=sc10x.Group.Epi,reduction.type="NE",dims.use=1),n=50)
+if (opt$u==FALSE){
+  d1 <- density(GetCellEmbeddings(object=sc10x.Group.Epi,reduction.type="NE",dims.use=1),n=50)
+} else {
+  d1 <- density(GetCellEmbeddings(object=sc10x.Group.Epi,reduction.type="NE",dims.use=1),n=100)
+}
 #d2 <- density(GetCellEmbeddings(object=sc10x.Group.Epi,reduction.type="NE",dims.use=2))
 #lines(d1,col="blue")
 ts_1 <- ts(d1$y)