sliceGammaMean: Generate a random sample of gamma for the R-BSL-M method of...

In ziwenan/BSL: Bayesian Synthetic Likelihood

Description

This function updates the gamma of the R-BSL-M method of \insertCiteFrazier2019;textualBSL with a slice sampler \insertCiteNeal2003BSL. Note this function is mainly designed for internal usage.

Usage

sliceGammaMean(
  ssy,
  ssx,
  loglike,
  gamma = numeric(length(ssy)),
  tau = 1,
  w = 1,
  std = NULL,
  maxit = 1000
)

Arguments

ssy	The observed summary statisic.
ssx	A matrix of the simulated summary statistics. The number of rows is the same as the number of simulations per iteration.
loglike	The current log synthetic likelihood. This is computed with function synLikeMisspec with the current gamma value.
gamma	The additional latent parameter to account for possible incompatability between the model and observed summary statistic. In “BSLmisspec” method, this is updated with a slice sampler \insertCiteNeal2003BSL. The default gamma implies no model misspecification and is equivalent to the standard gaussianSynLike estimator.
tau	Scale (or inverse rate) parameter of the Laplace prior distribution for gamma.
w	Step size for the stepping out in the slice sampler. The default step size is 1.
std	Standard deviation of the columns of ssx. If this is not NULL, it will be used to save computation effort.
maxit	The maximum number of iteration of the stepping out and shrink steps of slice sampler. The default is 1e3.

References

\insertAllCited

See Also

sliceGammaVariance for the slice sampler of the variance inflated target distribution.

ziwenan/BSL documentation built on April 25, 2021, 1:11 p.m.