Geneset=GO0035091 Size=77 ES=0.525 NES=2.913 NominalP=0.000 Geneset=GO0030163 Size=54 ES=0.463 NES=2.248 NominalP=0.013 Geneset=GO0007067 Size=44 ES=0.484 NES=1.975 NominalP=0.018 #### GO0046800 GO0046800 CD209 CD209L CD209L1 CD209_HUMAN CLC4M GO0032104 GO0032104 CART CARTPT CART_HUMAN GHRL GHRL_HUMAN ......... .......... ..... ...... .......... .... .......... ......... .......... ..... ...... .......... .... .......... GO0035091 GO0035091 41_HUMAN 5NTD_HUMAN 9804 A-388D4.1 ..... #### rs10904494 NP_817124 17881 rs7906287 NP_817124 39800 rs4881551 41_HUMAN 21567 rs5416721 5NTD_HUMAN 0 .................... .............. .... #### Marker CHI2 rs3749375 11.7268615355335 rs10499549 10.4656064706897 rs5416721 9.85374546064131 #### Geneset Genes SNP/Marker ES NES NominalP GO0035091 41_HUMAN rs4881551 .... ....... .... 5NTD_HUMAN rs5416721 .... ....... .... .................... .................... .................. #### open(OUTPUT, "; close(OUTPUT); open(GENESETS, "; close(GENESETS); @NewgeneSet; @genesInSet; $i=0; while($i < 10){ @outputLineSplit = split(/\t/,$output[$i]); #print "$outputLineSplit[0] \n"; $setName = $outputLineSplit[0]; $equalLoc = index($setName, "="); $setName = substr($setName,$equalLoc+1,length($setName)); #print "$setName\n"; @genesInSet[$i]= $setName; $i++; } foreach $genesInSet(@genesInSet){ print "$genesInSet\n"; foreach $geneSets(@geneSets){ if($geneSets=~m/$genesInSet/i){ #print "$geneSets\n"; } } }