C++ implementation of Zellner and Ando (2010) Direct Monte Carlo method for sampling from the posterior of a Bayesian SUR
sur_sample_cpp(Y, Xlist, y, X, XtX, pj, M)matrix \((y_1, \ldots y_J)\)
A list, each element a design matrix
vector of responses
design matrix
matrix giving crossprod(cbind(X1, ..., XJ))
vector giving number of covariates per endpoint
An integer giving the number of desired samples