Biocaml 0.4-dev : Biocaml_math.mad

let mad a =

 assert (Array.length a > 0);

 let med = median a in

 let a = Array.map (fun v -> Float.abs (v -. med)) a in

 median a