Hi! I have a quite hard question that have haunting me for months. I’m building a complicated model which looks like this: [image] In the red circle, I make the random variables “ARR” and “LARR” equal by soft zero difference: zerodiff = pm.Normal.dist(mu=0, sigma=0.01) pm.Potential(f"z…

I guess I would try to avoid defining a model that way. If you want y to be a function of X, link them directly via regression. Something like y = pm.Normal('y', mu = hyper_y_mu + beta_y * X, sigma = hyper_y_sigma). My advice would be that any time you run into a situation where you want to impose …

Make two random variables equal

Questions v5

jessegrabowski June 20, 2023, 4:14am 6

Can you use a deterministic for LARR then? Keep ARR as-is, and use LARR = pm.Deterministic(ARR[year_indices])?

Topic		Replies	Views
Notation for hierarchical models Questions	5	1061	April 20, 2018
Constrain one variable to be greater than another Questions	14	3519	January 30, 2018
MAP not working with shape>1 and VI on a pooled model Questions	11	837	November 25, 2017
Hierarchical linear model with two parent distributions Questions	16	1338	September 15, 2018
Imposing order constraint Questions	3	976	April 9, 2019

Make two random variables equal

Related topics