(ii) We are able to calculate non-Gaussian errors on the parameter thus taking into account the full shape of the posterior distribution. This is a feature which was studied for example in (but without including errors in the independent variable) but is ignored in the literature otherwise. This can lead to biased error estimates depending on the data (see appendix A for more details).