Beefy Boxes and Bandwidth Generously Provided by pair Networks
Don't ask to ask, just ask
 
PerlMonks  

Re^2: calculate average in sliding windows

by Eliya (Vicar)
on May 23, 2012 at 14:27 UTC ( #972045=note: print w/ replies, xml ) Need Help??


in reply to Re: calculate average in sliding windows
in thread calculate average in sliding windows

but it doesn't take into account the positions in column 1

Sorry, I'm still not sure what you want (and it's a little diffficult to infer it from code which - as you state - doesn't work).

Are you saying the positions in column 1 are meant to represent indices in some sparse array, with all unspecified positions (such as 12498249..12512574, etc.) implicitly having a value of zero, or what?

I think it would help if you could come up with some simplified (but sufficient/complete) sample input (say using a window size of 3) together with the desired output values.


Comment on Re^2: calculate average in sliding windows
Re^3: calculate average in sliding windows
by marfabe (Initiate) on May 23, 2012 at 14:43 UTC

    Thanks Eliya
    Here's the simplified sample set

    position value 1 1 10 3 30 1 40 2 60 2

    And here the output for window size of 30

    1-30 1.666666667 31-60 2

      Ok, thanks.  So maybe like this?

      #!/usr/bin/perl -lw use strict; use constant WINDOW => 30; my $sum = 0; my $n = 0; my $p = WINDOW; while (<DATA>) { my ($pos, $val) = split; if (eof) { # corner case $sum += $val; $n++; } if ($pos > $p or eof) { print $sum / $n if $n > 0; $sum = 0; $n = 0; $p += WINDOW while $pos > $p; } $sum += $val; $n++; } __DATA__ 1 1 10 3 30 1 40 2 60 2

      Output:

      1.66666666666667 2

      (Upd: fixed handling of corner case)

        Yes!! It works perfectly! Thank you!

Log In?
Username:
Password:

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

How do I use this? | Other CB clients
Other Users?
Others about the Monastery: (11)
As of 2015-07-02 22:07 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    The top three priorities of my open tasks are (in descending order of likelihood to be worked on) ...









    Results (46 votes), past polls