comparison tmp/igat.r @ 101:3cffb8a38bb1 draft

Uploaded
author davidvanzessen
date Fri, 17 Jun 2016 08:31:20 -0400
parents 86206431cbb0
children e6bc976760d4
comparison
equal deleted inserted replaced
100:ff5be711382b 101:3cffb8a38bb1
3 imgt.dir = args[1] 3 imgt.dir = args[1]
4 merged.file = args[2] 4 merged.file = args[2]
5 gene = args[3] 5 gene = args[3]
6 6
7 merged = read.table(merged.file, header=T, sep="\t", fill=T, stringsAsFactors=F) 7 merged = read.table(merged.file, header=T, sep="\t", fill=T, stringsAsFactors=F)
8
9 print(head(merged$best_match))
8 10
9 if(gene != "-"){ 11 if(gene != "-"){
10 merged = merged[grepl(gene, merged$best_match),] 12 merged = merged[grepl(gene, merged$best_match),]
11 } 13 }
12 14