Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

diff RScript.r @ 63:39ea37070e90 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 04 Nov 2015 11:02:01 -0500
parents b1ad6f515338
children 39fff0180d41
line wrap: on
line diff
--- a/RScript.r	Tue Nov 03 10:49:26 2015 -0500
+++ b/RScript.r	Wed Nov 04 11:02:01 2015 -0500
@@ -714,6 +714,10 @@
         hidden.clone.sequences = c(rows.1[-1,"Clone_Sequence"], rows.3[rows.3$Clone_Sequence != first.clone.sequence,"Clone_Sequence"])
         merge.list[["second"]] = append(merge.list[["second"]], hidden.clone.sequences)
 
+      } else if(nrow(rows.1) > 1){
+        patient1 = patient1[!(patient1$Clone_Sequence %in% rows.1$Clone_Sequence),]
+        patient1 = rbind(patient1, sum.1)
+        patient.fuzzy = patient.fuzzy[-match.filter.1,]
       } else {
         patient.fuzzy = patient.fuzzy[-1,]
       }
@@ -745,9 +749,9 @@
   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]]),]
+  patient1 = patient1[!(patient1$Clone_Sequence %in% merge.list[["second"]]),]
+  patient2 = patient2[!(patient2$Clone_Sequence %in% merge.list[["second"]]),]
+  patient3 = patient3[!(patient3$Clone_Sequence %in% merge.list[["second"]]),]
 
   if(F){
     patientMerge = merge(patient1, patient2, by="merge")