Using a bounded variable within a Mixture

Questions
5

@aloctavodia
Yes, that second one works. Thank you very much.

However, how do I refer to the variable by name? So I can grab it from trace, like so trace[‘myvar’]. Here is the code in case someone else needs to use similar tricks:

with pm.Model() as model:
    alpha = pm.Gamma('alpha', 1,1)
    beta = pm.Beta('beta', 1, alpha, shape=K)
    w = pm.Deterministic('w', stick_breaking(beta))
    #a = pm.Uniform('a', 0, 50, shape=K)
    #b = pm.Uniform('b', 0, 50, shape=K)
    mu = pm.Uniform('mu', 0, 50, shape=K)
    bp = [pm.Bound(pm.Poisson, lower=1).dist( mu=mu[i]) for i in range(K)] 
    #x = pm.Mixture('obs', w, pm.Poisson.dist(mu), observed=df['x1'].values)
    x = pm.Mixture('obs', w, bp, observed=df['x1'].values)