This is a test version of Biostars. For the public version, visit https://www.biostars.org.
grouping the genes based on the length of coding sequence

I have a list of genes with coding sequence length. I want to group them into 10 groups. do you know what the best way is to do so?

sequencing
• 1,497 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by ashkan ashkan

1 answer

if your data is:

gene1(tab)length
gene2(tab)length
(...)
geneN(tab)length

I would try:

sort -t $'\t' -k2,2n data.tsv | split -l ` awk 'END {print  NR/10}' data.tsv `  - OUT
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Pierre always give us shell magic :)

• 0 views • link
written by shenwei356 shenwei356

Log in to answer this question.

More posts like this