Mercurial > repos > davidvanzessen > combined_immune_repertoire_imgt
diff RScript.r @ 8:141b0362603e draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Fri, 10 Jan 2014 03:42:43 -0500 |
| parents | c134f2d0a484 |
| children |
line wrap: on
line diff
--- 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)