Simple estimation using ABC-SMC
Create a reference trajectory
sum_stat_obs <- c(2.0, 0.75)
Run abc smc procedure
res <- abcsmc(model_list = model_list,
prior_dist = prior_dist,
ss_obs = sum_stat_obs,
max_number_of_gen = 15,
nb_acc_prtcl_per_gen = 2000,
new_threshold_quantile = 0.8,
experiment_folderpath = "smpl",
max_concurrent_jobs = 5,
verbose = FALSE)
Plot results
all_accepted_particles <- res$particles
all_thresholds <- res$thresholds
plot_abcsmc_res(data = all_accepted_particles, prior = prior_dist,
filename = "smpl/res/figs/smpl_pairplot_all.png", colorpal = "YlGnBu")
#> [1] "Plot saved as 'png'."
#> Registered S3 method overwritten by 'GGally':
#> method from
#> +.gg ggplot2
plot_thresholds(data = all_thresholds, nb_threshold = 1,
filename = "smpl/res/figs/smpl_thresholds.png", colorpal = "YlGnBu")
#> [1] "Plot saved as 'png'."
plot_ess(data = all_accepted_particles,filename = "smpl/res/figs/smpl_ess.png", colorpal = "YlGnBu")
#> [1] "Plot saved as 'png'."
#> gen ess
#> 1 1 2000.000
#> 2 2 1881.592
#> 3 3 1905.910
#> 4 4 1930.612
#> 5 5 1924.049
#> 6 6 1934.766
#> 7 7 1932.013
#> 8 8 1913.087
#> 9 9 1921.480
#> 10 10 1911.549
#> 11 11 1911.130
#> 12 12 1891.114
#> 13 13 1880.753
#> 14 14 1863.416
#> 15 15 1873.565
plot_densityridges(data = all_accepted_particles, prior = prior_dist,
filename = "smpl/res/figs/smpl_densityridges.png", colorpal = "YlGnBu")
#> [1] "Plot saved as 'png'."

Pairplot of all iterations

Threshold evolution over iterations

ESS evolution over iterations

Density estimates for alpha

Density estimates for beta