Gibbs sampler for corrected parametric likelihood + Bernstein-Dirichlet mixture, including possibility of using time series as mere nuisance parameter
gibbs_nuisance( data, mcmc_params, corrected, prior_params, model_params, full_lik = NULL )
Useful links