comparison tmp/igat.r @ 102:e6bc976760d4 draft

Uploaded
author davidvanzessen
date Tue, 21 Jun 2016 03:32:50 -0400
parents 3cffb8a38bb1
children e4957ad476a2
comparison
equal deleted inserted replaced
101:3cffb8a38bb1 102:e6bc976760d4
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))
10 8
11 if(gene != "-"){ 9 if(gene != "-"){
12 merged = merged[grepl(gene, merged$best_match),] 10 merged = merged[grepl(gene, merged$best_match),]
13 } 11 }
14 12