Mercurial > repos > davidvanzessen > plotting_merged
changeset 53:6ab4c67af3b7 draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Fri, 10 Jan 2014 03:40:06 -0500 |
| parents | 907e2a756f85 |
| children | 3fb05bdd5299 |
| files | RScript.r |
| diffstat | 1 files changed, 1 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/RScript.r Wed Jan 08 06:26:43 2014 -0500 +++ b/RScript.r Fri Jan 10 03:40:06 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)