Mercurial > repos > davidvanzessen > mutation_analysis

comparison mutation_analysis.r @ 33:ac9a4307861a draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Thu, 16 Apr 2015 09:36:19 -0400
parents 2433a1e110e1
children e022c21f8c47
comparison
equal deleted inserted replaced
32:2a7343e4be5a 33:ac9a4307861a
282 cat(length(tmp$Sequence.ID), file="total_n.txt") 282 cat(length(tmp$Sequence.ID), file="total_n.txt")
283 283
284 284
285 285
286 result = data.frame(matrx) 286 result = data.frame(matrx)
287 row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR S/R (ratio)", "CDR S/R (ratio)") 287 row.names(result) = c("Number of Mutations (%)", "Transition (%)", "Transversions (%)", "Transitions at G C (%)", "Targeting of C.G (%)", "FR R/S (ratio)", "CDR R/S (ratio)")
288 288
289 write.table(x=result, file="mutations.txt", sep=",",quote=F,row.names=T,col.names=F) 289 write.table(x=result, file="mutations.txt", sep=",",quote=F,row.names=T,col.names=F)
290 290
291 291
292 if (!("ggplot2" %in% rownames(installed.packages()))) { 292 if (!("ggplot2" %in% rownames(installed.packages()))) {