- Compute 100 $(\rho^2, \kappa, \gamma$) points for the computer
experiment and check whether MLE asymptotically exists for each
one, and produce Figure 1 of the manuscript.
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary/";
ncores <- 10;
npoints <- 100;
seed <- 0;
source(file.path(project_path, "code/0-design.R"))'
- Training phase
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary/";
ncores <- 10;
nobs <- 2000;
psi <- 0;
beta_star_setting <- "u2";
repetitions <- 100;
seed <- 101;
source(file.path(project_path, "code/1-train-set.R"))'
- Predict conjectured approximation and produce Figure 2 and Table 1
of the manuscript.
Rscript --no-init-file -e \
'project_path <- "~/Repositories//mJPL-conjecture-supplementary";
nobs <- 2000;
psi <- 0;
beta_star_setting <- "u2";
source(file.path(project_path, "code/2-conjecture.R"))'
- Produce test sets
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 1000;
beta_star_setting <- "s1";
seed <- 103;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 1000;
beta_star_setting <- "s2";
seed <- 104;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 1000;
beta_star_setting <- "u1";
seed <- 105;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 1000;
beta_star_setting <- "u2";
seed <- 106;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 2000;
beta_star_setting <- "s1";
seed <- 107;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 2000;
beta_star_setting <- "s2";
seed <- 108;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 2000;
beta_star_setting <- "u1";
seed <- 109;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 2000;
beta_star_setting <- "u2";
seed <- 110;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 3000;
beta_star_setting <- "s1";
seed <- 111;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 3000;
beta_star_setting <- "s2";
seed <- 112;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 3000;
beta_star_setting <- "u1";
seed <- 113;
source(file.path(project_path, "code/3-test-set.R"))'
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
nobs <- 3000;
beta_star_setting <- "u2";
seed <- 114;
source(file.path(project_path, "code/3-test-set.R"))'
- Compute phase transition curves for the setting where the test sets
have been produced
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ncores <- 10;
seed <- 123
source(file.path(project_path, "code/4-test-set-pt.R"))'
- Reproduce Figure 3 of the manuscript
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary";
ns <- c(1000, 2000, 3000)
beta_star <- c("s1", "s2", "u1", "u2")
conjecture_model <- "conjecture-n-2000-beta-u2-psi-0.rda"
source(file.path(project_path, "code/5-test-summaries.R"))'
- Reproduce Figure 4 and Figure 5 in manuscript
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary"
nobs <- c(3000, 1000)
beta_star_setting <- c("u1", "s1")
psi <- c(0.3, 0.6)
rhosq <- c(0.6, 0.1)
conjecture_model <- "conjecture-n-2000-beta-u2-psi-0.rda"
source(file.path(project_path, "code/6-test-inset-plots.R"))'
- Compute reported summaries on elapsed time in Table 2
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary"
nobs <- 2000
beta_star_setting <- "s2"
source(file.path(project_path, "code/7-computational-performance.R"))'
- Reproduce Figure 6 in manuscript
Rscript --no-init-file -e \
'project_path <- "~/Repositories/mJPL-conjecture-supplementary"
nobs <- 2000
ncores <- 10
beta_star_setting <- "s1"
seed <- 111
rsqs <- 0.3
conjecture_model <- "conjecture-n-2000-beta-u2-psi-0.rda"
source(file.path(project_path, "code/8-bernoulli.R"))'