toSSChol: Convert to Non-Innovation State Space Model
Description
This function may not be working properly.
Convert to a non-innovations state space representation using
the given matrix (Om) as the measurement noise covariance.
Om would typically be an estimate of the output noise, such as returned
in $estimates$cov
of the function l
(l.SS
or l.ARMA
).
This assumes that the noise processes in the arbitrary SS representation
are white and uncorrelated.
Usage
toSSChol(model, ...)
# S3 method for TSmodel
toSSChol(model, Om=diag(1,nseriesOutput(model)), ...)
# S3 method for TSestModel
toSSChol(model, Om=NULL, ...)
Arguments
model
An object of class TSmodel.
Om
a matrix to be used as the measurement noise covariance. If Om is
not supplied and model is of class TSestModel then
model$estimates$cov
is used. Otherwise, Om is set to the
identity matrix.
...
arguments to be passed to other methods.
Value
An object of class 'SS' 'TSmodel' containing a state space model which is
not in innovations form.
Details
Convert to a non-innovations SS representation using a Cholesky
decomposition of Om as the coefficient matrix of the output noise.
Examples
Run this code# NOT RUN {
data("eg1.DSE.data.diff", package="dse")
model <- estVARXls(eg1.DSE.data.diff)
model <- toSSChol(model)
# }
Run the code above in your browser using DataLab