# HG changeset patch
# User davidvanzessen
# Date 1446565766 18000
# Node ID b1ad6f515338086818bad078859ee703d1abb303
# Parent  77a090cd0e028dce8b4db0b5c43b63139d3c12bb
Uploaded

diff -r 77a090cd0e02 -r b1ad6f515338 RScript.r
--- a/RScript.r	Fri Oct 30 04:39:33 2015 -0400
+++ b/RScript.r	Tue Nov 03 10:49:26 2015 -0500
@@ -745,6 +745,10 @@
   patientMerge13$thresholdValue = pmax(patientMerge13[,onx], patientMerge13[,ony])
   patientMerge23$thresholdValue = pmax(patientMerge23[,onx], patientMerge23[,ony])
 
+  patient1 = patient1[!(patient1$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
+  patient2 = patient2[!(patient2$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
+  patient2 = patient3[!(patient3$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
+
   if(F){
     patientMerge = merge(patient1, patient2, by="merge")
     patientMerge = merge(patientMerge, patient3, by="merge")