Python rbeta Examples

Example #1

File: infinite_multinomials.py Project: JohnReid/infpy

def rBeta(alpha, beta):
    from rpy2.robjects import r
    x = r.rbeta(1, alpha, beta)
    if N.isnan(x):
        raise RuntimeError(
            'Cannot draw from beta distribution with parameters: (%f, %f)' % (alpha, beta))
    return x[0]

Example #2

File: test_scipy_beta.py Project: pombredanne/biopsy

import scipy.stats as S, numpy as N
from rpy2.robjects import r

alpha, beta = 0.0710, 0.4222
for i in xrange(20):
    x_from_scipy = S.beta.rvs(alpha, beta)
    x_from_R = r.rbeta(1, alpha, beta)
    print 'Alpha=%.2e; Beta=%.2e; scipy.stats.beta.rvs=%.2e; R.rbeta=%.2e' % (
        alpha, beta, x_from_scipy, x_from_R[0])
    alpha /= 10.
    beta /= 10.

Example #3

File: test_scipy_beta.py Project: JohnReid/biopsy

import scipy.stats as S, numpy as N
from rpy2.robjects import r

alpha, beta = 0.0710, 0.4222
for i in xrange(20):
    x_from_scipy = S.beta.rvs(alpha, beta)
    x_from_R = r.rbeta(1, alpha, beta)
    print "Alpha=%.2e; Beta=%.2e; scipy.stats.beta.rvs=%.2e; R.rbeta=%.2e" % (alpha, beta, x_from_scipy, x_from_R[0])
    alpha /= 10.0
    beta /= 10.0