Can I make out-of-sample predictions using posterior means? Or how to use scipy.optimize with pymc predictions?

[jessegrabowski]

It’s still the same point as above, that for a non-linear function f, E[f(x)] \neq f(E[x]).