diff --git a/r.scripts/sc_PC.Score.NE.R b/r.scripts/sc_PC.Score.NE.R
index bf657dcd3f1fd7534f2892ca4104f1b7adbd61ae..a6c5e0012adfe4790a78a0ce4ce041ad9703cb44 100755
--- a/r.scripts/sc_PC.Score.NE.R
+++ b/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)