pub fn iqr_py(data: &Bound<'_, PyArray1<f64>>) -> PyResult<f64>
Calculate interquartile range (IQR) using optimized partial sort