Mercurial > repos > davidvanzessen > mutation_analysis

changeset 87:fd63f5307962 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 30 May 2016 06:49:13 -0400
parents 0011a0597685
children d57c624a9aa9
files sequence_overview.r
diffstat 1 files changed, 4 insertions(+), 3 deletions(-) [+]
line wrap: on
line diff
--- a/sequence_overview.r	Mon May 30 05:36:53 2016 -0400
+++ b/sequence_overview.r	Mon May 30 06:49:13 2016 -0400
@@ -23,7 +23,8 @@
 
 dat = dat[dat$Sequence.ID %in% merged$Sequence.ID,]
 
-dat$seq_conc = paste(dat$CDR1.IMGT, dat$CDR2.IMGT, dat$CDR3.IMGT, dat$FR2.IMGT, dat$FR3.IMGT)
+#dat$seq_conc = paste(dat$CDR1.IMGT, dat$CDR2.IMGT, dat$CDR3.IMGT, dat$FR2.IMGT, dat$FR3.IMGT)
+dat$seq_conc = paste(dat$CDR1.IMGT, dat$CDR2.IMGT, dat$CDR3.IMGT)
 
 IDs = dat[,c("Sequence.ID", "seq_conc", "best_match", "Functionality")]
 IDs$best_match = as.character(IDs$best_match)
@@ -32,7 +33,7 @@
 
 dat = data.frame(table(dat$seq_conc, dat$Functionality))
 
-dat = dat[dat$Freq > 1,]
+#dat = dat[dat$Freq > 1,]
 
 names(dat) = c("seq_conc", "Functionality", "Freq")
 
@@ -65,7 +66,7 @@
 	
 	classes.sum = sum(classes)
 	
-	if(!any(classes != 0 & classes != classes.sum)){
+	if(classes.sum == 1){
 		next
 	}