# 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])