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

Guys, thank you for links. It’s a bit hard for me to understand everything that happening there, nevertheless I will spend more time this week. Seems like @ricardoV94 explained my idea clearly in his last message (this is exactly what I thought I needed).

@jessegrabowski can you please share an example of non-linear function, just to understand the difference compare to linear?