Beefy Boxes and Bandwidth Generously Provided by pair Networks
Welcome to the Monastery

Re: Weighted random numbers generator

by blokhead (Monsignor)
on Mar 13, 2003 at 17:05 UTC ( #242751=note: print w/ replies, xml ) Need Help??

in reply to Weighted random numbers generator

I'm no probabalist, but your code seems on the ball to me. Here's an alternative method I use from time to time that doesn't require the creation of an additional accumulating-weight array:

sub choose_weighted { my $weights = shift; my $total = 0; $total += $_ for @$weights; my $rand_val = $total * rand; my $i = -1; $rand_val -= $weights->[++$i] while ($rand_val > 0); return $i; }

  1. Pick a random number less than the total
  2. At each position of the array, subtract the weight at that position
  3. When the result is negative, return that position

I have found this algorithm easier to adapt for multiple selections, i.e., calling choose_weighted( [1.75, 2, 3.6, 4], 2 ) and getting back two (different) numbers randomly chosen with this weighting. However, once you start doing this, you have to modify the array of weights (each time you select an item, set its weight to zero and subtract accordingly from $total), which may be undesired here.



Comment on Re: Weighted random numbers generator
Select or Download Code

Log In?

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

How do I use this? | Other CB clients
Other Users?
Others surveying the Monastery: (7)
As of 2015-11-25 15:20 GMT
Find Nodes?
    Voting Booth?

    What would be the most significant thing to happen if a rope (or wire) tied the Earth and the Moon together?

    Results (681 votes), past polls