Multiple priors on the same parameter?

6

Okay the below works:

theta = pm.Flat("theta")
pm.Potential("term1", pm.logp(pm.Normal.dist(2, 4), theta))
pm.Potential("term2", pm.logp(pm.Normal.dist(0, 10), theta))

now, what if we wanted to do a non-central parameterization? logp doesn’t accept additions so the following cannot be done:

theta = pm.Flat("theta")
pm.Potential("term1", pm.logp(2 + 4*pm.Normal.dist(0,1), theta))
pm.Potential("term2", pm.logp(10*pm.Normal.dist(0, 1), theta))

is there a way around this?