pub fn autocorrelation(data: &[f64], lag: usize) -> f64
Compute the empirical autocorrelation at lag k.