Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

diff RScript.r @ 62:b1ad6f515338 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 03 Nov 2015 10:49:26 -0500
parents 77a090cd0e02
children 39ea37070e90
line wrap: on
line diff
--- 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")