Towards optimal dynamical models of the Milky Way halo