# HG changeset patch
# User davidvanzessen
# Date 1444402548 14400
# Node ID c5c2a790d4761a4384455c04e23e676b5d77e50e
# Parent 17e677c72e49cbc60dbdaa89d1f44d46cc5a4607
Uploaded
diff -r 17e677c72e49 -r c5c2a790d476 RScript.r
--- 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("
", nrow(patient1), " in ", oneSample, " and ", nrow(patient2), " in ", twoSample, ", ", nrow(patientMerge), " in both (finding both took ", (proc.time() - start.time)[[3]], "s) | ", 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])