Library: | times |
See also: | genarch |
Quantlet: | genarma | |
Description: | generates an autoregressive moving average process (ARMA) in deviations from the mean form y_t = a(B)y_t + eps_t + b(B)eps_t B denotes the backshift (or lag) operator. You have to deliver the AR coefficients vector a, the MA coefficients vector b and the (T x 1) white noise series eps. The output y is a (T x 1) ARMA series. The simulation routine sets the necessary initial values to 0 (that is y_t=eps_t=0 for t<1). |
Usage: | y = genarma(a,b,eps) | |
Input: | ||
a | (p x 1) vector of AR coefficients | |
b | (q x 1) vector of MA coefficients | |
eps | (T x 1) vector of white noise disturbances | |
Output: | ||
y | (T x 1) vector with the generated ARMA process |
library("times") ; loads the quantlets from times library randomize(0) ; sets the seed y = genarma(0.5,0.5,normal(1000)) pacfplot(y,12) ; display of the sample pacf for y_t acfplot(y,12) ; display of the sample acf for y_t
Two displays that show the acf and the pacf of the generated ARMA(1,1) process y_t. The series y has T=1000 observations. The displays resemble the exhibits in Mills (cited above), p. 89.