Description Usage Arguments Value Author(s) References See Also Examples
View source: R/lineqGPsamplers.R
Sampler for truncated multivariate normal distributions
via Hamiltonian Monte Carlo using the package tmg
(Pakman and Paninski, 2014).
1 2 |
object |
an object with |
nsim |
an integer corresponding to the number of simulations. |
control |
extra parameters required for the MC/MCMC sampler. |
... |
further arguments passed to or from other methods. |
A matrix with the simulated samples. Samples are indexed by columns.
A. F. Lopez-Lopera.
Pakman, A. and Paninski, L. (2014), "Exact Hamiltonian Monte Carlo for truncated multivariate Gaussians". Journal of Computational and Graphical Statistics, 23(2):518-542. [link]
1 2 3 4 5 6 7 8 9 | n <- 100
x <- seq(0, 1, length = n)
Sigma <- kernCompute(x1 = x, type = "gaussian", par = c(1,0.2))
tmgPar <- list(mu = rep(0,n), Sigma = Sigma + 1e-9*diag(n), lb = rep(-1,n), ub = rep(1,n))
class(tmgPar) <- "HMC"
y <- tmvrnorm(tmgPar, nsim = 10)
matplot(x, y, type = 'l', ylim = c(-1,1),
main = "Constrained samples using Hamiltonian MC")
abline(h = c(-1,1), lty = 2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.