remove genes which are Unknown or "Not Imprinted"

[imprinted_genes.git] / combine_imprinted_genes.R

diff --git a/combine_imprinted_genes.R b/combine_imprinted_genes.R
index 7b75d523fd5cd0bcc440435b5b1ab73491d13e1b..c1e9271716df4b266abf62b69e182d11db04c4a8 100644 (file)
--- a/combine_imprinted_genes.R
+++ b/combine_imprinted_genes.R
@@ -4,6 +4,10 @@ args <- commandArgs(trailingOnly=TRUE)
 
 geneimprint <- fread(args[1])
 geneimprint <- geneimprint[!grepl(" ",Gene),]
+
+geneimprint <- geneimprint[Status!="Not Imprinted",]
+geneimprint <- geneimprint[Status!="Unknown",]
+
 parent <- fread(args[2])
 ### fix up the 0 prefixed chromosomes
 parent[,chr:=gsub("^0","",chromosome)]