| mix_se-methods | R Documentation | 
Compute standard errors of estimates of MixAR models.
mix_se(x, model, fix_shift)
x | 
 time series.  | 
model | 
 MixAR model, an object inheriting from class “MixAR”.  | 
fix_shift | 
 
  | 
For formulas used in the computation, see \insertCiteWongPhD;textualmixAR.
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.  | 
signature(x = "ANY", model = "list")signature(x = "ANY", model = "MixAR")signature(x = "ANY", model = "MixARGaussian")Davide Ravagli
WongPhDmixAR
## 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'
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.