Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison RScript.r @ 45:570fef3d876f draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 16 Sep 2015 05:31:09 -0400
parents 07278582b735
children bec38379cf69
comparison
equal deleted inserted replaced
44:07278582b735 45:570fef3d876f
182 sample.filter = patient.fuzzy[1,"Sample"] != patient.fuzzy$Sample 182 sample.filter = patient.fuzzy[1,"Sample"] != patient.fuzzy$Sample
183 183
184 #match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & length.filter & sample.filter 184 #match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & length.filter & sample.filter
185 match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & sample.filter 185 match.filter = merge.filter & grepl(first.clone.sequence, patient.fuzzy$Clone_Sequence) & sample.filter
186 186
187 if(sum(match.filter) == 2){ 187 if(sum(match.filter) == 1){
188 second.match = which(match.filter)[2] 188 second.match = which(match.filter)[1]
189 second.clone.sequence = patient.fuzzy[second.match,"Clone_Sequence"] 189 second.clone.sequence = patient.fuzzy[second.match,"Clone_Sequence"]
190 first.sample = patient.fuzzy[1,"Sample"] 190 first.sample = patient.fuzzy[1,"Sample"]
191 second.sample = patient.fuzzy[second.match,"Sample"] 191 second.sample = patient.fuzzy[second.match,"Sample"]
192 192
193 if(((nchar(second.clone.sequence) - nchar(first.clone.sequence)) <= 9) & (first.sample != second.sample)){ 193 if(((nchar(second.clone.sequence) - nchar(first.clone.sequence)) <= 9) & (first.sample != second.sample)){