Locally fit multivariate autoregressive models to non-stationary time series by a Bayesian procedure.
blomar(y, max.order = NULL, span)
mean.
variance.
Bayesian weight.
AIC with respect to the present data.
AR coefficients. arcoef[[m]][i,j,k]
shows the value of
\(i\)-th row, \(j\)-th column, \(k\)-th order of \(m\)-th model.
innovation variance.
equivalent AIC of Bayesian model.
start point of the data fitted to the current model.
end point of the data fitted to the current model.
A multivariate time series.
upper limit of the order of AR model, less than or equal to
\(n/2d\) where \(n\) is the length and \(d\) is the dimension of the
time series y
.
Default is \(min(2 \sqrt{n}, n/2d)\).
length of basic local span. Let \(m\) denote max.order
,
if \(n-m-1\) is less than or equal to span
or \(n-m-1-\)span
is less than \(2md\), span
is \(n-m\).
The basic AR model is given by
$$y(t) = A(1)y(t-1) + A(2)y(t-2) + \ldots + A(p)y(t-p) + u(t),$$
where \(p\) is order of the AR model and \(u(t)\) is innovation variance
v
.
G.Kitagawa and H.Akaike (1978) A Procedure for the Modeling of Non-stationary Time Series. Ann. Inst. Statist. Math., 30, B, 351--363.
H.Akaike (1978) A Bayesian Extension of The Minimum AIC Procedure of Autoregressive Model Fitting. Research Memo. NO.126. The institute of Statistical Mathematics.
H.Akaike, G.Kitagawa, E.Arahata and F.Tada (1979) Computer Science Monograph, No.11, Timsac78. The Institute of Statistical Mathematics.
data(Amerikamaru)
blomar(Amerikamaru, max.order = 10, span = 300)
Run the code above in your browser using DataLab