# HG changeset patch
# User davidvanzessen
# Date 1448026033 18000
# Node ID ef3ac4f431bb6fc84df3f8ca40b2466579aac014
# Parent  7a63e90cfc3db3c75ad2895ed2c301f5a1f85266
Uploaded

diff -r 7a63e90cfc3d -r ef3ac4f431bb RScript.r
--- a/RScript.r	Mon Nov 09 08:19:53 2015 -0500
+++ b/RScript.r	Fri Nov 20 08:27:13 2015 -0500
@@ -17,7 +17,7 @@
 #require(xtable)
 cat("<tr><td>Reading input</td></tr>", file=logfile, append=T)
 dat = read.table(inFile, header=T, sep="\t", dec=".", fill=T, stringsAsFactors=F)
-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", "Related_to_leukemia_clone", "Clone_Sequence")]
+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$dsPerM = 0
 dat = dat[!is.na(dat$Patient),]
 dat$Related_to_leukemia_clone = F