Function partial_correlation 

pub fn partial_correlation(
    data: &[Vec<f64>],
    x: usize,
    y: usize,
    z: &[usize],
) -> f64

Compute the partial correlation of X and Y given the set Z.

Uses recursive formula via the Gram matrix.