# HG changeset patch
# User davidvanzessen
# Date 1453478041 18000
# Node ID 2e0a7c35082ef86c77518c2b9046c90a0182b9ab
# Parent  539618b939bc1970350f79d3dcd1b428db56b6ce
Uploaded

diff -r 539618b939bc -r 2e0a7c35082e report_clonality/RScript.r
--- a/report_clonality/RScript.r	Fri Jan 22 10:44:50 2016 -0500
+++ b/report_clonality/RScript.r	Fri Jan 22 10:54:01 2016 -0500
@@ -564,7 +564,8 @@
     ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "clonaltype")])
     ReplicateReads = data.frame(data.table(ReplicateReads)[, list(Reads=.N), by=c("Sample", "Replicate")])
     clonalFreqCount$Reads = as.numeric(clonalFreqCount$Reads)
-    ReplicateReads$squared = ReplicateReads$Reads * ReplicateReads$Reads
+    ReplicateReads$Reads = as.numeric(ReplicateReads$Reads)
+    ReplicateReads$squared = as.numeric(ReplicateReads$Reads * ReplicateReads$Reads)
     
     ReplicatePrint <- function(dat){
       write.table(dat[-1], paste("ReplicateReads_", unique(dat[1])[1,1] , ".csv", sep=""), sep=",",quote=F,na="-",row.names=F,col.names=F)