Mercurial > repos > davidvanzessen > change_o

annotate define_clones.r @ 4:c75b30eb0b07 draft

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