pub fn correlation(x: &Array1<f64>, y: &Array1<f64>) -> Option<f64>
Compute correlation coefficient (Pearson)