This paper has two interrelated foci: (i) obtaining stable and efficient
data-driven closure models by using a multivariate time series of partial
observations from a large-dimensional system; and (ii) comparing these closure
models with the optimal closures predicted by the Mori-Zwanz