Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

diff RScript.r @ 52:c5c2a790d476 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 09 Oct 2015 10:55:48 -0400
parents 17e677c72e49
children 8ebe57feecd6
line wrap: on
line diff
--- a/RScript.r	Fri Oct 09 06:58:17 2015 -0400
+++ b/RScript.r	Fri Oct 09 10:55:48 2015 -0400
@@ -282,7 +282,8 @@
     cat(paste("<td>", nrow(patient1), " in ", oneSample, " and ", nrow(patient2), " in ", twoSample, ", ", nrow(patientMerge), " in both (finding both took ", (proc.time() - start.time)[[3]], "s)</td></tr>", sep=""), file=logfile, append=T)
   }
 
-  print(names(patient.merge.list))
+  patient1 = patient1[!(patient1$Clone_Sequence %in% patient.merge.list.second[[patient]]),]
+  patient2 = patient2[!(patient2$Clone_Sequence %in% patient.merge.list.second[[patient]]),]
 
   
   patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony])