pub fn cor(x: &[f64], y: &[f64]) -> Result<f64, String>
Pearson correlation coefficient between two arrays. cor(x, y) = cov(x,y) / (sd(x) * sd(y))