Mercurial > repos > davidvanzessen > clonal_sequences_in_paired_samples

comparison RScript.r @ 62:b1ad6f515338 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 03 Nov 2015 10:49:26 -0500
parents 77a090cd0e02
children 39ea37070e90
comparison
equal deleted inserted replaced
61:77a090cd0e02 62:b1ad6f515338
743 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony], patientMerge[,onz]) 743 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony], patientMerge[,onz])
744 patientMerge12$thresholdValue = pmax(patientMerge12[,onx], patientMerge12[,ony]) 744 patientMerge12$thresholdValue = pmax(patientMerge12[,onx], patientMerge12[,ony])
745 patientMerge13$thresholdValue = pmax(patientMerge13[,onx], patientMerge13[,ony]) 745 patientMerge13$thresholdValue = pmax(patientMerge13[,onx], patientMerge13[,ony])
746 patientMerge23$thresholdValue = pmax(patientMerge23[,onx], patientMerge23[,ony]) 746 patientMerge23$thresholdValue = pmax(patientMerge23[,onx], patientMerge23[,ony])
747 747
748 patient1 = patient1[!(patient1$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
749 patient2 = patient2[!(patient2$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
750 patient2 = patient3[!(patient3$Clone_Sequence %in% patient.merge.list.second[[label1]]),]
751
748 if(F){ 752 if(F){
749 patientMerge = merge(patient1, patient2, by="merge") 753 patientMerge = merge(patient1, patient2, by="merge")
750 patientMerge = merge(patientMerge, patient3, by="merge") 754 patientMerge = merge(patientMerge, patient3, by="merge")
751 colnames(patientMerge)[which(!grepl("(\\.x$)|(\\.y$)|(merge)", names(patientMerge)))] = paste(colnames(patientMerge)[which(!grepl("(\\.x$)|(\\.y$)|(merge)", names(patientMerge), perl=T))], ".z", sep="") 755 colnames(patientMerge)[which(!grepl("(\\.x$)|(\\.y$)|(merge)", names(patientMerge)))] = paste(colnames(patientMerge)[which(!grepl("(\\.x$)|(\\.y$)|(merge)", names(patientMerge), perl=T))], ".z", sep="")
752 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony], patientMerge[,onz]) 756 patientMerge$thresholdValue = pmax(patientMerge[,onx], patientMerge[,ony], patientMerge[,onz])