Mercurial > repos > davidvanzessen > change_o

comparison define_clones.r @ 6:3ddd933dd7a2 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author davidvanzessen
date Thu, 15 Sep 2016 03:54:33 -0400
parents 3d97839fda21
children
comparison
equal deleted inserted replaced
5:3d97839fda21 6:3ddd933dd7a2
6 change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F) 6 change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F)
7 7
8 freq = data.frame(table(change.o$CLONE)) 8 freq = data.frame(table(change.o$CLONE))
9 freq2 = data.frame(table(freq$Freq)) 9 freq2 = data.frame(table(freq$Freq))
10 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)) 11 freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))
18
19 print(freq2)
20 12
21 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences") 13 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")
22 14
23 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T) 15 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)