powered by
Gibbs sampling algorithm for VAR model
gibbs_VAR_nuisance_intern( data, mcmc_params, prior_params, model_params, full_lik = F )