Develop

r.scripts/sc-TissueMapper_RUN.D17_FACS.R

@@ -167,13 +167,12 @@ if (opt$st==TRUE){
   sc10x <- scCluster(sc10x,pc.use=pc.use.poststress,res.use=0.5,folder="post.stress",red="pca")
 }
 
-gene.set1 <- read_delim("./genesets/genes.deg.Epi.csv",",",escape_double=FALSE,trim_ws=TRUE,col_names=TRUE)
-gene.set1 <- gene.set1[1]
+
+gene.set1 <- read_delim("./genesets/DEG_Epi_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE)
 gene.set1 <- as.list(gene.set1)
 names(gene.set1) <- "Epi"
 gene.set <- c(gene.set1)
-gene.set1 <- read_delim("./genesets/genes.deg.St.csv",",",escape_double=FALSE,trim_ws=TRUE,col_names=TRUE)
-gene.set1 <- gene.set1[1]
+gene.set1 <- read_delim("./genesets/DEG_FMSt_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE)
 gene.set1 <- as.list(gene.set1)
 names(gene.set1) <- "St"
 gene.set <- c(gene.set,gene.set1)