Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

diff RScript.r @ 43:2eef2c438500 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 14 Sep 2015 11:01:26 -0400
parents 7d4eda663c29
children 07278582b735
line wrap: on
line diff
--- a/RScript.r	Mon Sep 14 10:56:05 2015 -0400
+++ b/RScript.r	Mon Sep 14 11:01:26 2015 -0400
@@ -179,7 +179,9 @@
       
       length.filter = nchar(patient.fuzzy$Clone_Sequence) - nchar(first.clone.sequence) <= 9
       
-      match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & length.filter
+      sample.filter = patient.fuzzy[1,"Sample"] != patient.fuzzy$Sample
+      
+      match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & length.filter & sample.filter
       
       if(sum(match.filter) == 2){
         second.match = which(match.filter)[2]