Why not use a Gaussian process with the squared exponential cov function? The two is equivalent.
If you want to recreate the RBF regression, you can try rewriting the example from Edward into PyMC3: http://edwardlib.org/tutorials/supervised-classification with the rbf function.