comparison change_o/define_clones.r @ 0:c33d93683a09 draft

Uploaded
author davidvanzessen
date Thu, 13 Oct 2016 10:52:24 -0400
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:c33d93683a09
1 args <- commandArgs(trailingOnly = TRUE)
2
3 input=args[1]
4 output=args[2]
5
6 change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F)
7
8 freq = data.frame(table(change.o$CLONE))
9 freq2 = data.frame(table(freq$Freq))
10
11 freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))
12
13 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")
14
15 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)