Mercurial > repos > davidvanzessen > combined_immune_repertoire_pipeline
changeset 14:b07346fe54b4 draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Fri, 10 Jan 2014 03:43:04 -0500 |
| parents | c22b8c7e40d3 |
| children | c69d9262c66d |
| files | RScript.r |
| diffstat | 1 files changed, 1 insertions(+), 1 deletions(-) [+] |
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)