This is a test version of Biostars. For the public version, visit https://www.biostars.org.
select the highest value for the same gene

Hi every body,

I have the coverage for different genes as:

GENE COVERAGE
A     0.7
A     0.2
A     0.9
B     0.5
B     1.2
B     0.3
B     0.5
B     0.6
C     0.1

and I want to get ONLY the highest coverage for each gene as follow

GENE COVERAGE
  A     0.9
  B     1.2
  C     0.1

Because I need the most representative to calculate the density of each one. Any suggestion?

PS. I can't install mySQL in the computer, I know that it will be the best option...

Thanks!

gene chip-seq
• 1,815 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by Lila M Lila M

1 answer

assuming that tabulation is the delimiter:

sed 's/^GENE/#GENE/' input.txt | LC_ALL=C sort -t $'\t' -k1,1 -k2,2rg  | LC_ALL=C sort -t $'\t' -k1,1 --stable -u
  • the first sed is just here to be sure that the header will be the first row in the final output
  • sort on first column , and then second column, numeric, reverse order
  • sort on first column, keep unique, don't mess the input order
 #GENE  COVERAGE
 A  0.9
 B  1.2
 C  0.1
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

can you explain a bit the code please?

Thank you!

• 0 views • link
written by Lila M Lila M

I am not very good using this language, it seems quite easy to use, for example, if I want to order based on different colunm for example 7 and 11 I can try: 

sed 's/^GENE/#GENE/' gene_body_annot_coverage | LC_ALL=C sort -t $'\t ' -k7,7 -k11,11rg  | LC_ALL=C sort -t $'\t' -k7,7 --stable -u > output

But can you explain what LC_ALL=C means?

Thanks!

• 0 views • link
written by Lila M Lila M

this 'language' is bash.

The column are specified by the '-k' option. So it would be (....) -k7,7 -k11,11rg (...)

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

yes, sorry about that, I've just figure it out very quickly ! Nice approach to handling files :)

• 0 views • link
written by Lila M Lila M

Log in to answer this question.

More posts like this