A methodology for Bayesian inference in binary stellar systems based on the No-U-Turn sampler Markov chain Monte Carlo algorithm is presented, providing a precise and efficient estimation of the joint posterior distribution ...