# HG changeset patch # User davidvanzessen # Date 1462350482 14400 # Node ID 51d92233fb5d7b4ea6b376422c237823635471ba # Parent 9165bec41c028fdc11346278a17c2333c5e27acd Uploaded diff -r 9165bec41c02 -r 51d92233fb5d mutation_analysis.r --- a/mutation_analysis.r Wed May 04 03:48:57 2016 -0400 +++ b/mutation_analysis.r Wed May 04 04:28:02 2016 -0400 @@ -229,7 +229,7 @@ matrx[10,y] = round(f(tmp$VRegionNucleotides, na.rm=T), digits=1) matrx[10,z] = round(matrx[10,x] / matrx[10,y], digits=1) - matrx[11,x] = round(f(rowSums(tmp[,c("CDR1.IMGT.Nb.of.nucleotides", "CDR3.IMGT.Nb.of.nucleotides")], na.rm=T)), digits=1) + matrx[11,x] = round(f(rowSums(tmp[,c("CDR1.IMGT.Nb.of.nucleotides", "CDR2.IMGT.Nb.of.nucleotides")], na.rm=T)), digits=1) matrx[11,y] = round(f(tmp$VRegionNucleotides, na.rm=T), digits=1) matrx[11,z] = round(matrx[11,x] / matrx[11,y], digits=1) }