Beefy Boxes and Bandwidth Generously Provided by pair Networks
No such thing as a small change
 
PerlMonks  

Re^2: calculate average in sliding windows

by Eliya (Vicar)
on May 23, 2012 at 14:27 UTC ( [id://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

Replies are listed 'Best First'.
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
Domain Nodelet?
Node Status?
node history
Node Type: note [id://972045]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this?Last hourOther CB clients
Other Users?
Others cooling their heels in the Monastery: (7)
As of 2025-06-18 07:57 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    Voting Booth?

    No recent polls found

    Notices?
    erzuuliAnonymous Monks are no longer allowed to use Super Search, due to an excessive use of this resource by robots.