The following code executes without problems. import pymc3 as pm alpha1_, beta1_ = pm.Exponential.dist(0.1).random(size=2) alpha2_, beta2_ = pm.Exponential.dist(0.1).random(size=2) w_ = pm.Dirichlet.dist(a=np.ones(2), shape=2).random(size=1) p0 = pm.Gamma.dist([alpha1_, alpha2_], [beta1_, beta2_], …

Yeah this is a bug as the mixture class is wrongly rely on the mode to exist - could you file a bug on Github?

Mixture of Beta fails (because no mode exists) while Gamma is ok

junpenglao February 9, 2020, 8:24am 5

Thanks for reporting back - I think the issue could be leave open as is.

Topic		Replies	Views
Mixture of two beta geometric with beta shape > 1 Questions	6	1198	January 18, 2018
Zero and One inflated Beta Distribution Questions	5	1380	July 8, 2020
Struggling with Beta mixture models Questions	3	1874	February 2, 2019
Could you give me some advice about this modle? Questions	1	691	August 9, 2018
Is there an example on how to work with generalized mixture models? Questions	15	3593	March 16, 2019