Monte Carlo Critical Values

radf_mc_cv computes Monte Carlo critical values for the recursive unit root tests. radf_mc_distr computes the distribution.

radf_mc_cv(n, minw = NULL, nrep = 1000L, seed = NULL)

radf_mc_distr(n, minw = NULL, nrep = 1000L, seed = NULL)

Arguments

n

A positive integer. The sample size.

minw

A positive integer. The minimum window size (default = \((0.01 + 1.8/\sqrt(T))T\), where T denotes the sample size).

nrep

A positive integer. The number of Monte Carlo simulations.

seed

An object specifying if and how the random number generator (rng) should be initialized. Either NULL or an integer will be used in a call to set.seed before simulation. If set, the value is saved as "seed" attribute of the returned value. The default, NULL, will not change rng state, and return .Random.seed as the "seed" attribute. Results are different between the parallel and non-parallel option, even if they have the same seed.

Value

For radf_mc_cv a list that contains the critical values for ADF, BADF, BSADF and GSADF test statistics. For radf_mc_distr a list that contains the ADF, SADF and GSADF distributions.

See also

radf_wb_cv for wild bootstrap critical values and radf_sb_cv for sieve bootstrap critical values

Examples

# \donttest{
# Default minimum window
mc <- radf_mc_cv(n = 100)

tidy(mc)
#> # A tibble: 3 × 4
#>   sig       adf  sadf gsadf
#>   <fct>   <dbl> <dbl> <dbl>
#> 1 90    -0.451   1.03  1.64
#> 2 95    -0.0189  1.31  1.92
#> 3 99     0.616   1.81  2.58

# Change the minimum window and the number of simulations
mc2 <- radf_mc_cv(n = 100, nrep = 600, minw = 20)

tidy(mc2)
#> # A tibble: 3 × 4
#>   sig       adf  sadf gsadf
#>   <fct>   <dbl> <dbl> <dbl>
#> 1 90    -0.422  0.938  1.58
#> 2 95    -0.0849 1.28   1.79
#> 3 99     0.410  1.85   2.46

mdist <- radf_mc_distr(n = 100, nrep = 1000)

autoplot(mdist)

# }