T_stationary | R Documentation |
A hypothesis test for stationarity of functional time series.
T_stationary(sample, L = 49, J = 500, MC_rep = 1000, cumulative_var = .90,
Ker1 = FALSE, Ker2 = TRUE, h = ncol(sample)^.5, pivotal = FALSE,
use_table = FALSE, significance)
sample |
A matrix of discretised curves of dimension (p by n), where p represents the dimensionality and n represents sample size. |
L |
Number of Fourier basis functions. |
J |
Truncation level used to approximate the distribution of the squared integrals of Brownian bridges that appear in the limit distribution. |
MC_rep |
Number of replications. |
cumulative_var |
Amount of variance explained. |
Ker1 |
Flat top kernel in (4.1) of Horvath et al. (2014). |
Ker2 |
Flat top kernel in (7) of Politis (2003). |
h |
Kernel bandwidth. |
pivotal |
If |
use_table |
If |
significance |
Level of significance. Possibilities are “10%”, “5%”, “1%”. |
As in traditional (scalar and vector) time series analysis, many inferential procedures for functional time series assume stationarity. Stationarity is required for functional dynamic regression models, for bootstrap and resampling methods for functional time series and for the functional analysis of volatility.
p-value |
When |
Greg. Rice and Han Lin Shang
L. Horvath and Kokoszka, P. (2012) Inference for Functional Data with Applications, Springer, New York.
L. Horvath, P. Kokoszka, G. Rice (2014) "Testing stationarity of functional time series", Journal of Econometrics, 179(1), 66-82.
D. N. Politis (2003) "Adaptive bandwidth choice", Journal of Nonparametric Statistics, 15(4-5), 517-533.
A. Aue, G. Rice, O. S\"onmez (2018) "Detecting and dating structural breaks in functional data without dimension reduction", Journal of the Royal Statistical Society: Series B, 80(3), 509-529.
farforecast
result = T_stationary(sample = pm_10_GR_sqrt$y)
result_pivotal = T_stationary(sample = pm_10_GR_sqrt$y, J = 100, MC_rep = 5000,
h = 20, pivotal = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.