uGMAR (version 3.1.0)

get_regime_means: Calculate regime specific means ${\mu }_m$

Description

get_regime_means calculates regime means ${\mu }_m={\varphi }_{m,0}/(1-\sum {\varphi }_{i,m})$ of the given GMAR, StMAR or G-StMAR model

Usage

get_regime_means(gsmar)

Arguments

gsmar

object of class 'gsmar' created with the function fitGSMAR or GSMAR.

Value

Returns a length M vector containing regime mean ${\mu }_m$ in the m:th column, $m=1,..,M$.

References

• Kalliovirta L., Meitz M. and Saikkonen P. 2015. Gaussian Mixture Autoregressive model for univariate time series. Journal of Time Series Analysis, 36, 247-266.

• Meitz M., Preve D., Saikkonen P. 2018. A mixture autoregressive model based on Student's t-distribution. arXiv:1805.04010 [econ.EM].

• There are currently no published references for the G-StMAR model, but it's a straightforward generalization with theoretical properties similar to the GMAR and StMAR models.

See Also

condMoments, uncondMoments, get_regime_vars, get_regime_autocovs

Other moment functions: condMoments, get_regime_autocovs, get_regime_vars, uncondMoments

Examples

# NOT RUN {
# GMAR model
params13 <- c(1.4, 0.88, 0.26, 2.46, 0.82, 0.74, 5.0, 0.68, 5.2, 0.72, 0.2)
gmar13 <- GSMAR(p=1, M=3, params=params13, model="GMAR")
get_regime_means(gmar13)

# StMAR model
params12t <- c(1.38, 0.88, 0.27, 3.8, 0.74, 3.15, 0.8, 100, 3.6)
stmar12t <- GSMAR(p=1, M=2, params=params12t, model="StMAR")
get_regime_means(stmar12t)

# G-StMAR model (similar to the StMAR model above)
params12gs <- c(1.38, 0.88, 0.27, 3.8, 0.74, 3.15, 0.8, 3.6)
gstmar12 <- GSMAR(p=1, M=c(1, 1), params=params12gs, model="G-StMAR")
get_regime_means(gstmar12)
# }