#!/usr/bin/perl use strict; use warnings; # prints out inhibited object, directly activated ohjects, directly deactivated objects, indirectly activated objects, indirectly deacivated objects print "Modified Object\tActivated Objects\tDeactivated Objects\n"; open(DAT,'<',"./useful_dat.tab")||die "Canot open file\n"; my @file=; close DAT; shift @file;##remove header my %true; foreach my $line(@file){ my @line=split(/\t/,$line); if($true{$line[0]}){ next;#if element has already come up, skip it } else{ $true{$line[0]}=1;#if object has already come up, store true } my @activated; my @inhibited; print "$line[0]\t";#inhibited object &inhibited_rec($line[0]); foreach(@activated){ print "$_,"; } print "\t"; foreach(@inhibited){ print "$_,"; } print "\n"; } sub inhibited_rec{ my @a=`awk '\$1 ~ /$_[0]/' useful_dat.tab |cut -f 2,3|sort|uniq`; foreach (@a){ my @b=split(/\t/,$_); if($b[1] eq 'Activation'){ } else{ push (@activated,$b[0]); &activated_rec($b[0]); } } } sub activated_rec{ my @a=`awk '\$1 ~ /$_[0]/' useful_dat.tab |cut -f 2,3|sort|uniq`: foreach (@a){ my @b=split(/\t/,$_); if($b[1] eq 'Activation'){ push (@activated,$b[0]); &activated_rec($b[0]); } else{ push (@inhibited,$b[0]); &inhibited_rec($b[0]); } }