diff --git a/r.scripts/sc-TissueMapper_RUN.Pd.R b/r.scripts/sc-TissueMapper_RUN.Pd.R index 81cbf9a7a8ee33bf0dba5a209fdd5a4eefc06968..396cf44665e71c349c817667729f5b8b72afacea 100644 --- a/r.scripts/sc-TissueMapper_RUN.Pd.R +++ b/r.scripts/sc-TissueMapper_RUN.Pd.R @@ -377,7 +377,7 @@ gene.set2 <- as.list(gene.set2) names(gene.set2) <- "Ionocyte" gene.set <- c(gene.set,gene.set2) rm(gene.set2) -gene.set1 <- read_csv("./genesets/SupTab6_Krt13_Hillock.csv",skip=5) +gene.set1 <- read_csv("./genesets/SupTab6_Krt13_Hillock.csv") gene.set1 <- gene.set1[gene.set1$FDR<=0.05 & gene.set1$'log2 fold-change (MAST)'>=1.5,1] colnames(gene.set1) <- "genes" gene.set1 <- as.data.frame(merge(gene.set1,gene.orthog,by.x="genes",by.y="Gene.name")[,4])