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.