Mercurial > repos > davidvanzessen > mutation_analysis

comparison sequence_overview.r @ 89:480fdd383fdb draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 31 May 2016 08:30:50 -0400
parents d57c624a9aa9
children f0e8dac22c6e
comparison
equal deleted inserted replaced
88:d57c624a9aa9 89:480fdd383fdb
141 #NToverview = rbind(NToverview, NTsum) 141 #NToverview = rbind(NToverview, NTsum)
142 142
143 NTresult = data.frame(nt=c("A", "C", "T", "G")) 143 NTresult = data.frame(nt=c("A", "C", "T", "G"))
144 144
145 for(clazz in gene.classes){ 145 for(clazz in gene.classes){
146 NToverview.sub = NToverview[grepl(clazz, NToverview$best_match),] 146 NToverview.sub = NToverview[grepl(clazz, paste("^", NToverview$best_match, sep="")),]
147 new.col.x = c(sum(NToverview.sub$A), sum(NToverview.sub$C), sum(NToverview.sub$T), sum(NToverview.sub$G)) 147 new.col.x = c(sum(NToverview.sub$A), sum(NToverview.sub$C), sum(NToverview.sub$T), sum(NToverview.sub$G))
148 new.col.y = sum(new.col.x) 148 new.col.y = sum(new.col.x)
149 new.col.z = round(new.col.x / new.col.y * 100, 2) 149 new.col.z = round(new.col.x / new.col.y * 100, 2)
150 150
151 tmp = names(NTresult) 151 tmp = names(NTresult)