pub fn calculate_mad(observed: &[f64], expected: &[f64]) -> f64
Calculate Mean Absolute Deviation (MAD)