pub fn interquartile_range(data: &[f64]) -> Result<f64, String>
Calculate interquartile range (IQR = Q3 - Q1)