Mercurial > repos > davidvanzessen > mutation_analysis

comparison sequence_overview.r @ 78:b523ce95d857 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Wed, 11 May 2016 10:29:33 -0400
parents c5c86d15cb94
children 0513b46178c4
comparison
equal deleted inserted replaced
77:c5c86d15cb94 78:b523ce95d857
14 14
15 dat = merge(sequences, genes, by="Sequence.ID") 15 dat = merge(sequences, genes, by="Sequence.ID")
16 16
17 dat$seq_conc = paste(dat$CDR1.IMGT, dat$CDR2.IMGT, dat$CDR3.IMGT, dat$FR2.IMGT, dat$FR3.IMGT) 17 dat$seq_conc = paste(dat$CDR1.IMGT, dat$CDR2.IMGT, dat$CDR3.IMGT, dat$FR2.IMGT, dat$FR3.IMGT)
18 18
19 IDs = dat[,c("Sequence.ID", "seq_conc", "best_match")] 19 IDs = dat[,c("Sequence.ID", "seq_conc", "best_match", "Functionality")]
20 IDs$best_match = as.character(IDs$best_match) 20 IDs$best_match = as.character(IDs$best_match)
21 21
22 #dat = data.frame(data.table(dat)[, list(freq=.N), by=c("best_match", "seq_conc")]) 22 #dat = data.frame(data.table(dat)[, list(freq=.N), by=c("best_match", "seq_conc")])
23 23
24 dat = data.frame(table(dat$best_match, dat$seq_conc)) 24 dat = data.frame(table(dat$best_match, dat$seq_conc))