Beefy Boxes and Bandwidth Generously Provided by pair Networks
"be consistent"
 
PerlMonks  

Re: stupid/simple mistake

by TomDLux (Vicar)
on Oct 18, 2011 at 15:55 UTC ( #932185=note: print w/replies, xml ) Need Help??


in reply to stupid/simple mistake

#!/usr/bin/perl use warnings; use strict; use autodie; use List::Util; my $file = 'sequence.txt'; my $goal = 'GATC'; my $stats = process_file( $file, $goal ); # do something with stats ... # ---------------------------------------- # Subroutines # sub process_file { my ( $file, $goal ) = @_; open my $infh, '<', $file; my $stats = process_lines( $infh, $goal ); close $infh; return $stats } sub process_lines { my ( $infh, $goal ) = @_; my @stats; while ( my $line - <$infh> ) { chomp $line; my $linestats = process_one_line( $line, $goal ); push @stats, $linestats || 0; } return \@stats; } sub process_one_line { my ( $line, $goal ) = @_ my @occurences; my ( $offset ) = ( 0 ); SEEK: while ( 1 ) { $idx = index( $line, $goal, $offset ); last SEEK if $idx == -1; # no more occurences push @occurences, $idx; $offset = $idx; } return calc_avg_distance( \@occurences, length $goal ); } sub calc_avg_distance { my ( $occurences, $len ) = @_; return unless $occurences and scalar @$occurences; my $start = shift @$occurences; my @distances; while ( my $end = shift @$occurences ) { push @distances, ( $end - $start ) - $len; $start = $end; } my $sum = reduce { $a + $b }, @distances; my $n = scalar @distances; return $sum / $n; }

As Occam said: Entia non sunt multiplicanda praeter necessitatem.

Log In?
Username:
Password:

What's my password?
Create A New User
Node Status?
node history
Node Type: note [id://932185]
help
Chatterbox?
[Corion]: Oh yay. $project has faffed around for 2 years and now pushed a hasty 2 weeks solution into production without telling me. Of course, the new data also needs changes on my side to be processed correctly. Now they try to raise incidents against my ...
[Corion]: ... input systems to make the change outside of the organized project because I don't have ressources for the project to make the change.

How do I use this? | Other CB clients
Other Users?
Others lurking in the Monastery: (6)
As of 2017-08-17 09:12 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?
    Who is your favorite scientist and why?



























    Results (285 votes). Check out past polls.

    Notices?