Mercurial > repos > davidvanzessen > mutation_analysis

diff mutation_analysis.r @ 40:e022c21f8c47 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 02 Nov 2015 07:10:36 -0500
parents ac9a4307861a
children 9afd8430de2c
line wrap: on
line diff
--- a/mutation_analysis.r	Mon Nov 02 04:52:26 2015 -0500
+++ b/mutation_analysis.r	Mon Nov 02 07:10:36 2015 -0500
@@ -145,6 +145,9 @@
 CDR_nonSilentMutations_columns = paste(CDRRegions, ".IMGT.Nb.of.nonsilent.mutations", sep="")
 dat$nonSilentMutationsCDR = apply(dat, FUN=sum_by_row, 1, columns=CDR_nonSilentMutations_columns)
 
+mutation.sum.columns = c("Sequence.ID", "VRegionMutations", "VRegionNucleotides", "transitionMutations", "transversionMutations", "transitionMutationsAtGC", "silentMutationsFR", "nonSilentMutationsFR", "silentMutationsCDR", "nonSilentMutationsCDR")
+
+write.table(dat[,mutation.sum.columns], "mutation_by_id.txt", sep="\t",quote=F,row.names=F,col.names=T)
 
 setwd(outputdir)