diff sequence_overview.r @ 89:480fdd383fdb draft

Uploaded
author davidvanzessen
date Tue, 31 May 2016 08:30:50 -0400
parents d57c624a9aa9
children f0e8dac22c6e
line wrap: on
line diff
--- a/sequence_overview.r	Mon May 30 10:11:20 2016 -0400
+++ b/sequence_overview.r	Tue May 31 08:30:50 2016 -0400
@@ -143,7 +143,7 @@
 NTresult = data.frame(nt=c("A", "C", "T", "G"))
 
 for(clazz in gene.classes){
-	NToverview.sub = NToverview[grepl(clazz, NToverview$best_match),]
+	NToverview.sub = NToverview[grepl(clazz, paste("^", NToverview$best_match, sep="")),]
 	new.col.x = c(sum(NToverview.sub$A), sum(NToverview.sub$C), sum(NToverview.sub$T), sum(NToverview.sub$G))
 	new.col.y = sum(new.col.x)
 	new.col.z = round(new.col.x / new.col.y * 100, 2)