adapt_ISMCMC: Independence Sampler/RWM MCMC Adaptation

In JMacDonaldPhD/REpi: Pseudo-Marginal Analysis of partially observed Epidemics

Independence Sampler/RWM MCMC Adaptation

Description

Adapts proposal parameters of INdependence Sampler/RWM MCMC scheme to give near-optimal performance.

Usage

adapt_ISMCMC(
  init,
  U_init,
  epiModel,
  obsFrame,
  epiSample,
  I0,
  alpha,
  logPrior,
  lambda0,
  V0,
  S0,
  noIts,
  delta = 0.05
)

Arguments

init	initial value of epidemic parameters
U_init	initial value of random variables used to construct epidemic trajectory.
epiModel	epidemic model
obsFrame	Generator function for observational model.
epiSample	Observed epidemic data.
I0	Initial state of the epidemic, which assumed to be known.
alpha	Observational parameters to be passed on to the log-likelihood function generated by obsFrame(X_sim)
logPrior	Functions which calculate log-density of the assumed priors of the epidemic parameters.
lambda0	Initial value of the Random Walk Metropolis (RWM) proposal scale parameter
V0	A square matrix containing the initial covariance values for RWM proposal. Shape of matrix should match the length of the init parameter.
S0	Initial number of non-centered random variables to be 'refreshed' in Independence Sampler proposal step.
noIts	Number of iterations of MCMC sampler scheme to carry out.
delta	Probability that the initial proposal parameters are used as opposed to using the current adapted proposal parameters.

Value

A list of the adapted proposal parameters.

JMacDonaldPhD/REpi documentation built on Aug. 2, 2022, 2:09 p.m.