# HG changeset patch
# User davidvanzessen
# Date 1389343206 18000
# Node ID 6ab4c67af3b7ffe8b537dfebb766b75857cdc4cd
# Parent  907e2a756f85d8f61f50a1946b85e64087127abe
Uploaded

diff -r 907e2a756f85 -r 6ab4c67af3b7 RScript.r
--- a/RScript.r	Wed Jan 08 06:26:43 2014 -0500
+++ b/RScript.r	Fri Jan 10 03:40:06 2014 -0500
@@ -289,7 +289,7 @@
 
 	clonalFreqCountSum = data.frame(data.table(clonalFreqCount)[, list(Numerator=sum(WeightedCount, na.rm=T)), by=c("Sample")])
 	clonalFreqCount = merge(clonalFreqCount, clonalFreqCountSum, by.x="Sample", by.y="Sample", all.x=T)
-
+	clonalFreqCount$ReadsSum = as.numeric(clonalFreqCount$ReadsSum) #prevent integer overflow
 	clonalFreqCount$Denominator = (((clonalFreqCount$ReadsSum * clonalFreqCount$ReadsSum) - clonalFreqCount$ReadsSquaredSum) / 2)
 	clonalFreqCount$Result = (clonalFreqCount$Numerator + 1) / (clonalFreqCount$Denominator + 1)