view define_clones.r @ 1:326165da9ece draft

Uploaded
author davidvanzessen
date Thu, 04 Aug 2016 08:59:47 -0400
parents
children 5173675d0d10
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)