| TAR.test | R Documentation |
Heteroskedasticity robust supremum Lagrange Multiplier test for a AR specification versus a TAR specification. Includes the classic (non robust) AR versus TAR test.
TAR.test(x, pa = 0.25, pb = 0.75, ar.ord, d = 1)
x |
A univariate time series. |
pa |
Real number in |
pb |
Real number in |
ar.ord |
Order of the AR part. |
d |
Delay parameter. Defaults to |
Implements a heteroskedasticity robust asymptotic supremum Lagrange Multiplier test to test an AR specification versus a TAR specification.
This is an asymptotic test and the value of the test statistic has to be compared with the critical values tabulated
in \insertCiteGor21tseriesTARMA or \insertCiteAnd03tseriesTARMA.
Both the non-robust supLM and the robust supLMh statistics are returned.
An object of class TARMAtest with components:
statisticA named vector with the values of the classic supLM and robust supLMh statistics.
parameterA named vector: threshold is the value that maximises the Lagrange Multiplier values.
test.vMatrix of values of the LM statistic for each threshold given in thd.range.
thd.rangeRange of values of the threshold.
fitThe null model: AR fit over x.
sigma2Estimated innovation variance from the AR fit.
data.nameA character string giving the name of the data.
propProportion of values of the series that fall in the lower regime.
p.valueThe p-value of the test. It is NULL for the asymptotic test.
methodA character string indicating the type of test performed.
dThe delay parameter.
paLower threshold quantile.
dfreeEffective degrees of freedom. It is the number of tested parameters.
Simone Giannerini, simone.giannerini@uniud.it
Greta Goracci, greta.goracci@unibz.it
Gor23tseriesTARMA
And03tseriesTARMA
TAR.test.B for the bootstrap version of the test.
TARMA.test for the (robust) ARMA vs TARMA asymptotic version of the test, which includes also the AR vs TAR test, with different defaults.
TARMAGARCH.test for the robust version of the ARMA vs TARMA test that assumes GARCH innovations.
TARMA.sim to simulate from a TARMA process.
set.seed(123)
## a TAR(1,1) ---------------
x1 <- TARMA.sim(n=100, phi1=c(0.5,-0.5), phi2=c(0.0,0.8), theta1=0, theta2=0, d=1, thd=0.2)
TAR.test(x1, ar.ord=1, d=1)
## a AR(1) ----------------
x2 <- arima.sim(n=100, model=list(order=c(1,0,0), ar=0.5))
TAR.test(x2, ar.ord=1, d=1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.