# HG changeset patch
# User davidvanzessen
# Date 1389343363 18000
# Node ID 141b0362603ec8c8bb813711734f6daab2ea81fe
# Parent  c134f2d0a484210c282315a3796c845845e63c0c
Uploaded

diff -r c134f2d0a484 -r 141b0362603e RScript.r
--- a/RScript.r	Wed Jan 08 06:26:03 2014 -0500
+++ b/RScript.r	Fri Jan 10 03:42:43 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)