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:
statistic
A named vector with the values of the classic supLM
and robust supLMh
statistics.
parameter
A named vector: threshold
is the value that maximises the Lagrange Multiplier values.
test.v
Matrix of values of the LM statistic for each threshold given in thd.range
.
thd.range
Range of values of the threshold.
fit
The null model: AR fit over x
.
sigma2
Estimated innovation variance from the AR fit.
data.name
A character string giving the name of the data.
prop
Proportion of values of the series that fall in the lower regime.
p.value
The p-value of the test. It is NULL
for the asymptotic test.
method
A character string indicating the type of test performed.
d
The delay parameter.
pa
Lower threshold quantile.
dfree
Effective 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.