Mercurial > repos > davidvanzessen > combined_immune_repertoire_pipeline

diff RScript.r @ 14:b07346fe54b4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Fri, 10 Jan 2014 03:43:04 -0500
parents a15d060a6248
children c69d9262c66d
line wrap: on
line diff
--- 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)