Changes in plot function: sampling vs numerical approximation
Added **likelihood_params to predictive_values