Mercurial > repos > davidvanzessen > argalaxy_tools

annotate change_o/define_clones.r @ 4:5ffd52fc35c4 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Mon, 12 Dec 2016 05:22:37 -0500
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
1 args <- commandArgs(trailingOnly = TRUE)
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
2
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
3 input=args[1]
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
4 output=args[2]
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
5
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
6 change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F)
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
7
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
8 freq = data.frame(table(change.o$CLONE))
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
9 freq2 = data.frame(table(freq$Freq))
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
10
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
11 freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
12
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
13 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
14
5ffd52fc35c4 Uploaded
davidvanzessen
parents:
diff changeset
15 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)