Rao_test: Perform Rao's score test for a GSMVAR model

In saviviro/gmvarkit: Estimate Gaussian and Student's t Mixture Vector Autoregressive Models

Perform Rao's score test for a GSMVAR model

Description

Rao_test performs Rao's score test for a GSMVAR model

Usage

Rao_test(gsmvar)

Arguments

gsmvar

an object of class 'gsmvar' generated by fitGSMVAR or GSMVAR, containing the model specified by the null hypothesis (i.e., the constrained model).

Details

Tests the constraints imposed in the model given in the argument GSMVAR. This implementation uses the outer product of gradients approximation in the test statistic.

Value

A list with class "hypotest" containing the test results and arguments used to calculate the test.

References

• Buse A. (1982). The Likelihood Ratio, Wald, and Lagrange Multiplier Tests: An Expository Note. The American Statistician, 36(3a), 153-157.

See Also

LR_test, Wald_test, fitGSMVAR, GSMVAR, diagnostic_plot, profile_logliks

Examples

## These are long running examples that use parallel computing!
## The below examples take around 30 seconds to run.

# Structural GMVAR(2, 2), d=2 model with recursive identification
# with the AR matrices  restricted to be the identical across the regimes:
W22 <- matrix(c(1, NA, 0, 1), nrow=2, byrow=FALSE)
C_mat <- rbind(diag(2*2^2), diag(2*2^2))
fit22sc <- fitGSMVAR(gdpdef, p=2, M=2, constraints=C_mat,
                     structural_pars=list(W=W22), ncalls=1, seeds=1)

# Test the null:
Rao_test(fit22sc)

saviviro/gmvarkit documentation built on March 8, 2024, 4:15 a.m.