fonnesbeck
10/26/2016 - 9:38 PM

Sampling from the product of two betas

Sampling from the product of two betas

from pymc3 import Beta, Deterministic, Model, sample

with Model():
    a = Beta('a', 3, 1)
    b = Beta('b', 2, 5)
    a_times_b = Deterministic('a_times_b', a*b)

    trace = sample(2000)