pub fn population_mean_median( theta: &Array2<f64>, w: &Array1<f64>, ) -> Result<(Array1<f64>, Array1<f64>)>