5

に統計ブートストラップ ライブラリはありPythonますか?

で提供されているものと同様の機能が必要ですR bootstrap

http://statistics.ats.ucla.edu/stat/r/library/bootstrap.htm

私が見つけた検索:

http://mjtokelly.blogspot.com/2006/04/bootstrap-statistics-in-python.html (コードへのリンクが壊れています)

http://adorio-research.org/wordpress/?p=9048

https://github.com/cgevans/scikits-bootstrap

しかし、これらはすべての機能 (特に確率の重み) を提供しているようには見えません。

ポインタはありますか?

これは最近numpy.randomに追加されました

ありがとう

4

1 に答える 1

4

R のサンプル関数の Python バージョンを探しているだけなら、これを試してみてください。

import collections
import random
import bisect

def sample(xs, sample_size = None, replace=False, sample_probabilities = None):
    """Mimics the functionality of http://statistics.ats.ucla.edu/stat/r/library/bootstrap.htm sample()"""

    if not isinstance(xs, collections.Iterable):
        xs = range(xs)
    if not sample_size:
        sample_size = len(xs)            

    if not sample_probabilities:
        if replace:
            return [random.choice(xs) for _ in range(sample_size)]
        else:
            return random.sample(xs, sample_size)
    else:
        if replace:
            total, cdf = 0, []
            for x, p in zip(xs, sample_probabilities):
                total += p
                cdf.append(total)

            return [ xs[ bisect.bisect(cdf, random.uniform(0, total)) ] 
                    for _ in range(sample_size) ]
        else:            
            assert len(sample_probabilities) == len(xs)
            xps = list(zip(xs, sample_probabilities))           
            total = sum(sample_probabilities)
            result = []
            for _ in range(sample_size):
                # choose an item based on weights, and remove it from future iterations.
                # this is slow (N^2), a tree structure for xps would be better (NlogN)
                target = random.uniform(0, total)
                current_total = 0                
                for index, (x,p) in enumerate(xps):
                    current_total += p
                    if current_total > target:
                        xps.pop(index)
                        result.append(x)
                        total -= p
                        break
            return result
于 2013-03-04T00:50:05.883 に答える