Adding constraint to the parameter

[junpenglao]

@gaddamanil16

I want the constraint to be p1+p2+p3 = 1. How can I do it?

What you are describing is a simplex, which is usually modelled using a Dirichlet distribution:

p = pm.Dirichlet('p', a=np.ones(3))