Hello,
This is a basic mixture example : 2 beta components to estimate the rate of an observed binomial. It samples well, and the posterior of the rate (r) looks ok.
Now I’d like to see the individual component posterior (c), but they are not in the trace.
How can I see or produce them ?
with pm.Model() as model:
n = 2
c = pm.Beta.dist(alpha=1,beta=1,shape=n)
w = [1/2,1/2]
r = pm.Mixture('r',w=w,comp_dists=c)
obs = pm.Binomial('obs',n=[10000],p=r,observed=[5000])
idata = pm.sample()
az.plot_trace(idata)
Any help will be appreciated.