Setting prior over initial values for Euler-Mayurama model

Questions
#1

Hi, is there a way to define a prior over initial values of an SDE used in Euler-Mayurama distribution?

For example I have defined the following SDE

def ou_sde(x, theta, alpha, s):
    "Ornstein-Uhlenbeck process"
    return alpha - theta * x, s

and passed it to the Euler-Mayurama distribution from timeseries models.

However, I can not figure out how to pass initial conditions for the integration process.
For example,

with pm.Model() as model:
    x0 = pm.Normal('x', 0., 1.)
    x = EulerMayurama('x', ou_sde, *args, init=x0)

The reason I am asking this is that by default both prior and posterior initial values make no sense, as I and up with some very large values for x0, far away from the mean of the process.