Mercurial > repos > davidvanzessen > combined_immune_repertoire_pipeline

diff RScript.r @ 15:c69d9262c66d draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Tue, 14 Jan 2014 04:33:44 -0500
parents b07346fe54b4
children
line wrap: on
line diff
--- a/RScript.r	Fri Jan 10 03:43:04 2014 -0500
+++ b/RScript.r	Tue Jan 14 04:33:44 2014 -0500
@@ -267,6 +267,7 @@
 
 	ReplicateReads = data.frame(data.table(clonalityFrame)[, list(Type=.N), by=c("Sample", "Replicate", "VDJCDR3")])
 	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
 
 	ReplicatePrint <- function(dat){