cov.match.sample: Sample with covariance-matching slice sampling

Description Usage Arguments Details Value References See Also

Description

Generate a sample from a probability distribution with the covariance-matching slice sampling method.

Usage

1
2
cov.match.sample(target.dist, x0, sample.size, tuning=1,
                 theta=1, limit=length(x0)*100)

Arguments

target.dist

Target distribution; see make.dist.

x0

Initial coordinates.

sample.size

Sample size to draw.

tuning

A tuning parameter; corresponds to σ_c in sec. 4 of Thompson and Neal (2010).

theta

A factor to scale the crumb standard deviation in every direction after a proposal is rejected. So, after k proposals, crumbs have standard deviation theta^k*tuning in directions orthogonal to all the proposal gradients.

limit

A limit on the number of log-density evaluations per observation before sampling is aborted.

Details

This function implements the covariance-matching method of slice sampling, as described by Thompson and Neal (2010). It can be passed to compare.samplers in the samplers list argument.

Value

A list with elements X, evals, grads, and adapt.rate. adapt.rate indicates the fraction of crumb draws that resulted in adaptation. This sampler follows the calling convention of compare.samplers.

References

Thompson, M. B. and Neal, R. M. (2010). Covariance-adaptive slice sampling. Technical Report TR-1002, Dept. of Statistics, University of Toronto.

See Also

compare.samplers, shrinking.rank.sample



Search within the SamplerCompare package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.