Campionamento serbatoio

Question

Per recuperare i numeri casuali k da una matrice di dimensioni indeterminate si utilizza una tecnica chiamata campionamento del serbatoio. Qualcuno può evidenziare brevemente come accade con un codice di esempio?

Answer 1

Io in realtà non avevo capito c'era un nome per questo, così ho dimostrato e realizzato questo da zero:

import random 
def random_subset(iterator, K): 
    result = [] 
    N = 0 

    for item in iterator: 
     N += 1 
     if len(result) < K: 
      result.append(item) 
     else: 
      s = int(random.random() * N) 
      if s < K: 
       result[ s ] = item 

    return result 

Da: http://web.archive.org/web/20141026071430/http://propersubset.com:80/2010/04/choosing-random-elements.html

Con una prova vicino alla fine.

Answer 2

Java

import java.util.Random; 

public static void reservoir(String filename,String[] list) 
{ 
    File f = new File(filename); 
    BufferedReader b = new BufferedReader(new FileReader(f)); 

    String l; 
    int c = 0, r; 
    Random g = new Random(); 

    while((l = b.readLine()) != null) 
    { 
     if (c < list.length) 
      r = c++; 
     else 
      r = g.nextInt(++c); 

     if (r < list.length) 
      list[r] = l; 

     b.close();} 
} 

Answer 3

seguendo (1981) descrizione di Knuth più da vicino, Reservoir campionamento (Algoritmo R) potrebbe essere implementato come segue:

import random 

def sample(iterable, n): 
    """ 
    Returns @param n random items from @param iterable. 
    """ 
    reservoir = [] 
    for t, item in enumerate(iterable): 
     if t < n: 
      reservoir.append(item) 
     else: 
      m = random.randint(0,t) 
      if m < n: 
       reservoir[m] = item 
    return reservoir