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