Bayesian parameter estimation through variational methods