In case you need any further examples I bring up an explicit problem. I included one above in a reply btw hopefully that adds some clarification. For clarify though here's the explicit form of the exact question I'm seeking an answer to:

I) Solve the following initial-value difference equations. Plot the time path of Yt and discuss the convergence or divergence of the solutions as $ t \rightarrow \infty $

$$ y_t = 1.6 + 0.75 y_{t-1} + \epsilon_t $$ $$y_0 = 2 $$

At first it was an issue with the error term. But now that you showed me the syntax for that -- all that would remain is how to specify the initial conditions correctly. What's odd is I'm not getting an error and your code is directly from the references so I can't see it being incorrect. But the solution it gives me when I use your code snippet on these problems with the initial conditions specified I get an unfamiliar answer as you see above.

Hope this clears some things up.

Thanks again.

Ali,

Glad you got it working again by reopening MMA.

As to the other part of your post, I'm not sure I understand what you are trying to do. You say

But as far as then providing an functional solution to the associated recurrence equation, I haven't had any luck

I thought you were looking for a way to enter the recurrence equation as an ARMAProcess and compute stationarity (with WeakStationarity[]) and invertibility with (TimeSeriesInvertibility[]) as shown in the examples.

Can you show a reference or, better yet, give a simple example and the form of the solution you are trying to obtain?

Regards,

Neil

So after trying the code again for some odd reason it's now outputting an error:

Inputting RSolve[{a[n + 1] - 2 a[n] + eps[n] == 1, a[0] == 1}, a[n], n]

Out: RSolve::deqn: Equation or list of equations expected instead of False in the first argument {-2 a[n]+a[1+n]+eps[n]==1,False}.

I've looked at the ARMA reference in the past and its the closest I've gotten to getting my actual specification correctly. Primarily since it includes the stochastic error term correctly. But as far as then providing an functional solution to the associated recurrence equation, I haven't had any luck. Clearly I think there is a function required to transform an ARMA process to it's associated functional form. If I could do that then I imagine I could bracket everything in RSolve[] and find the functional solution.

That aside -- the original code you gave me seemed to work I'm really struggling to see what went wrong.

Best, Ali

Hi Neil,

I'm so sorry I don't know why I wasn't getting notified to a response. Hopefully you're able to see my messages; but I appreciate your help it's still very much an inquiry, as I'm not able to do the calculation in Mathematica. Fortunately I've been using R with no problems. That being said I know once I get it in Mathematica there'll be some really nice functions.

From the example it seems eps is correct thanks for that. But I was really trying to avoid creating seeded random distributions; I was hoping for the explicit analytic solutions. So is there no direct way to transform an ARMA specification to an objective function, to then have solved?

For instance instead of specifying the recurrence equation:

$$ y_t = \alpha_0 + \beta_1 y_{t-1} + \beta_2 y_{t-2} + \epsilon_t + \epsilon_{t-1} $$

Then having it solve it as done above, is there a way to instead specify that equation as an: ARMA{2, 1}? Aside from that obviously inputting the coefficients as well. The use to this would be to test for stationarity, invertibility of the time series.

Thanks again

Ali,

I do not understand from your post what you want to do. I am guessing from your post that you have either a dynamic system with a random input or a random process. You really need to post some example code for anyone to help because it is not clear what "solve" means in your context. For example, are you trying to simulate the process or are you trying to obtain statistics (ie correlations or covariances) or are you trying to use ARMAProcess[] or a related function to create a model? If you post an example, people here can be more helpful.

Regards,

Neil