# HG changeset patch # User davidvanzessen # Date 1464076611 14400 # Node ID e1114c791eb7a90bd85502e9a01d6197c40d2cca # Parent 564c4f6da20301b40c3e8a3e963959577b821e5a Uploaded diff -r 564c4f6da203 -r e1114c791eb7 merge_and_filter.r --- a/merge_and_filter.r Tue May 24 03:52:14 2016 -0400 +++ b/merge_and_filter.r Tue May 24 03:56:51 2016 -0400 @@ -66,7 +66,7 @@ summ[!higher_than,"best_match"] = paste("unmatched,", summ[!higher_than,"best_match"]) } -if(any(higher_than)){ +if(any(higher_than, na.rm=T)){ #summ = summ[higher_than,] } print(paste("Number of matched sequences:", sum(!grepl("^unmatched", summ$best_match))))