I'm guessing you are looking for the genes for a certain id, and then a looking for all the id's that have those genes. If so try this
#!perl
use strict;
use warnings;
my $match = 'DOID:2055';
my $filename = 'do_human_mapping.gmt';
open (my $fh, '<', $filename) or die "open: $!";
my @genes=();
my %gene2id=();
while (<$fh>){
my ($name,$id,@temp) = split /\s+/;
if ($id eq $match) {
@genes = @temp;
} else {
for my $gene (@temp){
push @{$gene2id{$gene}},$id
}
}
}
for my $gene (@genes){
if (exists $gene2id{$gene}){
print join ' ',$gene,@{$gene2id{$gene}},"\n";
}
}
__DATA__
DOID:00001 APOE IL4 RTG5
DOID:00002 FG6 CRH APOE
DOID:00003 RTG5 HUTN CRH
DOID:2055 APOE FKBP5 CRH
poj
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.