Function order_statistic 

pub fn order_statistic(values: &[f64], rank: usize) -> f64

Exact 1-based order statistic from an unsorted slice.

This centralizes the sort+select path for code that needs an observed sample value rather than quantile_from_sorted’s linear interpolation.