pcts (version 0.15.8)

test_piar: Test for periodic integration

Description

Test if a time series is periodically integrated.

Usage

test_piar(x, d, p, sintercept = FALSE, sslope = FALSE, homoschedastic = FALSE)

Value

a list with the following components:

p

autoregressive order.

spec

values of sintercept, sslope, and homoschedastic, a named logical vector.

statistics

a matrix containing the test statistics (first row) and the corresponding p-values (second row). "LR" is not normalised, so its p-value is NA.

Arguments

x

time series.

d

period.

p

autoregressive order, a positive integer.

sintercept

if TRUE, include seasonal intercept.

sslope

if TRUE, include seasonal slope.

homoschedastic

if TRUE, assume the innovations variance is the same for all seasons.

Author

Georgi N. Boshnakov

Details

Computes test statistics for Franses (1996) test for periodic integration of order 1. The test is based on periodic autoregression of order p, where p can be any positive integer.

References

Boswijk HP and Franses PH (1996). ``Unit roots in periodic autoregressions.'' Journal of Time Series Analysis, 17(3), pp. 221--245.

See Also

pclspiar, pclsdf

Examples

ts1 <- window(dataFranses1996[ , "CanadaUnemployment"],
              start = c(1960, 1), end = c(1987, 4))
test_piar(ts1, 4, 1, sintercept = TRUE)
pcTest(ts1, "piar", 4, 1, sintercept = TRUE) # same

test_piar(ts1, 4, 1, sintercept = TRUE, sslope = TRUE)
test_piar(ts1, 4, 1)              
test_piar(ts1, 4, 1, homoschedastic = TRUE)