Mercurial > repos > davidvanzessen > shm_csr
comparison change_o/define_clones.r @ 81:b6f9a640e098 draft
Uploaded
| author | davidvanzessen |
|---|---|
| date | Fri, 19 Feb 2021 15:10:54 +0000 |
| parents | c33d93683a09 |
| children |
comparison
equal
deleted
inserted
replaced
| 80:a4617f1d1d89 | 81:b6f9a640e098 |
|---|---|
| 1 args <- commandArgs(trailingOnly = TRUE) | |
| 2 | |
| 3 input=args[1] | |
| 4 output=args[2] | |
| 5 | |
| 6 change.o = read.table(input, header=T, sep="\t", quote="", stringsAsFactors=F) | |
| 7 | |
| 8 freq = data.frame(table(change.o$CLONE)) | |
| 9 freq2 = data.frame(table(freq$Freq)) | |
| 10 | |
| 11 freq2$final = as.numeric(freq2$Freq) * as.numeric(as.character(freq2$Var1)) | |
| 12 | |
| 13 names(freq2) = c("Clone size", "Nr of clones", "Nr of sequences") | |
| 14 | |
| 15 write.table(x=freq2, file=output, sep="\t",quote=F,row.names=F,col.names=T) |