# HG changeset patch # User davidvanzessen # Date 1471269984 14400 # Node ID 3d97839fda21128bda04189e260a558dc6c52eaa # Parent c75b30eb0b07194160d807820878cebbf322af46 Uploaded diff -r c75b30eb0b07 -r 3d97839fda21 define_clones.r --- 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")