# HG changeset patch # User davidvanzessen # Date 1471268126 14400 # Node ID c75b30eb0b07194160d807820878cebbf322af46 # Parent 5173675d0d10abf24248d2f8df0554cb2de7636f Uploaded diff -r 5173675d0d10 -r c75b30eb0b07 define_clones.r --- a/define_clones.r Mon Aug 15 08:33:53 2016 -0400 +++ b/define_clones.r Mon Aug 15 09:35:26 2016 -0400 @@ -7,8 +7,8 @@ freq = data.frame(table(change.o$CLONE)) freq2 = data.frame(table(freq$Freq)) -freq2$final = round(freq2$Freq * as.numeric(freq2$Var1), 2) +freq2$final = freq2$Freq * as.numeric(freq2$Var1) -names(freq2) = c("freq", "freq.of.freq", "divided") +names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences") write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)