powered by
A function for sampling from conditional multivariate normal distributions with mean A^-1b and covariance matrix A^-1.
rmvn_arma(A, b)
A A \(d \times d\) matrix for the Gaussian full conditional distribution precision matrix.
A
matrix
b A \(d\) vector for the Gaussian full conditional distribution mean.
b
vector
# NOT RUN { set.seed(111) A <- diag(4) b <- rnorm(4) sample <- rmvn_arma(A, b) # }
Run the code above in your browser using DataLab