Mercurial > repos > rhpvorderman > shm_csr

view change_o/define_clones.r @ 5:495a521cf9f2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload commit 40d62ec3d2fea3c3bac605c976941e1f3b7e2cd4"
author rhpvorderman
date Tue, 16 Nov 2021 11:16:06 +0000
parents 64d74ba01a7c
children
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 = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1))

names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences")

write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T)