# HG changeset patch
# User davidvanzessen
# Date 1389692049 18000
# Node ID 3fb05bdd5299b86be0f2afe10ff57467d668c089
# Parent  6ab4c67af3b7ffe8b537dfebb766b75857cdc4cd
Uploaded

diff -r 6ab4c67af3b7 -r 3fb05bdd5299 RScript.r
--- a/RScript.r	Fri Jan 10 03:40:06 2014 -0500
+++ b/RScript.r	Tue Jan 14 04:34:09 2014 -0500
@@ -267,6 +267,7 @@
 
 	ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "VDJCDR3")])
 	ReplicateReads = data.frame(data.table(ReplicateReads)[, list(Reads=.N), by=c("Sample", "Replicate")])
+	clonalFreqCount$Reads = as.numeric(clonalFreqCount$Reads)
 	ReplicateReads$squared = ReplicateReads$Reads * ReplicateReads$Reads
 
 	ReplicatePrint <- function(dat){