I have some data r and want to run this Multinomial model, I get an initialisation error, although the prior predictive works. T = 12, g = 3, I = 50 with pm.Model() as model: alpha = pm.Gamma('alpha', 1, 1, shape=g) M = pm.Gamma('M', 1, 1) K = pm.Gamma('K', 1, 1) mu = pm.Gamma('mu', M * T,…

I think the problem is likely: n = pm.Poisson('n', mu, shape=I) Since you observed the multinomial, n is known n = r.sum(axis=-1). I guess you are assigning it as random variable so you can infer the parameter mu and its hyper priors. So suggestion here is to do n = pm.Poisson('n', mu, observed=r…

MultiNomial init error

Questions

Dirk_Nachbar1 February 8, 2022, 8:00am 4

Junpeng, you are right, it dawned on me later, n was actually given.

Topic		Replies	Views
pm.Multinomial gives Sampling Error v5 bug	3	462	October 20, 2022
Multinomial sampling init val errors out v5	5	348	January 18, 2023
Multinomial sampling crashes Questions	2	501	December 23, 2020
Multinomial hierarchical regression with multiple observations per group ("Bad energy issue") Questions	22	3450	March 14, 2019
Dirichlet-Multinomial Model does not appear to be sampling at all Questions	8	2697	November 9, 2018

MultiNomial init error

Related topics