Function that calculates the empirical cross-covariance of order h for a bivariate time series.
sigma_zf(h, vec1, vec2, N)
the lag value.
observations for the first variable of the bivariate time series.
observations for the second variable of the bivariate time series.
the common length of vec1 and vec2.
The value of lambda that corresponds to a smoothing level s.