pub fn quantile_fast(data: &[f64], p: f64) -> Result<f64, String>
O(n) quantile using introselect (R type 7 interpolation).