pub fn chi_squared_test(observed: &[f64], expected: &[f64]) -> (f64, f64)
Chi-squared goodness-of-fit test. Returns (chi2-statistic, p-value).