Demonstrates three RcppArmadillo implementations of sampling random draws from a multivariate normal distribution specified by a precision matrix of special types.
We illustrate the development process of creating code to estimate the parameters of a Probit regression model using the EM algorithm sequentially and in parallel.
Detects how to detect a change point in a time series using Bayesian methodology. Compares naive and vectorized R solutions to one implemented using Rcpp and RcppArmadillo's sample function.