View source: R/TARMAur.test.B.R
| TARMAur.test.B | R Documentation |
Implements a supremum Lagrange Multiplier unit root test for the null hypothesis of a integrated MA process versus a stationary TARMA process.
TARMAur.test.B(
x,
B = 1000,
pa = 0.25,
pb = 0.75,
thd.range,
method = "ML",
btype = c("wb.r", "wb.n", "iid"),
...
)
x |
A univariate vector or time series. |
B |
Integer. Number of bootstrap resamples. Defaults to 1000. |
pa |
Real number in |
pb |
Real number in |
thd.range |
Vector of optional user defined threshold range. If missing then |
method |
Fitting method to be passed to |
btype |
Bootstrap type, can be one of |
... |
Additional arguments to be passed to |
Implements the bootstrap version of TARMAur.test the supremum Lagrange Multiplier test
to test an integrate MA(1) specification versus a stationary TARMA(1,1) specification.
The option btype specifies the type of bootstrap as follows:
wb.rResidual wild bootstrap with Rademacher auxiliary distribution. See \insertCiteGia22tseriesTARMA.
wb.nResidual wild bootstrap with Normal auxiliary distribution. See \insertCiteGia22tseriesTARMA.
iidResidual iid bootstrap. See \insertCiteGor21btseriesTARMA.
An object of class TARMAtest with components:
statisticThe value of the supLM statistic.
parameterA named vector: threshold is the value that maximises the Lagrange Multiplier values.
test.vVector of values of the LM statistic for each threshold given in thd.range.
thd.rangeRange of values of the threshold.
fit.ARMAThe null model: IMA(1) fit over x.
sigma2Estimated innovation variance from the IMA fit.
data.nameA character string giving the name of the data.
p.valueThe bootstrap p-value of the test.
methodA character string indicating the type of test performed.
dThe delay parameter.
paLower threshold quantile.
TbThe bootstrap null distribution.
Simone Giannerini, simone.giannerini@uniud.it
Greta Goracci, greta.goracci@unibz.it
Cha24tseriesTARMA
TARMAur.test for the asymptotic version of the test. print.TARMAtest for the print method.
## a TARMA(1,1,1,1)
set.seed(123)
x1 <- TARMA.sim(n=100, phi1=c(0.5,-0.5), phi2=c(0.0,0.8), theta1=0.5, theta2=0.5, d=1, thd=0.2)
TARMAur.test.B(x1, B=100) # B=100 for speedup
## a IMA(1,1)
x2 <- arima.sim(n=100, model=list(order = c(0,1,1),ma=0.6))
TARMAur.test.B(x2, B=100) # B=100 for speedup
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.