# Load ABC input data from example simulation
abc_input <- readRDS(
system.file("extdata", "rdm_minimal", "abc", "abc_input.rds", package = "eam")
)
# Perform ABC resampling
results <- abc_resample(
target = abc_input$target,
param = abc_input$param,
sumstat = abc_input$sumstat,
n_iterations = 100,
n_samples = 100,
tol = 0.5,
method = "rejection"
)
# plot forest plots showing parameter ranges
plot_resample_forest(results, ci_level = 0.95)
Run the code above in your browser using DataLab