diff --git a/r.scripts/sc_QuSAGE.Lineage.R b/r.scripts/sc_QuSAGE.Lineage.R index 076e212d0fae5cc1bedc576d6f4d817e3b7da2cc..faeba901ffb275b3f0c0423d7377595547bebdd2 100755 --- a/r.scripts/sc_QuSAGE.Lineage.R +++ b/r.scripts/sc_QuSAGE.Lineage.R @@ -73,11 +73,11 @@ for (i in 1:Number.Clusters){ cC[i] <- paste0("Cluster_",i,"-REST") } rm(labels) -gene.set2 <- read_delim("./genesets/DEG_FMSt_2FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) +gene.set2 <- read_delim("./genesets/DEG_FMSt_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) gene.set2 <- as.list(gene.set2) names(gene.set2) <- "St" gene.set <- c(gene.set2) -gene.set2 <- read_delim("./genesets/DEG_Epi_2FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) +gene.set2 <- read_delim("./genesets/DEG_Epi_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) gene.set2 <- as.list(gene.set2) names(gene.set2) <- "Epi" gene.set <- c(gene.set,gene.set2) diff --git a/r.scripts/sc_QuSAGE_EpiSubClust.R b/r.scripts/sc_QuSAGE_EpiSubClust.R index be4edef44c00f2c9110b656f6486373db3658239..dbfcc066da4462ff7bcc8d4c453249d1fba6c868 100755 --- a/r.scripts/sc_QuSAGE_EpiSubClust.R +++ b/r.scripts/sc_QuSAGE_EpiSubClust.R @@ -90,15 +90,15 @@ if (opt$lgea==TRUE){ sets=c(sets,"lgea") } -gene.set2 <- read_delim("./genesets/DEG_BE_2FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) +gene.set2 <- read_delim("./genesets/DEG_BE_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) gene.set2 <- as.list(gene.set2) names(gene.set2) <- "BE" gene.set <- c(gene.set2) -gene.set2 <- read_delim("./genesets/DEG_LE_2FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) +gene.set2 <- read_delim("./genesets/DEG_LE_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) gene.set2 <- as.list(gene.set2) names(gene.set2) <- "LE" gene.set <- c(gene.set,gene.set2) -gene.set2 <- read_delim("./genesets/DEG_OE_2FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) +gene.set2 <- read_delim("./genesets/DEG_OE_5FC.txt","\t",escape_double=FALSE,trim_ws=TRUE,col_names=FALSE) gene.set2 <- as.list(gene.set2) names(gene.set2) <- "OE" gene.set.dws <- c(gene.set,gene.set2)