diff --git a/r.scripts/Pd_convert_to_UMAP.R b/r.scripts/Pd_convert_to_UMAP.R
new file mode 100644
index 0000000000000000000000000000000000000000..a5f96189072b576cc3e605cf1b2c5e1b85fc9d8d
--- /dev/null
+++ b/r.scripts/Pd_convert_to_UMAP.R
@@ -0,0 +1,32 @@
+###RUN TO START
+gc()
+library(methods)
+library(optparse)
+library(Seurat)
+library(RColorBrewer)
+library(viridis)
+options(bitmapType="cairo")
+
+setwd("../")
+
+source("./r.scripts/sc-TissueMapper_functions.R")
+source("./r.scripts/sc-TissueMapper_process.R")
+
+load("./analysis/sc10x.Rda")
+try(
+  if (as.numeric(substring(sc10x@version,1,1))<3){
+    sc10x <- UpdateSeuratObject(sc10x)
+  }
+)
+
+sc10x <- RunUMAP(sc10x,dims=1:30,reduction="cca.aligned")
+
+Idents(object=sc10x) <- "Merge_Epi.dws_St.go_NE"
+sc10x <- RenameIdents(object=sc10x, "OE1" = "Club")
+sc10x <- RenameIdents(object=sc10x, "OE2" = "Hillock")
+
+postscript(paste0("./analysis/UMAP.aggr.eps"))
+DimPlot(sc10x,reduction="umap",label=FALSE,order=rev(c("BE","LE","Hillock","Club","NE","Fib","SM","Endo","Leu")),cols=rev(c(brewer.pal(5,"Reds")[2:5],"grey50",brewer.pal(5,"Blues")[2:5])))
+dev.off()
+
+write.table(Embeddings(object = sc10x, reduction = "umap"),file=paste0("./analysis/UMAP.aggr.csv"),row.names=TRUE,col.names=NA,append=FALSE,sep=",")