Mercurial > repos > davidvanzessen > change_o
view define_clones.r @ 3:5173675d0d10 draft
Uploaded
author | davidvanzessen |
---|---|
date | Mon, 15 Aug 2016 08:33:53 -0400 |
parents | 326165da9ece |
children | c75b30eb0b07 |
line wrap: on
line source
args <- commandArgs(trailingOnly = TRUE) input=args[1] output=args[2] change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F) freq = data.frame(table(change.o$CLONE)) freq2 = data.frame(table(freq$Freq)) freq2$final = round(freq2$Freq * as.numeric(freq2$Var1), 2) names(freq2) = c("freq", "freq.of.freq", "divided") write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)