I think your centering of mu_phi in phi2 = pm.Deterministic('phi', mu_phi - tt.mean(mu_phi)) makes the model unidentifiable.
I think your centering of mu_phi in phi2 = pm.Deterministic('phi', mu_phi - tt.mean(mu_phi)) makes the model unidentifiable.