Skip to content
Snippets Groups Projects
Commit 1debf570 authored by Gervaise Henry's avatar Gervaise Henry :cowboy:
Browse files

Fix run_Pd for Hillock geneset

parent 97e6ecfc
Branches
Tags
1 merge request!6Develop
......@@ -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])
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment