Gibbs sampler for corrected parametric likelihood + Bernstein-Dirichlet mixture, including possibility of using time series as mere nuisance parameter Multivariate case
gibbs_multivariate_nuisance(data, mcmc_params, corrected, prior_params,
model_params)