Mercurial > repos > davidvanzessen > test_plotting_merged

comparison RScript.r @ 7:ba6c6725b12f draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 28 Oct 2013 11:11:46 -0400
parents f9a657db7af5
children
comparison
equal deleted inserted replaced
6:f9a657db7af5 7:ba6c6725b12f
31 31
32 NONPROD = test[test$VDJ.Frame == "In-frame with stop codon" | test$VDJ.Frame == "Out-of-frame" | test$CDR3.Found.How == "NOT_FOUND" , ] 32 NONPROD = test[test$VDJ.Frame == "In-frame with stop codon" | test$VDJ.Frame == "Out-of-frame" | test$CDR3.Found.How == "NOT_FOUND" , ]
33 33
34 PRODF = PROD[ -1] 34 PRODF = PROD[ -1]
35 35
36 PRODF = unique(PRODF) 36 #PRODF = unique(PRODF)
37 PRODF = PRODF[!duplicated(PRODF$VDJCDR3), ]
38
37 39
38 uniqueCount = split(PRODF, f=PRODF[,"Sample"]) 40 uniqueCount = split(PRODF, f=PRODF[,"Sample"])
39 41
40 for(i in 1:length(uniqueCount)) { 42 for(i in 1:length(uniqueCount)) {
41 dat = data.frame(uniqueCount[i]) 43 dat = data.frame(uniqueCount[i])