# HG changeset patch # User davidvanzessen # Date 1444901657 14400 # Node ID 36e307208f1befde927bc777f061d5ae2922caa9 # Parent 3ed7878f75c333cde045c5e274dff844ac6c468a Uploaded diff -r 3ed7878f75c3 -r 36e307208f1b RScript.r --- a/RScript.r Wed Oct 14 10:21:13 2015 -0400 +++ b/RScript.r Thu Oct 15 05:34:17 2015 -0400 @@ -276,6 +276,21 @@ } } + } else if(nrow(first.rows) > 1) { + if(patient1[1,"Sample"] == first.sample){ + patient1 = patient1[!(patient1$Clone_Sequence %in% first.rows$Clone_Sequence),] + patient1 = rbind(patient1, first.sum) + } else { + patient2 = patient2[!(patient2$Clone_Sequence %in% first.rows$Clone_Sequence),] + patient2 = rbind(patient2, first.sum) + } + + hidden.clone.sequences = c(first.rows[-1,"Clone_Sequence"]) + merge.list[["second"]] = append(merge.list[["second"]], hidden.clone.sequences) + + patient.fuzzy = patient.fuzzy[-first.match.filter,] + + cat(paste("", patient, " row ", 1:nrow(first.rows), "", first.rows$Sample, ":", first.rows$Clone_Sequence, "", first.rows$normalized_read_count, "", sep=""), file="single_matches.html", append=T) } else { patient.fuzzy = patient.fuzzy[-1,] } diff -r 3ed7878f75c3 -r 36e307208f1b wrapper.sh