test_piar: Test for periodic integration

In GeoBosh/pcts: Periodically Correlated and Periodically Integrated Time Series

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)

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.

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.

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.

Note

Currently only the case p = 1 is handled, for p > 1 the statistics are set to NA. :TODO: handle this.

All statistics are computed but some p-values are not computed yet.

Author(s)

Georgi N. Boshnakov

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)              

GeoBosh/pcts documentation built on Dec. 8, 2023, 9:57 p.m.