comparison define_clones.r @ 3:5173675d0d10 draft

Uploaded
author davidvanzessen
date Mon, 15 Aug 2016 08:33:53 -0400
parents 326165da9ece
children c75b30eb0b07
comparison
equal deleted inserted replaced
2:6f57d43a1f4c 3:5173675d0d10
5 5
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 freq2$final = round(freq2$Freq / as.numeric(freq2$Var1), 2) 10 freq2$final = round(freq2$Freq * as.numeric(freq2$Var1), 2)
11 11
12 names(freq2) = c("freq", "freq.of.freq", "divided") 12 names(freq2) = c("freq", "freq.of.freq", "divided")
13 13
14 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T) 14 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)