Mercurial > repos > davidvanzessen > prisca

diff RScript.r @ 7:7ce82833977c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 12 Dec 2017 04:53:20 -0500
parents a9d2ed661541
children eb2aa7cffca3
line wrap: on
line diff
--- a/RScript.r	Thu Oct 05 09:06:21 2017 -0400
+++ b/RScript.r	Tue Dec 12 04:53:20 2017 -0500
@@ -34,6 +34,7 @@
 }
 
 dat = dat[,c("Patient",  "Receptor", "Sample", "Cell_Count", "Clone_Molecule_Count_From_Spikes", "Log10_Frequency", "Total_Read_Count", "J_Segment_Major_Gene", "V_Segment_Major_Gene", "CDR3_Sense_Sequence", "Clone_Sequence")] 
+dat = dat[!(nchar(dat$Clone_Sequence) < 2),]
 
 dat$dsPerM = 0
 dat = dat[!is.na(dat$Patient),]