mix_se-methods {mixAR}

Compute standard errors of estimates of MixAR models

Description

Compute standard errors of estimates of MixAR models.

Usage

mix_se(x, model, fix_shift)

Arguments

x

time series.

model

MixAR model, an object inheriting from class “MixAR”.

fix_shift

logical. Should the shift paramters be fixed? Default is FALSE.

Details

For formulas used in the computation, see Wong (1998).

Value

a list with components:

standard_errors

Standard error of parameter estimates,

covariance_matrix

The covariance matrix, obtained as inverse of the information matrix,

Complete_Information

Complete information matrix,

Missing_Information

Missing information matrix.

Methods

signature(x = "ANY", model = "list")

signature(x = "ANY", model = "MixAR")

signature(x = "ANY", model = "MixARGaussian")

Author(s)

Davide Ravagli

References

Wong CS (1998). Statistical inference for some nonlinear time series models. Ph.D. thesis, University of Hong Kong, Hong Kong .

Examples

## Example with IBM data

## data(ibmclose, package = "fma")

moWLprob <- exampleModels$WL_ibm@prob    # 2019-12-15; was: c(0.5339,0.4176,0.0385)     
moWLsigma <- exampleModels$WL_ibm@scale  #                  c(4.8227,6.0082,18.1716)
moWLar <- list(-0.3208, 0.6711,0)        # @Davide - is this from some model?

moWLibm <- new("MixARGaussian", prob = moWLprob, scale = moWLsigma, arcoef = moWLar)

IBM <- diff(fma::ibmclose)
mix_se(as.numeric(IBM), moWLibm, fix_shift = TRUE)$'standard_errors'