Description Usage Arguments Value References
View source: R/quantileResidualTests.R
get_test_Omega
computes the covariance matrix Omega used in the
quantile residuals tests described by Kalliovirta and Saikkonen 2010.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18  get_test_Omega(
data,
p,
M,
params,
model,
conditional,
parametrization,
constraints,
same_means,
structural_pars = NULL,
g,
dim_g,
ncores = 2,
stat_tol = 0.001,
posdef_tol = 1e08,
df_tol = 1e08
)

data 
a matrix or class 
p 
a positive integer specifying the autoregressive order of the model. 
M 

params 
a real valued vector specifying the parameter values.
Above, φ_{m,0} is the intercept parameter, A_{m,i} denotes the ith coefficient matrix of the mth
mixture component, Ω_{m} denotes the error term covariance matrix of the m:th mixture component, and
α_{m} is the mixing weight parameter. The W and λ_{mi} are structural parameters replacing the
error term covariance matrices (see Virolainen, 2020). If M=1, α_{m} and λ_{mi} are dropped.
If In the GMVAR model, M1=M and ν is dropped from the parameter vector. In the StMVAR model, M1=0.
In the GStMVAR model, the first The notation is similar to the cited literature. 
model 
is "GMVAR", "StMVAR", or "GStMVAR" model considered? In the GStMVAR model, the first 
conditional 
a logical argument specifying whether the conditional or exact loglikelihood function should be used. 
parametrization 

constraints 
a size (Mpd^2 x q) constraint matrix C specifying general linear constraints
to the autoregressive parameters. We consider constraints of form
(φ_{1},...,φ_{M}) = C ψ,
where φ_{m} = (vec(A_{m,1}),...,vec(A_{m,p}) (pd^2 x 1), m=1,...,M,
contains the coefficient matrices and ψ (q x 1) contains the related parameters.
For example, to restrict the ARparameters to be the same for all regimes, set C=
[ 
same_means 
Restrict the mean parameters of some regimes to be the same? Provide a list of numeric vectors
such that each numeric vector contains the regimes that should share the common mean parameters. For instance, if

structural_pars 
If
See Virolainen (2020) for the conditions required to identify the shocks and for the Bmatrix as well (it is W times a timevarying diagonal matrix with positive diagonal entries). 
g 
function g specifying the transformation. 
dim_g 
output dimension of the transformation 
ncores 
the number of CPU cores to be used in numerical differentiation 
stat_tol 
numerical tolerance for stationarity of the AR parameters: if the "bold A" matrix of any regime
has eigenvalues larger that 
posdef_tol 
numerical tolerance for positive definiteness of the error term covariance matrices: if the error term covariance matrix of any regime has eigenvalues smaller than this, the model is classified as not satisfying positive definiteness assumption. Note that if the tolerance is too small, numerical evaluation of the loglikelihood might fail and cause error. 
df_tol 
the parameter vector is considered to be outside the parameter space if all degrees of
freedom parameters are not larger than 
Returns the covariance matrix Omega described by Kalliovirta and Saikkonen 2010.
Kalliovirta L., Meitz M. and Saikkonen P. 2016. Gaussian mixture vector autoregression. Journal of Econometrics, 192, 485498.
Kalliovirta L. and Saikkonen P. 2010. Reliable Residuals for Multivariate Nonlinear Time Series Models. Unpublished Revision of HECER Discussion Paper No. 247.
Virolainen S. 2020. Structural Gaussian mixture vector autoregressive model. Unpublished working paper, available as arXiv:2007.04713.
Virolainen S. 2021. Gaussian and Student's t mixture vector autoregressive model. Unpublished working paper, available as arXiv:2109.13648.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.