Displaying 1 result from an estimated 1 matches for "xsim1".
Did you mean:
sim1
2009 Jul 21
0
Specifying initial values for arima.sim
...it
> xarima = arima(xs,order=c(1,0,0))
> xarima
Call:
arima(x = xs, order = c(1, 0, 0))
Coefficients:
ar1 intercept
0.9895 8.6341
s.e. 0.0106 6.1869
I should then be able to simulate this ARIMA process, using the
residuals. Lets do this twice for comparison
> xsim1 <-arima.sim(n = 100,innov=residuals(xarima),list(ar = c(0.9895)), )
> xsim2 <-arima.sim(n = 100,innov=residuals(xarima),list(ar = c(0.9895)), )
> xsim1[1]
[1] -4.855137
> xsim2[1]
[1] 5.511827
> xs[1]
[1] 1.014863
Clearly these series are starting from different initial values....