CPAT (version 0.1.0)

HR.test: R<U+00E9>nyi-Type Test

Description

Performs the (univariate) R<U+00E9>nyi-type test for change in mean, as described in horvathricemiller19CPAT. This is effectively an interface to stat_Zn; see its documentation for more details. p-values are computed using pZn, which represents the limiting distribution of the test statistic under the null hypothesis, which represents the limiting distribution of the test statistic under the null hypothesis when kn represents a sequence \(t_T\) satisfying \(t_T \to \infty\) and \(t_T/T \to 0\) as \(T \to \infty\). (log and sqrt should be good choices.)

Usage

HR.test(x, kn = log, use_kernel_var = FALSE, stat_plot = FALSE,
  kernel = "ba", bandwidth = "and")

Arguments

x

Data to test for change in mean

kn

A function corresponding to the trimming parameter \(t_T\); by default, the square root function

use_kernel_var

Set to TRUE to use kernel methods for long-run variance estimation (typically used when the data is believed to be correlated); if FALSE, then the long-run variance is estimated using \(\hat{\sigma}^2_{T,t} = T^{-1}\left( \sum_{s = 1}^t \left(X_s - \bar{X}_t\right)^2 + \sum_{s = t + 1}^{T}\left(X_s - \tilde{X}_{T - t}\right)^2\right)\), where \(\bar{X}_t = t^{-1}\sum_{s = 1}^t X_s\) and \(\tilde{X}_{T - t} = (T - t)^{-1} \sum_{s = t + 1}^{T} X_s\); if custom_var is not NULL, this argument is ignored

stat_plot

Whether to create a plot of the values of the statistic at all potential change points

kernel

If character, the identifier of the kernel function as used in cointReg (see getLongRunVar); if function, the kernel function to be used for long-run variance estimation (default is the Bartlett kernel in cointReg)

bandwidth

If character, the identifier for how to compute the bandwidth as defined in cointReg (see getBandwidth); if function, a function to use for computing the bandwidth; if numeric, the bandwidth value to use (the default is to use Andrews' method, as used in cointReg)

Value

A htest-class object containing the results of the test

References

Examples

# NOT RUN {
HR.test(rnorm(1000))
HR.test(rnorm(1000), use_kernel_var = TRUE, kernel = "bo", bandwidth = "nw")
# }