Bootstrap resample ABC posterior samples

Usage

abc_posterior_bootstrap(abc_result, n_samples, replace = TRUE)

Arguments

abc_result: An abc object from abc
n_samples: Number of bootstrap samples to draw (default 1000)
replace: Logical, whether to sample with replacement (default TRUE)

Value

Data frame of bootstrapped parameter values

Examples

# Load an example abc output, you should generate it by applying ABC to your data
# check abc_abc for details on fitting ABC models
rdm_minimal_example <- system.file("extdata", "rdm_minimal", package = "eam")
abc_model <- readRDS(file.path(rdm_minimal_example, "abc", "abc_neuralnet_model.rds"))

# Bootstrap resample posterior parameters
posterior_params <- abc_posterior_bootstrap(
  abc_model,
  n_samples = 100
)

# View the first few rows of the bootstrapped posterior parameters
head(posterior_params)
#>        V_beta_1 V_beta_group
#> X150  0.2258720    0.1938078
#> X186  0.2188160    0.1983776
#> X61   0.2587934    0.1670919
#> X61.1 0.2587934    0.1670919
#> X470  0.2612701    0.1621439
#> X318  0.2608088    0.1624729