Beefy Boxes and Bandwidth Generously Provided by pair Networks
XP is just a number
 
PerlMonks  

Re^3: Got some problem with read write file

by choroba (Cardinal)
on Sep 03, 2016 at 21:03 UTC ( [id://1171123]=note: print w/replies, xml ) Need Help??


in reply to Re^2: Got some problem with read write file
in thread Got some problem with read write file

Why is the second while there? Do you want to find several occurences on the same line? Also, you probably don't want to print all the genes found so far after finding a gene, you want to print them once all of them have been found:
#!/usr/bin/perl use warnings; use strict; open my $HUBFILE, '<', '1048_undefined.tsv' or die $!; my @hubs; while (my $line = <$HUBFILE>) { push @hubs, $1 if $line =~ /\d \t (\w+) \t/x; } close $HUBFILE; open my $OUT, '>', 'hubs.txt' or die $!; for my $hub (@hubs) { print {$OUT} "HUB:$hub\n"; } close $OUT;

Notice I modified some other parts of the code, too: I switched to 3-argument open with lexical filehandles, foreach style loop instead of the C-style one, etc.

($q=q:Sq=~/;[c](.)(.)/;chr(-||-|5+lengthSq)`"S|oS2"`map{chr |+ord }map{substrSq`S_+|`|}3E|-|`7**2-3:)=~y+S|`+$1,++print+eval$q,q,a,

Replies are listed 'Best First'.
Re^4: Got some problem with read write file
by SilverWol (Novice) on Sep 03, 2016 at 21:15 UTC

    Thought wrong about the second while, I want to print in the second file a list with all the gene names but it just print only the 1st gene name. I didn't know about the lexical filehandles, thank you!

      I am not a bioinformatician, but I wanted to learn about bioinformatics. I read Beginning Perl for Bioinformatics. I really enjoyed it. It got me started on understanding bioinformatics. I think you will enjoy it. It explains how to code in Perl by just explaining to biologists: Here is the problem, here is the code to solve the problem, let me explain the code, move on to the next example program.

      (I got my copy from the library.)

Log In?
Username:
Password:

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

How do I use this?Last hourOther CB clients
Other Users?
Others meditating upon the Monastery: (4)
As of 2024-04-25 17:07 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found