## Solving stochastic differentials for time series forecasting

I am trying to reproduce results of a paper. The model is:

$dS = (v-y-\lambda_1)Sdt + \sigma_1Sdz_1 \\ dy = (-\kappa y - \lambda_2)dt + \sigma_2 dz_2 \\ dv = a((\bar{v}-v)-\lambda_3)dt + \sigma_3 dz_3 \\ dz_1dz_2 = \rho_{12}dt \\ dz_1dz_3 = \rho_{13}dt \\ dz_2dz_3 = \rho_{23}dt \\$
where S, y, and v are state variables and $z_1, z_2, z_3$ are Wiener processes.

The parameters $\kappa, a, \bar{v}, \lambda_1, \lambda_2, \lambda_3, \sigma_1, \sigma_2, \sigma_3, \rho_{12}, \rho_{13}, \rho_{23}$ are given. There is time series data such for each day I can estimate S, y, and v.

The goal is to estimate the next day's S, y, and v. It feels like an easy task, but am not familiar with differential equations or Wiener processes (let alone stochastic calculus). Without understanding the model, I was naively trying to "substitute" dt with Δt, etc, but the given parameters are such that $\rho_{12} > 0, \, \rho_{23} > 0, \, \rho_{13} < 0$, so I realized I should not do that.

Any help would be greatly appreciated.
 PhysOrg.com science news on PhysOrg.com >> New language discovery reveals linguistic insights>> US official: Solar plane to help ground energy use (Update)>> Four microphones, computer algorithm enough to produce 3-D model of simple, convex room
 Recognitions: Homework Help Are you trying to solve the equations numerically? I think people typically interpret the dt as a small timestep Δt. I guess the tricky part is the stochastic terms - it looks like they should be Gaussian distributed random variables with covariances $\rho_{12}$, etc? Generating those might be tricky? Otherwise everything else looks straightforward. The equations are simple enough to solve analytically, but if you're not familiar with differential equations that might not be the best avenue for you. But, if you want to try, you can solve the y and v equations using integrating factors, and then plug the solutions into the S equation, which you can solve by dividing out the S and writing $dS/S = d(\ln S)$, then it's a simple matter of integration. (I would write $dz_1/dt = \eta_1(t)$, etc.)

 Tags stochastic calculus, time series, wiener process

 Similar discussions for: Solving stochastic differentials for time series forecasting Thread Forum Replies Differential Equations 0 Calculus 4 Set Theory, Logic, Probability, Statistics 2 Set Theory, Logic, Probability, Statistics 2