Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison 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
comparison
equal deleted inserted replaced
51:17e677c72e49 52:c5c2a790d476
280 patient.merge.list[[patient]] <<- patientMerge 280 patient.merge.list[[patient]] <<- patientMerge
281 patient.merge.list.second[[patient]] <<- merge.list[["second"]] 281 patient.merge.list.second[[patient]] <<- merge.list[["second"]]
282 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) 282 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)
283 } 283 }
284 284
285 print(names(patient.merge.list)) 285 patient1 = patient1[!(patient1$Clone_Sequence %in% patient.merge.list.second[[patient]]),]
286 patient2 = patient2[!(patient2$Clone_Sequence %in% patient.merge.list.second[[patient]]),]
286 287
287 288
288 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony]) 289 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony])
289 res1 = vector() 290 res1 = vector()
290 res2 = vector() 291 res2 = vector()