Mercurial > repos > davidvanzessen > argalaxy_tools

comparison imgt_loader/imgt_loader.r @ 31:3a76faa53c59 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 14 Mar 2017 06:51:38 -0400
parents 798b62942b4b
children bad6a9a53ae7
comparison
equal deleted inserted replaced
30:b50965edac24 31:3a76faa53c59
29 out[,"CDR1.Length"] = summ[,"CDR1.IMGT.length"] 29 out[,"CDR1.Length"] = summ[,"CDR1.IMGT.length"]
30 30
31 out[,"CDR2.Seq"] = aa[,"CDR2.IMGT"] 31 out[,"CDR2.Seq"] = aa[,"CDR2.IMGT"]
32 out[,"CDR2.Length"] = summ[,"CDR2.IMGT.length"] 32 out[,"CDR2.Length"] = summ[,"CDR2.IMGT.length"]
33 33
34 out[,"CDR3.Seq"] = aa[,"CDR3.IMGT"] 34 out[,"CDR3.Seq"] = gapped.aa[,"CDR3.IMGT"]
35 out[,"CDR3.Length"] = summ[,"CDR3.IMGT.length"] 35 out[,"CDR3.Length"] = summ[,"CDR3.IMGT.length"]
36 36
37 out[,"CDR3.Seq.DNA"] = gapped.aa[,"CDR3.IMGT"] 37 out[,"CDR3.IMGT"] = out[,"CDR3.Seq"]
38
39 out[,"CDR3.Seq.DNA"] = sequences[,"CDR3.IMGT"]
38 out[,"CDR3.Length.DNA"] = nchar(as.character(out[,"CDR3.Seq.DNA"])) 40 out[,"CDR3.Length.DNA"] = nchar(as.character(out[,"CDR3.Seq.DNA"]))
39 out[,"Strand"] = summ[,"Orientation"] 41 out[,"Strand"] = summ[,"Orientation"]
40 out[,"CDR3.Found.How"] = "a" 42 out[,"CDR3.Found.How"] = "a"
41 43
42 out[,added_summary_columns] = summ[,added_summary_columns] 44 out[,added_summary_columns] = summ[,added_summary_columns]