|
1
|
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))
|
|
5
|
10
|
|
|
11 print(freq)
|
|
|
12
|
|
|
13 print(freq2$Freq)
|
|
|
14
|
|
|
15 print(freq2$Var1)
|
|
|
16
|
|
|
17 freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))
|
|
|
18
|
|
|
19 print(freq2)
|
|
1
|
20
|
|
4
|
21 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")
|
|
1
|
22
|
|
|
23 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)
|
