“…For the purpose of estimating the covariance matrix of the errors contained in a multivariate time series there can be used any noise reduction algorithm designed for multivariate time series (for instance that proposed in [48] and implemented in the ghkss routine of the TISEAN library [49]). However, as far as we know, the algorithm first proposed in [33] and later improved in [34] is the only one designed for multivariate time series that takes into account either the possible differences between the levels of uncertainty of different coordinates of the time series or the correlations between the coordinates of the error. The performance of the noise reduction algorithm depends strongly on the metric considered, because the metric determines which linear subspace T i is closest to the data points in the neighborhood U i , what is the orthogonal projection of linear models and, for Gaussian errors, these estimates are those of maximum likelihood [30].…”