Mercurial > repos > davidvanzessen > report_clonality_igg

comparison RScript.r @ 26:01f05258f672 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 09 Feb 2015 08:53:30 -0500
parents ea5c2a2cc1f3
children f919965e8816
comparison
equal deleted inserted replaced
25:ea5c2a2cc1f3 26:01f05258f672
69 #remove duplicates based on the clonaltype 69 #remove duplicates based on the clonaltype
70 if(clonaltype != "none"){ 70 if(clonaltype != "none"){
71 clonaltype = paste(clonaltype, ",Sample", sep="") #add sample column to clonaltype, unique within samples 71 clonaltype = paste(clonaltype, ",Sample", sep="") #add sample column to clonaltype, unique within samples
72 PRODF$clonaltype = do.call(paste, c(PRODF[unlist(strsplit(clonaltype, ","))], sep = ":")) 72 PRODF$clonaltype = do.call(paste, c(PRODF[unlist(strsplit(clonaltype, ","))], sep = ":"))
73 PRODF = PRODF[!duplicated(PRODF$clonaltype), ] 73 PRODF = PRODF[!duplicated(PRODF$clonaltype), ]
74
75 UNPROD$clonaltype = do.call(paste, c(UNPROD[unlist(strsplit(clonaltype, ","))], sep = ":"))
76 UNPROD = UNPROD[!duplicated(UNPROD$clonaltype), ]
74 77
75 #again for clonalityFrame but with sample+replicate 78 #again for clonalityFrame but with sample+replicate
76 clonalityFrame$clonaltype = do.call(paste, c(clonalityFrame[unlist(strsplit(clonaltype, ","))], sep = ":")) 79 clonalityFrame$clonaltype = do.call(paste, c(clonalityFrame[unlist(strsplit(clonaltype, ","))], sep = ":"))
77 clonalityFrame$clonality_clonaltype = do.call(paste, c(clonalityFrame[unlist(strsplit(paste(clonaltype, ",Replicate", sep=""), ","))], sep = ":")) 80 clonalityFrame$clonality_clonaltype = do.call(paste, c(clonalityFrame[unlist(strsplit(paste(clonaltype, ",Replicate", sep=""), ","))], sep = ":"))
78 clonalityFrame = clonalityFrame[!duplicated(clonalityFrame$clonality_clonaltype), ] 81 clonalityFrame = clonalityFrame[!duplicated(clonalityFrame$clonality_clonaltype), ]