# HG changeset patch # User davidvanzessen # Date 1464605353 14400 # Node ID fd63f5307962175f6c989dfdcb367f1bdec3201e # Parent 0011a05976854ee636c38fa0adc00a2587db967e Uploaded diff -r 0011a0597685 -r fd63f5307962 sequence_overview.r --- 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 }