--- a/define_clones.r	Mon Aug 15 09:35:26 2016 -0400
+++ b/define_clones.r	Mon Aug 15 10:06:24 2016 -0400
@@ -7,7 +7,16 @@
 
 freq = data.frame(table(change.o$CLONE))
 freq2 = data.frame(table(freq$Freq))
-freq2$final = freq2$Freq * as.numeric(freq2$Var1)
+
+print(freq)
+
+print(freq2$Freq)
+
+print(freq2$Var1)
+
+freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))
+
+print(freq2)
 
 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")