Beefy Boxes and Bandwidth Generously Provided by pair Networks
Syntactic Confectionery Delight
 
PerlMonks  

Re: Discard if present within other coordinate range

by kennethk (Abbot)
on Feb 17, 2017 at 22:58 UTC ( [id://1182246]=note: print w/replies, xml ) Need Help??


in reply to Discard if present within other coordinate range

What have you tried? What worked? What didn't? See How do I post a question effectively? In your previous posts, you've always put up code.

Does order in the output matter? The easiest solution I can see would be to sort your entries first by the first number ascending, and then by the second number descending. You can then cycle through the list, and any time the current second number is less than the biggest one you've seen so far, you know you've gotten a subset.

#!/usr/bin/perl use strict; use warnings; use 5.10.0; my @terms; while (<DATA>) { chomp; push @terms, [split /\s+/]; } my $biggest = 0; for my $term (sort sorter @terms) { if ($term->[2] > $biggest) { say join ' ', @$term; $biggest = $term->[2]; } } sub sorter { $a->[1] <=> $b->[1] || $b->[2] <=> $a->[2] } __DATA__ SEQ1 225 275 SEQ1 200 300 SEQ1 201 299 SEQ1 250 399 SEQ1 145 244 SEQ2 120 130 SEQ2 100 150 SEQ2 101 149 SEQ2 120 230 SEQ2 99 140

#11929 First ask yourself `How would I do this without a computer?' Then have the computer do it the same way.

Replies are listed 'Best First'.
Re^2: Discard if present within other coordinate range
by Cristoforo (Curate) on Feb 18, 2017 at 20:55 UTC
    I had a different take (that makes minor changes to your solution). I thought he wanted non-subsets within each seq identifier, (where your solution finds them for the whole file not considering the first column).

    I added 2 more lines to the for loop and 1 more line to the sorter sub.

    #!/usr/bin/perl use strict; use warnings; use 5.010; open my $fh, '<', \<<EOF; SEQ1 225 275 SEQ1 200 300 SEQ1 201 299 SEQ1 250 399 SEQ1 145 244 SEQ2 120 130 SEQ2 100 150 SEQ2 101 149 SEQ2 120 230 SEQ2 99 140 EOF my @terms; while (<$fh>) { chomp; push @terms, [split /\s+/]; } my $biggest = 0; my $id = ''; for my $term (sort sorter @terms) { $biggest = 0 if $id ne $term->[0]; if ($term->[2] > $biggest) { say join ' ', @$term; $biggest = $term->[2]; } $id = $term->[0]; } sub sorter { $a->[0] cmp $b->[0] || $a->[1] <=> $b->[1] || $b->[2] <=> $a->[2] }
    The output is:
    SEQ1 145 244 SEQ1 200 300 SEQ1 250 399 SEQ2 99 140 SEQ2 100 150 SEQ2 120 230
Re^2: Discard if present within other coordinate range
by Cow1337killr (Monk) on Feb 20, 2017 at 04:55 UTC

Log In?
Username:
Password:

What's my password?
Create A New User
Domain Nodelet?
Node Status?
node history
Node Type: note [id://1182246]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others having an uproarious good time at the Monastery: (7)
As of 2024-03-29 00:07 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found