Optimise beta values for MCMC rungs

Source: R/main.R

repeatedly run burn in phase of MCMC to find vector of beta values that achieve a required coupling acceptance rate between all rungs

optimise_beta(
  proj,
  K = 3,
  target_acceptance = 0.4,
  max_iterations = 100,
  beta_init = seq(0, 1, l = 10),
  silent = FALSE,
  ...
)

Arguments

proj

An rgeoprofile project
K

The value or values of K to optimise beta for
target_acceptance

The minimum acceptance probability to be reached between MCMC chains before the process terminates
max_iterations

Run analyses a maximum number of times before terminating
beta_init

An initial set of beta values to run - a sequence from zero to one
silent

print output from console
...

extra arguments to be sent to the run_mcmc function