Bayesian inference of chaotic dynamics by merging data assimilation, machine learning and expectation-maximization