data-raw/sim-crit.R

In kvasilopoulos/exuber: Econometric Analysis of Explosive Time Series

options("exuber.parallel" = TRUE)
options("ncores" = parallel::detectCores())
options("exuber.show_progress" = FALSE)

# Simulation -------------------------------------------------------------

nn <- 600
radf_crit <- vector(list, nn)
for (i in 6:nn) {
  radf_crit[[i]] <-  radf_mc_cv(i, seed = 123, nrep = 2000)
  print(i)
}
names(radf_crit) <- c(paste0("_nan", 1:5), paste0("n", 6:nn))
class(radf_crit) <- c("crit", "list")


# test and store ----------------------------------------------------------

all_equal <- all.equal(radf_crit[[100]], radf_mc_cv(100, nrep = 2000, seed = 123))

if (isTRUE(is_equal)) {
  usethis::use_data(radf_crit, overwrite = TRUE , compress = "xz")
}

waldo::compare(radf_crit[[100]], radf_mc_cv(100, nrep = 2000, seed = 123))
waldo::compare(radf_crit[[100]], exuberdata::radf_crit2[[100]])

# change ------------------------------------------------------------------

for (i in 6:600) {
  class(radf_crit[[i]]) <- c("radf_cv", "mc_cv", "list")
}