FARMASIM Call

CALL FARMASIM (series, d <*>, phi <*>, theta <*>, mu <*>, sigma <*>, n <*>, p <*>, q <*>, initial <*>, seed ) ;

The FARMASIM subroutine generates an ARFIMA($p,d,q$) process. The input arguments to the FARMASIM subroutine are as follows:

d

specifies a fractional differencing order. This argument is required; the value of $d$ should be in the open interval $(-1,1)$ excluding zero.

phi

specifies an $m_ p$-dimensional vector that contains the autoregressive coefficients, where $m_ p$ is the number of the elements in the subset of the AR order. The default is zero.

theta

specifies an $m_ q$-dimensional vector that contains the moving average coefficients, where $m_ q$ is the number of the elements in the subset of the MA order. The default is zero.

mu

specifies a mean value. The default is zero.

sigma

specifies a variance of the innovation series. The default is one.

n

specifies the length of the series. The value of $n$ should be greater than or equal to the AR order. The default is $n=100$ is used.

p

specifies the subset of the AR order. See the FARMACOV subroutine for additional details.

q

specifies the subset of the MA order. See the FARMACOV subroutine for additional details.

initial

specifies the initial values of random variables. The initial value is used for the nonstationary process. If initial$=a_0$, then $y_{-p+1},\ldots ,y_{0}$ take the same value $a_0$. If the initial option is not specified, the initial values are set to zero.

seed

is a scalar that contains the random number seed. At the first execution of the subroutine, the seed variable is used as follows:

If seed > 0, the input seed is used for generating the series.

If seed = 0, the system clock is used to generate the seed.

If seed < 0, the value ($-1$)$\times $(seed) is used for generating the series.

If the seed is not supplied, the system clock is used to generate the seed.

On subsequent calls of the subroutine in the DO-loop-like environment, the seed variable is used as follows: If seed > 0, the seed remains unchanged. In other cases, after each execution of the subroutine, the current seed is updated internally.

The FARMASIM subroutine returns the following value:

series

is an $n$ vector that contains the generated ARFIMA($p,d,q$) process.

As an example, consider the following ARFIMA($1,0.3,1$) process:

\[  (1-0.5B)(1-B)^{0.3}(y_ t-10) = (1+0.1B){\epsilon }_ t  \]

In this process, ${\epsilon }_ t \sim NID(0,1.2)$. The following statements generate this process:

d = 0.3;
phi = 0.5;
theta = -0.1;
mu    = 10;
sigma = 1.2;
call farmasim(yt, d, phi, theta, mu, sigma, 10) seed=1234;
print yt;

Figure 23.116: Data Simulated from a ARFIMA Process

yt
12.17358
13.954495
15.817231
15.94882
12.25926
13.641022
13.399623
11.930759
10.049435
9.1445036


The FARMASIM subroutine generates a time series of length $n$ from an ARFIMA($p,d,q$) model. If the process is stationary and invertible, the initial values $y_{-p+1}, \ldots , y_{0}$ are produced by using covariance matrices obtained from FARMACOV. If the process is nonstationary, the time series is recursively generated by using the user-defined initial value or the zero initial value.

To generate an ARFIMA($p,d,q$) process with $d\in [0.5,1)$, $x_ t$ is first generated for $d’\in (-0.5, 0)$, where $d’=d-1$ and then $y_ t$ is generated by $y_ t = y_{t-1} + x_ t$.

To generate an ARFIMA($p,d,q$) process with $d\in (-1,-0.5]$, a two-step approximation based on a truncation of the expansion $(1-B)^ d$ is used; the first step is to generate an ARFIMA($0,d,0$) process $x_ t=(1-B)^{-d}{\epsilon }_ t$, with truncated moving average weights; the second step is to generate $y_ t =\phi (B)^{-1}\theta (B)x_ t$.