pub fn wskew( w: &[f64], wstride: usize, data: &[f64], stride: usize, n: usize) -> f64
This function computes the weighted skewness of the dataset data.
skew = (\sum w_i ((x_i - \Hat x)/\Hat \sigma)^3) / (\sum w_i)