# HG changeset patch
# User davidvanzessen
# Date 1389343384 18000
# Node ID b07346fe54b45514c00e1c42b5eceda933276dc2
# Parent  c22b8c7e40d39ac47bab4d219b307b524c177840
Uploaded

diff -r c22b8c7e40d3 -r b07346fe54b4 RScript.r
--- a/RScript.r	Wed Jan 08 06:30:36 2014 -0500
+++ b/RScript.r	Fri Jan 10 03:43:04 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)