Bayesian inference of chaotic dynamics by merging data assimilation, machine learning and expectation-maximization (Foundations of data science, 2020)