From 52de10530cfd98c52f52fed69b65639431a650d0 Mon Sep 17 00:00:00 2001 From: Efrain Gonzalez Date: Thu, 2 Nov 2017 10:32:00 -0400 Subject: [PATCH] Updated: added |^GENE_SYMBOL$ to glossary --- RCleanDscret.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/RCleanDscret.R b/RCleanDscret.R index 17bd6fb..f05f745 100644 --- a/RCleanDscret.R +++ b/RCleanDscret.R @@ -258,7 +258,7 @@ geneIDNam <- paste0("Clean_GPL",gplnum,".txt") %>% .[IDLOCAL] geneIDNam <- genena %>% read_delim(delim="\t",col_names = TRUE, comment = "!", skip = idlocgpl) %>% - dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$",colnames(.))) + dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$|^GENE_SYMBOL$",colnames(.))) } else if(IDF == 0){ #No information on this particular GPL file idLOCGPL <- genena %>% @@ -270,7 +270,7 @@ geneIDNam <- paste0("Clean_GPL",gplnum,".txt") %>% cat(file="GPL_ID_LOC.txt",sep = "\t", fill = TRUE, append = TRUE) geneIDNam <- genena %>% read_delim(delim="\t",col_names = TRUE, comment = "!", skip = idLOCGPL) %>% - dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$",colnames(.))) + dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$|^GENE_SYMBOL$",colnames(.))) } } else if(fileex == 0){ #We must create a file that we can access for later use @@ -284,12 +284,12 @@ geneIDNam <- paste0("Clean_GPL",gplnum,".txt") %>% write.table(Firstval,file = "GPL_ID_LOC.txt", sep = "\t",row.names = FALSE, col.names = TRUE) geneIDNam <- genena %>% read_delim(delim="\t",col_names = TRUE, comment = "!", skip = idLOCGPL) %>% - dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$",colnames(.))) + dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$|^GENE_SYMBOL$",colnames(.))) } } else if(soft == FALSE){ geneIDNam <- genena %>% read_delim(delim="\t",comment = "#")%>% - dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$",colnames(.))) + dplyr::select(.,ID,grep("Symbol|^ORF\\s*$|^gene_assignment\\s*$|^Gene symbol$|^GENE_SYMBOL$",colnames(.))) } ##Labeling the gene IDs without names -- 2.29.0