pub fn box_cox_lambda( series: &Array1<f64>, lambda_range: (f64, f64), steps: usize, ) -> f64
Find optimal Box-Cox lambda using maximum likelihood