Accurate phase diagram calculation from molecular dynamics requires systematic treatment and convergence of statistical averages. In this work we propose a Gaussian process regression based framework for reconstructing the free energy functions using data various origin. Our allows propagating uncertainty finite trajectories to automatically performing with respect simulation parameters. Furthe...