Mercurial > repos > davidvanzessen > mutation_analysis

comparison mutation_analysis.r @ 94:e39176ccddc8 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 07 Jun 2016 08:39:08 -0400
parents b523ce95d857
children 5ffbf40cdd4b
comparison
equal deleted inserted replaced
93:53fb2948726e 94:e39176ccddc8
120 120
121 VRegionMutations_columns = paste(regions, ".IMGT.Nb.of.mutations", sep="") 121 VRegionMutations_columns = paste(regions, ".IMGT.Nb.of.mutations", sep="")
122 dat$VRegionMutations = apply(dat, FUN=sum_by_row, 1, columns=VRegionMutations_columns) 122 dat$VRegionMutations = apply(dat, FUN=sum_by_row, 1, columns=VRegionMutations_columns)
123 123
124 VRegionNucleotides_columns = paste(regions, ".IMGT.Nb.of.nucleotides", sep="") 124 VRegionNucleotides_columns = paste(regions, ".IMGT.Nb.of.nucleotides", sep="")
125 dat$FR3.IMGT.Nb.of.nucleotides = nchar(dat$FR3.IMGT.seq)
125 dat$VRegionNucleotides = apply(dat, FUN=sum_by_row, 1, columns=VRegionNucleotides_columns) 126 dat$VRegionNucleotides = apply(dat, FUN=sum_by_row, 1, columns=VRegionNucleotides_columns)
126 127
127 transitionMutations_columns = paste(rep(regions, each=4), c(".IMGT.a.g", ".IMGT.g.a", ".IMGT.c.t", ".IMGT.t.c"), sep="") 128 transitionMutations_columns = paste(rep(regions, each=4), c(".IMGT.a.g", ".IMGT.g.a", ".IMGT.c.t", ".IMGT.t.c"), sep="")
128 dat$transitionMutations = apply(dat, FUN=sum_by_row, 1, columns=transitionMutations_columns) 129 dat$transitionMutations = apply(dat, FUN=sum_by_row, 1, columns=transitionMutations_columns)
129 130